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Introduction 


During  the  period  of  August  1981  to  May  1982  DPSI  was 
under  contract  (F19628-81-C-0141)  to  the  Cloud  Physics  Branch 
(LYC) ,  Meteorology  Division  of  the  Air  Force  Geophysics  Lab¬ 
oratory  (AFGL) .  LYC  is  the  Air  Force  office  abbreviation 
for  Cloud  Physics  Branch  and  will  be  used  interchangeably 
throughout  this  report.  The  purpose  of  the  contract  was  to 
develop  and  apply  mathematical  procedures  to  a  variety  of 
standard  and  non-standard  cloud  physics  research  data. 

Work  t  under  this  latest  contract  t  was  performed  on  two 
distinct  computer  systems.  The  AFGL  inhouse  Cyber  74  Com¬ 
puter  Systems  (Control  Data  Corp.)  and  the  DEC  (Digital  Equip¬ 
ment  Corp. )  PDP-8/E.  This  report  is  prepared  with  this  dis¬ 
tinction  in  mind.  Chapter  2  describes  CDC  Cyber  74  programming 
while  chapters  3  and  4  detail  real-time  programming  for  the 
PDP-8/E  installed  on  an  AFGL  operated  MC-130E  instrumented 
aircraft-  and  inhouse  LYC  PDP-8/E  for  testing  of  flight  pro¬ 
gramming  and  post  flight  analysis. 

The  work  performed  under  this  contract  has  been  sub¬ 
mitted,  in  detail,  in  a  set  of  22  monthly  reports.  This 
document  will  summarize  changes  to  existing  programs;  in 
addition,  new  programs  and  procedures  will  be  fully  described 
in  this  document.  A  complete  set  of  updated  documentation 
and  operating  instructions  for  each  program  is  included. 
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AFGL  CYBER  74  Computer  System 

The  AFGL  computation  center  consists  of  two  Control  Data 
Corporation  6600  computers.  The  two  machines  are  referred  to 
as  system  I  and  system  II.  Each  computer  is  composed  of  cen¬ 
tral  memory,  a  high  speed  central  processor,  14  peripheral 
processors,  and  a  display  console.  Up  to  15  user  jobs  in  each 
system  can  be  in  central  memory  simultaneously,  sharing  the 
central  processor  at  scheduled  intervals.  Each  job  in  central 
memory  gains  access  to  the  central  processor  sequentially  un¬ 
til  execution  is  complete.  Each  computer  is  under  complete 
control  of  the  NOS/BE  operating  system. 

Both  systems  share  a  non-volatile  disk  package  for  per¬ 
manent  file  storage.  A  permanent  file  is  a  mass  storage  file 
cataloged  by  either  system,  so  that  its  location  and  identifi¬ 
cation  are  known  to  each  system  always.  Frequently  used  pro¬ 
grams,  subroutines  and  data  bases  are  immediately  available 
to  jobs  requesting  them  without  operator  intervention.  Perm¬ 
anent  files  cannot  be  destroyed  accidentally  during  normal 
system  operation,  including  normal  dead  start,  because  they 
are  protected  by  the  system  from  unauthorized  access  according 
to  the  privacy  controls  specified  when  they  are  created. 

The  plotting  facilities,  also  available  to  both  systems, 
consist  of  a  single  off-line  mode  of  operation.  This  mode 
produces  plots  on  one  of  two  media;  paper  or  film.  A  more 
detailed  description  of  these  facilities  appear  in  later 
paragraphs  of  this  document.  Batch  processing  is,  of  course, 
supported  by  both  systems.  This  allows  the  user  access  to 


AFGL  CYBER  72  Computer  System  (cont'd) 

any  peripheral  device  (tape  drives,  disks,  plotters,  etc.) 
that  may  be  necessary  for  program  execution. 

System  I,  the  more  heavily  used  machine,  supports  the 
time  sharing  system,  "INTERCOM",  as  well.  Time  sharing  ter¬ 
minals  are  available  to  LYC  using  Bell  system  terminals  via 
dedicated  telephone  lines  to  the  6600  at  300  baud  rate.  In 
addition,  there  are  dial-up  lines  available  for  use  v/ith 
portable  terminals.  Dial-up  terminals  must  be  connected 
through  an  acoustic  coupler  over  standard,  voice-grade  tele¬ 
phone  lines. 

INTERCOM  allows  the  user  to  create,  modify,  or  execute 
a  program  in  real  time.  In  addition  it  allows  for  data  base 
modification  such  as  permanent  file  creation  or  deletion. 

There  are  two  modes  of  operation  using  INTERCOM:  COMMAND 
mode  and  EDIT  mode.  The  EDIT  mode  is  more  often  used,  since 
the  COMMAND  mode  is  restrictive  and  anything  performed  in  the 
COMMAND  mode  can  be  done  in  EDIT  mode.  For  a  complete  list 
of  INTERCOM  commands  with  their  description,  the  reader  is 
referred  to  the  CDC  publication  "6000  Series  INTERCOM  v5 
Reference  Manual". 

As  with  all  programming  systems,  INTERCOM  will  require 
a  great  deal  of  experience  in  order  to  gain  expertise.  There 
are  many  "shortcuts"  that  may  be  employed,  and  a  rather  long 
list  of  extended  capabilities. 
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AFGL  CYBER  74  Computer  System  (cont'd) 

The  plotting  facilities  available  at  the  Computation 
Center  are  capable  of  producing  either  paper  or  film  plots. 
These  facilities  offer  a  wide  choice  for  data  display.  The 
pen  and  ink  plotters  will  produce  plots  of  large  dimensions, 
excellent  line  definition  and  multiple  colors.  The  film 
plotters  offer  high  speed  plots,  large  amount  of  shading  and 
line  flexibility. 

The  off-line  hardware  consists  of  a  Calcomp  780  Magnetic 
Tape  Unit  with  both  a  Drum-Type  Digital  Incremental  Zip  Mode 
30  Inch  Dual  Step  Size  Plotter  and  a  CRT  plotter.  Standard 
paper  is  available  in  either  12  or  30  inches.  Ballpoint 
colors  are  black,  red  or  blue,  in  addition  to  red  liquid  ink 
for  photo  reproduction. 

There  are  two  film  systems  available;  35mm  microfilm  and 
105mm  microfiche.  We  have  found  the  microfiche  system  to  be 
most  beneficial.  Programs  like  RAP*  produce  a  multitude  of 
plots  which  can  easily  be  plotted  on  only  a  few  sheet  or  film 
This  is  very  convenient  for  the  LYC  scientists.  They  can 
study  the  plots  at  a  viewer  and  make  hard  copies  of  just  the 
ones  they  find  most  useful.  This  eliminates  the  need  to  go 
through  roll  after  roll  of  pen  plots  produced  by  the  Calcomp 
plotters. 

There  is  an  interactive  graphics  system  at  AFGL.  Located 
in  the  systems  area  are  two  Tektronix  4014  graphics  terminals. 
The  Tektronix  terminals  are  available  for  fast  plotting  of 
small  (less  than  65K  words)  programs.  The  software  in  these 
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AFGL  CYBER  74  Computer  System  (cont'd) 
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programs  is  similar  to  that  used  above  except  with  the  added 
ability  for  the  user  to  interact  with  his  program.  For  ex¬ 
ample,  if  a  scientist  sits  at  the  terminal  and  produces  a 
plot;  he  can  replot  the  data  "blowing  up"  a  section  that  he 
find  most  interesting.  Hard  copies  are  produced  by  the 
pressing  of  a  button. 

The  primary  language  used  on  the  6600  is  Fortran  extended. 
It  is  by  far  more  efficient  than  the  conventional  Fortran  IV. 
Numerous  non-ANSI  (American  National  Standards  Institute)  ex¬ 
tensions  have  been  introduced,  including;  implied  typing, 
multiple  system  texts,  and  various  masking  and  shifting  func¬ 
tions.  For  a  complete  list  of  Fortran  extended  capabilities 
the  reader  is  referred  to  the  "6600  Series  Fortran  Extended 
v5  Reference  Manual" . 

The  assembler  language  used  with  the  6600  is  COMPASS  v4.0. 
In  some  instances  where  bit  manipulation  is  heavily  used,  or 
when  the  timing  of  an  inner  loop  is  of  paramount  importance, 
COMPASS  is  more  desirable  to  use  than  Fortran. 

DPSI  has  used  all  the  programming  capabilities  offered 
by  the  CDC  system.  The  user  of  our  programs  will  find  FORTRAN, 
INTERCOM,  COMPASS,  and  OFF-LINE  PLOTS,  etc.  utilized  through¬ 
out  the  program  library  and  thus  modification  of  these  pro¬ 
grams  require  all  personnel  to  be  well-versed  in  all  the  CDC 
systems. 
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2*1  LYC  CYBER  74  Programming 


LYC  provides  research  and  advanced  development  to  the 
USAF  in  the  area  of  cloude  physics.  Practical  applications, 
rather  than  pure  research  are  emphasized  as  the  mission  of 
the  branch.  At  present  these  practical  applications  are  in 
the  following  areas: 

-  aircraft  icing  rates 

-  melting  layer  research 

-  radar  correlation 

-  Passarelli  spiral 

-  Mathematical -Physical  cloud  modeling 

-  snow  rate  study 


The  following  chapters  describe  the  contributions  that  DPSI 
has  made  in  the  above  areas  during  the  length  of  this  contract 

Specific  programs  are  listed  under  each  general  applica¬ 
tion.  For  every  program  the  inherent  mathematical  modeling 
will  be  fully  described  or  referenced.  In  addition  a  sample 
output  description  will  be  included. 

Full  operating  instructions  are  presented  for  each  pro¬ 
gram.  The  user  of  this  information  should  be  familiar  with 
the  contents  of  the  "AFGL  USER'S  GUIDE".  That  document  des¬ 
cribes  most  of  the  conventions  that  must  be  compiles  with  in 
order  to  run  jobs  on  AFGL  CYBER  systems.  Although  the  user's 
guide  is  a  much  simplified  version  of  the  set  of  CDC  cyber 
manuals  it  gives  enough  information  for  the  running  of  simple 


2.1  LYC  CYBER  74  Programming  (cont'd) 


jobs.  This  knowledge  plus  the  detailed  operating  instruc- 
,  tions  presented  in  this  document  are  quite  sufficient  for 

the  proper  running  of  the  programs. 

DPSI  maintains  all  CYBER  programming  on  a  disk  pack  (LYCPOF) 
owned  by  the  computer  center  (SOD) .  In  the  succeeding  docu¬ 
mentation  DPSI  uses  the  convention  of  attaching,  in  the  in¬ 
structions,  the  word  BIN  to  the  program  "ATTACH"  to  make  the 
compiled  program  local  to  a  user  job. 

,  ex.  ATTACH , LGO , PLTEXTRACTB IN , ID-KAPLANF , MR= 1 . 

LGO. 

The  binary  of  program  PLTEXTRACT  is  made  local  to  the 
|  job  as  a  file  LGO.  Repetition  of  the  name  tells  the  system 

to  load  it  and  begin  execution.  DPSI  maintains  several 
binaries  on  the  shared  system  (ID=MILLERP , ROBERTSK , KAPLANF ) 
using  this  convention.  However  not  all  programs  are  so  saved. 

|  The  computer  center  does  not  allow  files  to  be  stored  on  the 

shared  disk  system  without  their  continued  usage.  For  this 
reason  DPSI  maintains  only  the  most  frequently  used  binaries 
j  on  the  disk.  If  the  user  wants  a  program  not  currently  on 

|  the  shared  system  (determined  by  use  of  the  "audit"  command) 

there  are  two  options.  INTERCOM  can  be  used  to  attach  the 
program.  This  will  compile  and  save  the  binary  on  the  shared 
system  as  above.  Or,  the  control  deck  can  be  altered  by 

,  a 

I  placing  a  "PK"  parameter  on  the  first  card  and  replacing  the 

original  attach  with  the  following  cards: 

ex.  PAUSE.  PLS  MOUNT  DISK  LYCP07. 

MOUNT , VSN=LYCP  0 7 , SN=L YC  P  0  7 . 

I  ATTACH, P, PLTEXTRACT, ID=LALLY,SN=LYCP07. 

FTN,I=P,PL=999999. 

•  LGO. 
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2.1  LYC  CYBER  74  Programming  (cont'd) 


In  this  example  the  compilation  is  actually  done  during 
the  job.  In  the  succeeding  documentation  the  single  attach 
convention  will  always  be  used.  It  is  left  to  the  user  to 
determine  if  the  particular  program  is  or  is  not  currently 
on  the  shared  system. 


Airborne  Data  Collection  Systems  Description 

During  the  last  data  collection  period,  ending  in  1981, 
two  aircraft  were  used  to  perform  the  atmospheric  sampling 
necessary  for  the  desired  analysis  of  the  independent  research 
projects.  A  USAF  owned  and  operated  MC-130E  was  used  for 
data  sampling  to  approximately  30,000  feet.  High  altitude 
observations  to  45,000  were  achieved  with  a  USAF  leased  Lear 
jet  model  36C. 

Each  aircraft  had  on-board  equipment  with  the  capability 
to  derive  a  liquid  water  content  (LWC)  value,  as  well  as  other 
meteorological  parameters,  in  real  time.  This  was  accomplished 
by  utilizing  the  respective  on-board  computing  facilities  of 
each  aircraft  in  conjunction  with  the  data  obtained  from  the 
PMS-lD  systems.  In  addition  each  aircraft  has  three  indepen¬ 
dent  means  of  measuring  and  recording  LWC. 

The  following  sections  of  this  chapter  delineate  a  des¬ 
cription  of  the  equipment  used  and  their  functions.  See  fig¬ 
ures  3.1  and  3.2  for  a  data  flow  diagram  of  the  on-board 
systems. 
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3.1  PMS-1D  sizing  system 

The  PMS-lD  particle  sizing  system  consists  of  three  dif¬ 
ferent  probes  that  record  particle  counts  in  overlapping  size 
ranges.  The  Axial  Scattering  Probe  detects  particles  in  the 
2-30y  range.  The  20-300y  particles  are  measured  with  the 
cloud  probe.  The  Precip  (or  Precipitation)  Probe  is  used  for 
particles  in  the  200-4500u  range.  Actual  size  ranges  for  the 
probes  in  each  aircraft  are  shown  in  the  following  chart. 


Axial  Scatter 
Cloud 

Precipitation 
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2-30y 

20-300y 

300-4500y 


LEAR 
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20-300y 

200-3000y 


The  axial  scattering  pro fee  is  considerably  different  than 
the  other  two.  It  measures  optical  forward  scattering  from 
small  particles  in  a  constant  size  sampling  volume.  Dual  photo¬ 
diode  detectors  are  used  to  verify  that  the  particles  are 
within  the  sampling  volume.  There  is  pulse  height  detection 
circuitry  to  classify  the  particles  into  fifteen  size  cate¬ 
gories.  The  size  classes  are  then  read  out  by  the  data  acqui¬ 
sition  system  at  one  second  intervals. 

This  probe  is  specifically  designed  for  water  particles 
only.  Since  the  scattering  function  of  ice  crystals  is  poorly 
understood,  the  probe  is  used  only  to  indicate  relative  num¬ 
bers  but  is  not  normally  relied  upon  for  determining  mass  of 
ice  crystals. 


The  cloud  and  precip  probes  utilize  a  laser  beam  condensed 
and  mirrored  to  a  zoom  lens  which  distributes  light  to  a  row 
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3.1  PMS-lD  sizing  System  (cont'd) 

of  diode  sensors;  the  cloud  probe  has  22  sensors  and  the  precip 
24.  As  the  aircraft  flies,  particles  appear  between  the  zoom 
lens  and  the  sensors,  interrupting  the  light.  A  shadow  is 
cast,  shutting  off  some  of  the  diode  sensors.  The  device  is 
read  when  a  diode  is  turned  off,  and  the  sampling  is  continued 
as  any  diode  state  changes  until  all  diodes  are  back  on.  At 
the  conclusion,  then,  a  particle  of  known  diode  length  has 
been  counted;  the  output  consists  of  the  count  of  particles 
seen  for  lengths  of  one  to  fifteen  diodes  for  each  second. 

As  an  example,  if  a  particle  traces  the  following  states  in 
the  diodes  (0  =  diode  on;  X  =  diode  off) 

0000000000000000000000 

oooooooxxooooooooooooo 

ooooooxxoooooooooooooo 

ooooxxxooooooooooooooo 

0000000000000000000000 

Figure  3.3:  Particle  trace  through  PMS-1D  system 

the  result  would  be  1  particle  of  diode  length  5.  The  "5" 
results  from  the  maximum  different  number  of  diodes  turned 
off  from  start  to  finish  of  the  sampling  (all  within  a  frac¬ 
tion  of  a  second) .  The  reader  should  note  that  this  ID  PMS 
device  will  suffer  errors  when  two  particles  are  seen  at  the 
same  time,  for  then  the  maximum  lengths  will  be  added  together. 
Of  course  the  probability  of  the  occurrence  is  low. 


3.2  PMS-2D  Sizing  System 

The  PMS-2D  hardware  configuration  is  described  below. 

C130 

Two  PMS-2D  particle  display  systems  with  two  size  ranges 
(25-800y  and  200-6400y)  tied  to  dual  Pertec  (model  F5640-9) 
digital  recorders.  The  recent  PMS-2D  data  collection  system 
aboard  a  Corvair  990  consisted  of  60-1920y  and  240-F680y 
particle  size  ranges. 

LEAR 

Two  PMS-2D  particle  display  systems  with  two  size  ranges 
(40-1280y  and  160-512y)  tied  to  a  Pertec  (model  T7640-9) 
digital  recorders. 

The  2D  Knollenberg  has  some  important  advantages  over 
the  earlier  ID  model.  Firstly,  there  are  32  sensors  each 
exactly  25y  in  diameter.  The  second  dimension  is  achieved 
by  taking  readings  over  time  so  that  a  two  dimensional  pic¬ 
ture  of  the  shadow  is  made.  The  sampling  rate  is  adjusted 
to  the  speed  of  the  aircraft  so  that  a  reading  would  be  taken 
every  25y  of  length.  That  is,  if  the  aircraft  flies  at  100 
meters/sec,  the  sampling  rate  would  have  to  be  4  megahertz. 
This  exact  ratio  cannot  be  maintained  perfectly,  so  the  re¬ 
sults  are  modified  slightly  in  the  computer  according  to  the 
true  airspeed  of  the  aircraft.  Like  the  ID  device,  the  out¬ 
put  is  turned  on  when  a  sensor  is  shut  off,  and  continues  un¬ 
til  all  sensors  are  back  on,  but  this  device  will  output  the 
status  of  each  of  the  32  sensors  every  four-millionth  of  a 
second  until  all  the  sensors  are  back  on.  Thus  the  2D  device 
gives  a  picture  of  the  particle(s)  as  subsequent  readouts  are 
placed  together,  and  will  i>ot  give  incorrect  results  when  two 


3.2  PMS-2D  sizing  system  (cont'd) 

particles  are  seen  simultaneously. 

An  additional  advantage  of  the  2D  system  is  the  end  re¬ 
jection  feature.  If  a  particle  occludes  either  ending  diode, 
it  is  still  recorded  on  magnetic  tape.  On  the  ID  system  these 
particles  are  not  counted  and  the  data  is  lost. 

The  reason  for  this  rejection  is  the  philosophy  of  the 
ID  system;  if  the  ending  diode  is  occluded,  there  is  no  way 
to  estimate  the  true  particle  length  and  it  was  felt,  at  the 
engineering  design  level,  that  it  would  be  better  to  eliminate 
the  particle  rather  than  counting  it  as  one  with  a  lesser 
diameter.  With  that  in  mind,  the  end  diode  rejection  feature 
was  incorporated  into  the  ID  system.  An  area  of  concern,  how¬ 
ever,  is  the  number  of  particles  being  rejected;  with  the  pre¬ 
sent  ID  system  there  is  no  way  to  determine  this.  It  may  be 
of  considerable  consequence  because  the  larger  particles  have 
the  greatest  probability  of  being  rejected,  and  it  is  precisely 
these  particles  which  will  contribute  heavily  to  the  liquid 
water  content  and  radar  reflectivity. 


3.3  Liquid  Water  Content  Measuring  Systems 

Each  aircraft  had  the  capability  to  derive  a  calculated 
liquid  water  content  value,  (as  well  as  numerous  other  meteor¬ 
ological  parameters)  LWC,  in  real  time.  This  was  accomplished 
by  utilizing  the  respective  on-board  computing  facilities  of 
each  aircraft  in  conjunction  with  the  data  obtained  from  the 
PMS-lD  systems.  In  addition  each  aircraft  had  two  independent 
means  of  measuring  and  recording  the  LWC. 

On  each  aircraft  there  was  a  Johnson-Williams  (model  LWH) 
liquid  water  content  meter  that  measured  the  concentration 
of  water  in  cloud  sized  droplets.  The  device  used  calibrated 
resistance  wires,  both  normal  and  parallel  to  the  airflow,  to 
measure  the  amount  of  cooling  due  to  the  evaporation  of  the 
impinging  water  dropletj$f  It  had  a  specified  range  of  0-6 
g/m3 ,  however,  it  was  considered  accurate  only  for  particles 
in  the  10-50p  class. 

On  the  C130  there  was  an  additional  device  to  measure 
LWC  known  by  the  acronym  EWER  (Evaporate  the  Water  that  ag¬ 
gravates  Erosion  on  Reentry) .  This  device  took  two  samples 
of  the  free-stream  airflow  simultaneously.  One  sample  was 
passed  through  a  high-powered  heater  to  melt  and  evaporate 
all  the  condensed  water.  In  the  other  sample  the  condensate 
was  mechanically  separated  from  the  air  by  a  centrifuging  ac¬ 
tion  just  inside  the  air  inlet.  A  small  portion  of  the  evapo¬ 
rated  sample,  and  all  of  the  separated  sample,  were  passed 
through  small  heaters  to  adjust  the  final  temperature  before 
entering  the  detector  chamber . 


3.3  Liquid  Water  Content  Measuring  Systems  (cont'd) 

The  samples  were  passed  between  the  source  and  receiver 
of  a  Lyman-Alpha  humidiometer .  The  ultraviolet  energy  emitted 
by  the  source,  was  attenuated  by  the  water  vapor  in  the  path 
between  the  source  and  receiver.  The  relative  attenuation  of 
the  UV  was  a  direct  measure  of  the  water  vapor  content  in  g/m3. 
The  water  vapor  in  the  evaporated  sample,  and  also  the  back¬ 
ground  water  vapor,  was  optically  sampled.  A  provision  was 
also  made  for  varying  the  gap  between  the  source  and  receiver 
so  that  the  best  sensitivity  in  various  ranges  of  LWC  can  be 
adjusted  in  flight. 

The  EWER  system  was  under  the  complete  control  of  a  ROLM 
(model  1601)  computer*.  In  addition  to  performing  all  the 
sequencing  operations,  the  computer  also  performed  the  real 
time  calculations  of  liquid  water  content.  The  shortcoming 
of  this  system  was  the  method  of  data  recording 

The  EWER  LWC  analog  output  was  directly  tied  to  one  of 
the  digital  VCO  channels  on  the  PMS-1P  system.  From  there 
it  was  recorded  on  both  the  PMS-1D  and  TU-10  flight  tapes. 

By  applying  the  correct  calibration  equation  to  this  data, 

LWC  values  were  processed  as  a  function  of  time  in  both  a 
real  time  or  post  processing  mode.  This  data  was  either  pre¬ 
sented  in  tabular  or  graphical  form. 

There  was  an  additional  device  to  measure  LWC,  known  by 


*DPSI  recognizes  this  machine  as  a  rugqedized  version  of  the 
Data  General  NOVA  computer.  The  remainder  of  this  text  refers 
to  it  as  such. 


3.3  Liquid  Water  Content  Measuring  Systems  (cont'd) 

the  acronym  TWCI  (Total  Water  Content  Indicator) .  The  device 
was  designed  to  measure  total  water  content  in  both  liquid 
and  solid  cloud  particles. 

A  chamber  was  exposed  to  the  air-flow  where  water  and 
ice  crystals  are  collected.  Then  butyl  carbitol  was  intro¬ 
duced  to  the  chamber.  The  sensor  used  a  capacitor  in  a  tuned 
circuit  with  the  fluid  as  the  dielectric.  Water  in  the  fluid 
changed  the  capacitance  and  the  output  frequency  of  the  tuned 
circuit.  Electronic  circuits  detected  the  shift  by  comparing 
it  with  an  identical  cell  utilising  pure  butyl  carbitol  as 
the  dielectric.  The  output  of  this  system  was  calibrated  in 
terms  of  water  content. 


3.4  Meteorological  Parameters  -  VCO  Data 

In  addition  to  the  computers,  their  peripherals,  particle 
information  systems,  and  various  liquid  water  content  devices, 
there  were  a  number  of  other  data  sensors  on-board  each  air¬ 
craft.  These  devices  were  fed  directly  to  voltage  controlled 
oscillators  (VCO's)  and  were  hard  wired  into  the  output  stream 
by  analog  to  digital  converters.  These  sensors  contained  en¬ 
vironmental  information  of  pressure,  temperature,  and  dewpoint 
and  also  aircraft  information  of  heading  and  airspeed.  The 
data  were  multiplexed  into  both  the  PMS-lD  and  PMS-2D  buffer 
memory  systems  and  recorded  redundantly  on  their  respective 
digital  tapes. 

The  placement  of  the  output  VCO  data  which  were  combined 
with  the  PMS  data  varies  from  aircraft  to  aircraft  and  from 
time  to  time.  The  numbers  themselves  were  limited  to  four 
decimal  digits  each,  giving  a  range  of  0  to  9999.  These  VCO's 
were  counts  and  had  to  be  converted  to  proportionate  values. 


4.  Meteorological  &  Cloud  Physics  Data 
4.1  Production  Job  Stream  Overview 

This  section  graphically  illustrates  the  interrelation¬ 
ship  of  existing  programs  in  the  product  job  stream.  Figures 
4.1  through  4.5  depict  the  logical  processing  flow  within  the 
general  headings  presented.  Cross  purpose  use  of  programs 
for  different  headings  is  quite  often  done.  Programs  are 
listed,  however,  under  the  system  that  most  often  use  them. 


Analysis  And  Correlation  Program  Flow 


RADMCOM 
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4.2  PMS-1D  Data  Processing 


For  the  most  part,  the  programs  concerned  with  the  pro¬ 
cessing  of  LYC’s  one-dimensional  data  have  not  been  changed 
during  the  period  of  this  contract  and  hence  will  not  be 
reported  on  herein.  Descriptions,  operating  instructions, 
and  sample  output  for  the  programs  may  be  found  in  the  Final 
Report  for  Contract  F19628-78-C-0131. 


The  following  sections  describe  those  programs,  dealing 
with  the  one-dimensional  data,  which  have  been  revised  during 
the  contract  period  (August  1981-May  1983) , 
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4.2.1  Program  EXAUST 


At  the  request  of  AFSC/SD  (Air  Force  Systems  Command/ 

Space  Division) ,  the  cloud  physics  branch  of  AFGL  used  its 
PMS-1D  probes  to  make  measurements  of  particle  distributions 
for  the  recent  space  shuttle  launch. 

During  the  March  22,  1982  launch  of  NASA's  Space 
Shuttle  the  LYC  scatter,  cloud  and  precipitation  probes  were 
strategically  placed  100  meters  from  the  base  of  the  launch 
pad  atop  a  12  foot  pole.  When  the  solid  rocket  boosters  ig¬ 
nited  and  the  rocket  lifted  off,  LYC  probes  were  directly  in 
the  path  of  its  plume.  Raw  data  collected  by  the  probes  is 
written  to  a  PMS-lD  data  tape.  This  data  tape  will  be  returned 
to  AFGL/LYC  for  post  processing  by  the  AFGL  CDC  CYBER  Computer. 

A  new  program,  program  EXAUST,  has  been  written  to  process 
this  data.  Using  the  PMS-lD  tape  as  input  EXAUST  reads  in 
the  raw  channel  counts  from  each  probe.  The  following  informa¬ 
tion  is  derived  from  each  probe;  channel  number  densities 
(N/CM**3-viM)  ,  channel  normalized  number  densities  (N/M**3-MM) 
and  their  totals  are  calculated  and  printed  as  output.  No 
channel  number  adjustments  or  melting  equations  are  used  so 
the  diameters  are  a  straight  forward  calculation  (center  di- 
amter  =  channel  #  *  Diode  Size) . 

An  optional  output  tape  has  been  included  to  provide 
processed  data  to  NASA  for  use  in  their  modeling  studies. 
Complete  operating  instructions,  output  tape  format,  and 
sample  output  follow. 


4. 2. 1.1  Program  EXAUST  Operating  Instructions  and  Sample  Output 

COMMAND  DECK 

DPSI,TJ00,TP1,NT1.  2*16  ID 

ATTACH,L60,EXAUSTIIN,IP*  KAPLANF 
NAP, OFF. 

VSN,TAPEI*PHSXXX. 

RE QUEST,! APE  1 ,S, Ml ,HT,N0R1N6. 

VSN,TAPE2*LTCXXX/NT.  (OPTIONAL  OUTPUT  TAPE) 

REQUEST, TAP£2,N, RING, NT.  (OPTIONAL  OUTPUT  TAPE) 

flLE<Uni  ,RT*U,BT«K,HRL*1024,NBL*1024,RR«1  ,»FS*105) 

LPSET, PRESET* ZERO. 

L6D. 

EX1T(U) 

7/8/9 


DATA  CARDS 


ID  CARr 
1-10 
11-20 
21-26 


FLT  ID  IN  FORM  CXX-XX 

FIT  DATE  IN  FORM  DD  HON  YR 

PMS  TAPE  NUMBER  IN  FORM  PNSXXX 


OPTION  CARD 

5  ICLK  -  0  USE  A/C  CLOCK  (USED  TO  TEST) 
1  USE  PAS  CLOCK 

10  IAVE  -  AVERINS  INTERVAL  IN  15  FORMAT 
11-18  PMS  ON  TIME  IN  FORM  MH  MM  SS 
21  LYC  -  0  NO  OUTPUT  TAPE 

1  CREATE  OUTPUT  TAPE 


PASS  CARDS  -  AS  MANY  AS  DESIRED  IN  INCREASING  ORDER 
1-6  PASS  START  TIME  IN  FORM  HH  MM  SS 
8-13  PASS  STOP  TIME  IN  FORM  «  MM  SS 


A/7/6/9 


I 


OUTPUT  TAPE  FORMAT 
RDI  CONTENT 


1 

2 

3 

4-18 

1  f -33 

34-48 

49-43 

44-78 

79-93 

94-108 

109-123 

124-138 


START  TIME 
STOP  TIME 
VINO  SPEED 
RAU  COUNTS  SCATTER 
RAN  COUNTS  CLOUD 
RAU  COUNTS  PRECIP. 
NUMBER  DENSITIES  SCAT 
NUMBER  DENSITIES  CLOUD 
NUMBER  DENSITIES  PRECIP 
NORMALIZED  ND  SCAT 
NORMALIZED  ND  CLOUD 
NORMALIZED  ND  PRECIP 


FORMAT  TYPES 


A10  START, STOP  TINES 

E10.5  UIND  SPEED 

110  RAU  CHANNEL  COUNTS 

E10.5  NORMALIZED  AND  UNNDRMALIZE*  ND 


ONE  RECORD  *  138  UORDS  *  1380  CHARACTERS 


SSSS  TTTTT  SSSS  3333  EEEEE  *  X  A  U  0  SSSS  TTTtt 


ZZZZZ 
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4.2.2  Program  SOUND 


Program  SOUND  is  designed  to  provide  graphical  and  numeri¬ 
cal  comparisions  of  aircraft  vs  Radiosonde  temperature  data. 

When  comparing  aircraft  temperature  values  with  Radiosonde 
data  it  is  necessary  to  reduce  the  array  of  aircraft  data  to 
discrete  values  for  any  given  pressure  level.  Thus  after  the 
temperature  data  is  read  from  the  input  tape  it  is  internally 
sorted  by  pressure.  Median  values  can  then  be  obtained  for 
every  10MB  pressure  band.  It  is  these  median  values  which 
are  used  for  comparison  and  analysis.  Median  values,  as  op¬ 
posed  to  a  mean,  were  used  because  they  are  less  affected  by 
extraneous  data  points. 

The  first  plot  produced  by  program  SOUND,  is  the  raw 
data  plot.  The  points  will  always  be  plotted  in  ascending 
order,  visual  inspection  can  verify  the  sorting  procedure. 

After  this  (or  any  plot)  a  hard  copy  can  be  obtained  by  hit¬ 
ting  the  'copy'  button  on  the  tektronix  terminal.  The  next 
plot  will  not  be  displayed  until  a  space  followed  by  a  carriage 
return  is  typed  on  the  keyboard.  This  allows  time  to  produce 
the  hard  copy  between  plots.  If  the  user  enters  the  letter 
"S"  in  place  of  a  space  the  program  will  jump  to  the  ' STOP(Y/N) ' 
prompt.  A  'Y*  will  stop  the  job  and  return  to  intercom.  An 
'N'  will  return  the  user  to  the  beginning  of  SOUND.  In  a 
full  sequence  for  any  given  data  set,  four  plots  and  one  page 
of  output  are  produced.  In  order;  the  four  plots  consist  of 
(1)  raw  data,  (2)  raw  data  with  a  median  value  line  overlaid, 

(3)  raw  data  with  median  and  Radiosonde  data  overlayed  (Radio¬ 
sonde  data  distinguishable  as  a  dashed  line) ,  and  (4)  a  plot 


4.2.2  Program  SOUND  (cont’d) 


of  median  values  and  Radiosonde  data  without  raw  values.  In 
addition,  dry  adiabatic  lines  are  superimposed  over  the  plot 
to  aid  in  the  analysis  of  the  change  in  temperature  of  an 
air  parcel  due  only  to  the  change  in  pressure.  Operating 
instructions  and  sample  output  follow. 
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4. 2. 2.1  Program  SOUND  Operating  Instructions  and  Sample  Output 

turn  on  computer  and  hard  copy  unit 
follow  'LOGIN'  procedure 
USER  INPUT 

ATTACH , LGO , SOUNDBIN , ID=OTOOLE , MF=1 . 

ATTACH ,  TAPEl ,  TEMPCKDATA ,  lJ>=OTOOLE,  MR=  1 

ATTACH , TEK , TEKSIM 

LIBRARY , TEK 

REQUEST , TAPE 6 , *Q 

ETL, 500 

LGO 

refer  to  user  dialogue  on  next  page  for  examples 
of  prompts  and  replys. 

DISPOSE , TAPE6  ,  PR , IAC 
LOGOUT 


ENTER  IN  FIT  1DFLT  F79-31 
ENTER  IN  FLT  DATE 16  DEC  79 
CHOOSE  TEMP  TYPE  — 

1  .  .  TRUE  TEMP 

2  .  .  CALC  TEMP 

3  .  .  DEN  POINT 
ENTER 

DEFAULT  LIMITS  FOR  TRUE  TEMP 
ARE  *20.  TO  30. 

00  YOU  NISH  TO  ALTER  THE  L1MITSYES 
ENTER  NEM  LIMITS  (HIN.HAX>-lB7f"~ 
DEFAULT  LIMITS  FOR  PRESSURE 
ARE  1050.0  TO  400.0 
00  YOU  NISH  TO  ALTER  THE  LIHITSYE9 
ENTER  NEU  LIMITS  < M IN , MAX >900760% 
SHOULD  INPUT  FILE  BE  REMOUHDTVTKTT.  . .  Y 


DC  YOU  NISH  TO  INPUT  RAOB  DATA<Y,N>. . .  Y 
ENTER  IN  •  OF  LEVELS...  IX 

AFTER  EACH  LEVEL  i  ENTER  PRESS, TEMP, DEWPOINT. 
1 >  -->8?4.-2.3,-2.3 
27  ~>87?«-1.3,rl.3 


27 
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s?  ~>mi 

67  ~>Z2J 
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97  — >Z£ 
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CO 
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637.00 

<MB> 

-11.10 

<c> 

-26.70 

CC> 

LEVEL 

• 

13 

601.00 

(MB) 

-13.70 

<c> 

-27.70 

CO 

DO  YOU 

WISH 

TO  ALTER  ANY 

VALUES (Y, 

N>..«1L 

47 
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iunss3Ud 


sample  plot 
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DEW  POINT 


4.2.3  Program  FILTER 


Several  changes  have  been  made  to  program  FILTER.  Orig¬ 
inally  this  program  was  written  to  filter  out  persistance  (non 
independence)  of  LWC  icing  data. 

The  program  was  modified  to  include  the  option  of  plotting 
either  LWC  or  true  temperature  vs.  distance  (km)  and  time  (sec). 
This  modification  also  involved  changing  the  plot  axis  and 
output  to  be  labeled  with  appropriate  captions. 

The  program  was  modified  to  calculate  and  print  out  a 
smoothed  filtered  power  along  side  the  filtered  power.  Each 
value  is  printed  in  tabular  form  with  a  one  to  one  correspon¬ 
dence.  The  equation  for  calculating  smoothed  filtered  power 
was  provided  by  the  scientist  and  is  as  follows: 

DATA(I-l) * . 23+  DATA ( I ) * . 54  +  DATA(I+1) * . 23 

For  the  sake  of  consistency,  the  scientist  also  requested 
that  the  plots  of  the  output  be  put  on  a  static  scale  which 
eliminates  the  need  to  calculate  appropriate  scales  each  time. 
Consult  the  sample  plots  to  examine  the  limits  used. 

Operating  instructions  and  a  list  of  sample  output  follow. 


4. 2. 3.1  FILTER  Operating  Instructions 


COMMAND  CARDS 

DPSI,CM65000,T500,NTl.  '  ACCT  «  NAME 

ATTACH ,  LGO , F ILTERBIN ,  ID=1  MILLERP ,  MR=1 
VSN,TAP£l=LYCXXX/NT.  (KNOLLlD  OUTPUT  TAPE) 
REQUEST , TAPE 1,E, NORING , NT . 

ATTACH , CRT , CRTPLOTS . 

LIBRARY, CRT. 

REQUEST , TAPE 39 , *Q . 

LDSET , PRESET* ZERO. 

LGO. 

EXIT  (U) 

DISPOSE ,TAPE39 , PM. 

7/8/9 

DATA  CARDS 

CARD  #1  (REQUIRED) 

IPLOT,ZMIN,ZMAX  (FREE  FORMAT) 

IPLQT  -  IF  1  TRUE  TEMP,  0  IF  LWC 
ZMIN  -  MINIMUM  PLOT  LIMIT  FOR  CENTERED  VALUES 
ZMAX  -  MAXIMUM  PLOT  LIMIT  FOR  CENTERED  VALUES 
CARD  2-  (N+l)  (N  PASSES  IN  ASCENDING  TIME  ORDER) 
HH:MM:SS  COL  1-8  PASS  START  TIME 
HHsMMsSS  COL  10-17  PASS  STOP  TIME 


6/7/B/9 


4. 2. 3. 2  List  of  FILTER  Sample  Output  Plots 

1.  Centered  LWC  Values  vs  Distance  (KM)  and  Time  (Sec) 

2.  Autocorrelation  COEFS  vs  Distance  (KM)  and  Time  LAG 

3.  Unfiltered  Power  Density  vs  Frequency  (cycles/KM) 

4.  Filtered  Power  Density  vs  Frequency  (cycles/KM) 

5.  Centered  Temp  Values  vs  Distance  (KM)  and  Time  (Sec) 

Each  plot  frame  also  includes  as  a  heading 

1.  start  time 

2 .  stop  time 

3.  averaging  interval 

4.  number  of  data  points 

5 .  LWC  mean 

6.  average  altitude 

7.  average  true  air  speed 


Sample  charts  follow. 
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4.2.4  JW-LWC  Corrections 


During  this  contract  period  the  algorithm  for  making 
JW-LWC  corrections  has  been  modified  so  that  specific  cor¬ 
rections  can  be  made  over  a  specific  time  interval.  Prior 
to  this  only  one  set  of  corrections  was  allowed  throughout 
a  run.  Now  more  accurate  corrections  may  be  made  by  breaking 
them  up  into  time  segments. 


The  programs  in  the  production  job  stream  that  use  JW-LWC 
corrections  are;  KNOLL ID,  ICEDMP,  and  ICEEX. 


The  new  format  of  the  JWADJ  namelist  are  provided  below: 


NAMELIST  JWADJ 

CONTAINS  HEIGHT  PROFILES  FOR  A  JW-LWC  ADJUSTMENT. 

ELEMENTS  ARE: 

L  NUMBER  OF  LEVELS  (DEFAULT  0-  NO  CORRECTION)  (MAXIMUM  10) 
HT  HEIGHT  OF  LEVEL  IN  KM'S  (HT(l)GT  HT(2)GT...GT  HT (L) ) 

XA  ORIGIN  OF  LEVEL  (ONE  PER  LEVEL) 

SLA  SLOPE  FROM  LEVEL  (I)  TO  LEVEL  (1+1)  (L-l  SLOPES  REQUIRED) 
LTIM  NUMBER  OF  TIME  CORRECTION  INTERVALS  (DEFAULT  0- 
NO  CORRECTIONS)  (MAXIMUM  10) 

JWTIM  INTEGER  ARRAY  (2.10)  CONTAINING  START  AND  STOP  TIMES 
IN  SECONDS 

SLATIM  SLOPE  OF  TIME  CHANGE 
XATIM  INTERCEP  OF  TIME  CHANGE 


4.2.5  Derivation  and  Implementation  of  Expotential  Functions 


As  part  of  the  requirements  under  this  contract,  DPSI 
prepared  an  Interim  Scientific  Report  (ISR)  for  submittal  to 
the  contract  manager  ,  AFGL/LYC.  This  report  was  a  compre¬ 
hensive  discussion  of  the  derivation  and  implementation  of 
the  Weibull  and  Gamma  distribution  functions  to  model  axial 
scattering  spectrometer  data  distributions. 

These  data  were  obtained  from  26  discrete  flight  passes 
made  by  the  Air  Force  Geophysics  Laboratory's  instrumental 
C130-E  aircraft  during  the  1979-1980  winter.  The  data  were 
collected  in  supercooled  clouds  associated  with  aircraft 
icing.  The  instrument  probe  was  designed  to  collect  water 
particle  data  in  the  2-30  micron  range  at  2  micron  intervals. 
The  resultant  15  data  values  (channels)  for  each  flight  pass 
comprised  the  raw  data  sets.  However,  due  to  the  excessively 
noisy  data  signal  present  in  channel  1,  these  data  were  eli¬ 
minated  from  the  analysis.  Thus  data  channels  2-15  were  used 
as  input  to  the  Weibull  and  Gamma  functions. 

It  was  noted  that  the  distribution  of  the  data's  number 
density  spectra  seemed  to  be  describable  by  either  the  Weibull 
function  or  the  Khrgian-Magin  form  of  the  Gamma  distribution. 

The  main  body  of  the  ISR  is  concerned  with  the  deriva¬ 
tion  of  these  two  distribution  functions.  First,  the  Weibull, 
defined  as: 
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4.2.5  Derivation  and  Implementation  of  Exponential  Function  (cont'd) 


c  /  \  .T  w  b-1  -axb 
f (x)  =  Nabx  e 

is  derived  by  solving  for  'a'  and  1  b'  via  the  maximum  likelihood 
method  and  the  Newton-Raphson  method,  respectively.  Similarly, 
the  Gamma  function,  defined  as: 

c  /  \  vra^  2  -ax 
f  (x)  =  N— x  e 

is  derived  by  using  the  maximum  likelihood  method  to  solve  for 
•a' . 


As  part  of  the  derivation  for  both  the  Weibull  and  Gamma 
functions,  the  independent  factor  N  was  introduced  to  enable 
'a'  and  'b'  to  be  independent  of  the  magnitude  of  the  distri¬ 
bution's  range,  and  dependent  only  on  the  nature  of  the  re¬ 
spective  function.  However,  solving  for  N  by  the  way  of  the 
maximum  likelihood  method  would  result  in  the  undefined  ex¬ 
pression 


and  thus  another  method  had  to  be  employed.  Therefore,  utiliz¬ 
ing  the  method  of  least  square  analysis  and  partial  derivatives, 
the  value  of  N  was  determined. 

The  analysis  section  of  the  ISR  identifies,  with  the  aid 
of  four  figures,  the  shortcomings  of  the  Weibull  and  the  over¬ 
all  superiority  of  the  Gamma  function  to  model  cloud  droplet 
spectra  data  distributions  . 


56 


I 


4.2.5  Derivation  and  Implementation  of  Exponential  Function  (cont'd) 

Program  WBLGMA  utilizes  these  functions,  as  derived,  to 
model  and  analyze  the  raw  channel  counts  (usually  2-15) 
collected  by  the  axial  probe.  Section  4 . 2. 5* 1  contains  sample 
output  from  the  program  for  one  time  pass.  The  first  line  of 
output  simply  re-states  flight  and  probe  data  for  that  pass, 
as  well  as  which  data  channels  are  being  used  as  input. 

The  main  body  of  the  output  in  the  following  table  is  the 
channe)  by  channel  comparison  of  the  two  modeled  data  distributions 
(G  &  W)  to  the  raw  data.  The  overall  performance  of  each  distri¬ 
bution  function  is  provided  by  the  LOG  RMS  valve  of  each. 

See  scientific  report  no.  1 (AFGL-TR-82-0309) ,  published 
30  Septmeber  1982,  for  the  complete  discussion  detailing  the 
deriviations  and  analysis  techniques. 


4. 2. 5.1  WBLGMA  Operating  Instructions  &  Sample  Output 


ATTACH, W,WBLGA 


ATTACH, W, WBLGMA, ID=MILLERP 

CONNECT, INPUT 

CONNECT, OUTPUT 

SCREEN  132 

LGO 


89.98  2.4  .004  2  15 
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4.3  Ice  Detector  Analysis 


Data  collection  by  ice  detector  can  be  characterized  by  three 
distinct  modes  of  operation. 

1)  sensing  mode 

2)  detecting  mode,  and 

3)  standby  mode. 

These  modes  occur  in  the  sequence  listed.  The  standby  mode 
however,  can  be  triggered  manually  at  any  time. 

A.  In  the  sensing  mode  the  probe  is  signaling  that  the  ice 

is  accumulating  on  it.  The  counts  will  range  in  the  mode 
from  6250  to  9318.  (The  latter  value  being  equal  to  4.8 
volts  at  which  point  the  accumulation  of  ice  triggers 
the  detecting  mode) .  If  icing  rates  are  very  low  it  is 
possible  that  the  prohf!  will  remain  in  the  sensing  mode, 
i.e.,  never  accumulating  enough  ice  to  go  into  the  de¬ 
tecting  mode  or  the  ice  will  evaporate  and  counts  will 
decrease  in  value  to  approximately  6250. 

In  the  detecting  mode,  an  ice  signal  is  generated  and 
heating  cycle  is  initiated.  During  this  period  random 
counts  are  generated.  Typically  values  under  5000  are 
generated  (see  C)  for  a  number  of  seconds.  All  these 
data  (i.e.,  exceeding  9818  and  less  than  5000  counts) 
are  to  be  disregared  as  data,  but  to  be  used  as  the  in¬ 
dicator  of  a  change  from  the  detecting  to  standby  mode. 

The  standby  mode  follows  the  detecting  mode  and  has  a 
range  of  counts  from  5350  to  5150.  The  most  frequently 
observed  values  is  5250  counts.  In  this  mode  all  ice 
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has  been  removed  from  the  detector,  and  the  DC  current 
to  the  probe  has  been  cutoff.  Typically  this  mode  lasts 
under  2  seconds. 

The  time  needed  for  the  ice  detector  to  return  from 
standby  mode  to  the  sensing  mode  depends  on  icing  con¬ 
ditions,  ambient  temperature  and  airspeed.  With  moder¬ 
ate  icing  7-10  seconds  are  typical.  The  return  to  the 
sensing  mode  occurs  when  counts  equal  6250. 

It  is  possible  while  in  the  sensing  mode  for  apparently 
random  counts  to  be  generated;  that  is,  9818  counts  has 
not  been  reached  nor  exceeded  and  apparently  random 
signals,  (similiar  to  those  occuing  when  the  instrument 
is  in  the  detecting  mode)  are  observed.  Following  the 
occurrence  of  these  random  counts,  the  counts  recorded 
return  to  values  indicating  the  unit  is  still  in  the 
sensing  mode. 

Random  counts  are  sudden  large  excursion  in  magnitude. 
Since  these  occur  primarily  in  connection  with  the 
change  from  sensing  to  detecting  mode  and  frequently 
occur  with  voltages  exceeding  5V,  causing  a  rollover 
in  counts,  (i.e.,  numbers  >  9999  begin  at  0001)  it  is 
useful  to  add  10000  counts  to  all  values  <  5000. 


4.3.1  Program  ICEEX 


A  flight  was  made  on  06  DEC  79  in  order  to  collect  data 
using  the  Ice  Detector  device.  As  a  "quick  look"  method  to 
examine  the  data,  a  plotting  program  PLTICE  was  written  (see 
section  on  PLTICE).  This  plot  coupled  with  KNIUTIL's  VCO 
listing  for  this  flight  enabled  certain  empirical  relationships 
between  the  recorded  data  and  the  behavior  of  the  Ice  Detector 
device  to  be  derived. 

Mr.  Morton  Glass  formulated  these  relationships  on  19 
DEC  79  and  an  updated  version  of  his  memo  is  given  on  the 
following  pages.  Program  ICEEX  incorporates  these  concepts 
and,  as  can  be  seen  in  figure  4.6,  is  at  the  heart  of  the 
icing  data  processing  and  analysis  stream. 

Data  from  the  Ice  Detector  is  given  in  VCO  position  ten 
of  the  PMSlD  Kennedy  tape.  Program  ICEEX  reads  this  tape, 
identifies  the  modes,  eliminates  erroneous/redundant  points, 
and  produces  an  edited  output  tape. 

ICEEX  inputs  a  flight  identification,  a  namelist  card 
to  change  the  default  VCO  calibrations,  and  a  set  of  start 
stop  time  ranges. 

For  each  start  and  stop  time  range  the  first  standby 
mode  is  found  by  performing  a  three  point  running  mean  until 
this  average  is  within  200  counts  of  5250.  The  midpoint  of 
this  mean  is  then  put  onto  the  output  tape  with  a  code  word 


4.3.1  Program  ICEEX  (cont'd) 

equal  to  one.  The  next  step  is  to  find  the  "trigger  value". 
This  value  is  the  first  count  whose  value  exceeds  6000,  thus 
marking  the  transition  from  standby  to  sensing  mode  (by  ig¬ 
noring  the  detecting  mode  completely) .  When  this  time  is 
found  (no  averaging  is  used  -  only  the  raw  count)  the  data 
record  is  sent  with  a  code  word  equal  to  two.  Now  the  ice 
detector  is  considered  to  be  at  the  start  of  the  sensing 
mode.  Until  the  counts ,go  below  the  "trigger  value",  exceed 
the  maximum  9818  counts,  or  degrade  into  four  consecutive 
"bad  points",  every  data  point  is  written  to  the  output  tape 
with  a  code  word  equal  to  three. 

The  cycle  of  finding  the  standby  mode,  trigger  value, 
and  then  all  sensing  values  is  repeated  until  a  time  value 
exceeds  the  input  stop  time. 

This  inputting  of  start  and  stop  times  is  terminated  by 
the  end  of  input  cards  or  an  end  of  file  on  the  PMS1D  tape. 

Every  point  in  the  sensing  mode  is  verified  by  extra¬ 
polating  from  the  previously  read  two  points.  If  the  input 
value  differs  by  more  than  400  counts  from  the  extrapolated 
value,  the  extrapolated  value  is  used.  The  extrapolating 
method  is  designed  to  compromise  between  changing  direction 
and  having  the  data  rise  rapidly.  The  extrapolated  point 
C(i)  is  given  by 

C ( suggested) =( (C (i-2) +C (i-1) ) /2+2*C (i-1) -C (i-2) )/2 
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Program  ICEEX  (cont'd) 

Since  ICEEX  lists,  on  the  printer,  each  value  as  it  is 
being  written  to  tape,  an  asterisk  is  printed  to  denote  the 
fact  that  it  is  a  "calculated"  and  not  a  "true"  value.  Ever 
second  of  extracted  data  as  well  as  PMSlD  raw  counts  are 
written  in  groups  of  59  words  to  the  output  tape.  Operating 
instructions, a  description  of  the  output  tape,  and  sample 
output  follow. 


4. 3. 1.1  Program  ICEEX  Operating  Instructions  and  Sample  Output 


CONTROL  CARPS 

J3BID*  Tl 000*  TP1 *NT1 .  JOB  «  PASSWORD 

ATTACH. LGO* ICEEXBIN. I DMMILLERP , MR=1 • 

VSN*TApE1*®MSTAPEN'JMBER.  (KENNEDY  TAPE  -  PUS  ID  TAPE) 

REQUEST*  TAPE1*S*HI*NT. 

VSN.  TA°E2*0'JTPUTTAPEN'JMBER^NT.  (PROCESSED  DATA  TAPE  -  OUTPUT) 

REQUEST*  TAPE2* PE* N*  RINS» NT . 

FILE  (TAPE1  *RT*'J»  BT*K*  MRL*1  120»MBL*1  1£0»  RB=i  *  BFS*1 15) 

LDSET » PRESET*2ER0. 

LSO. 

7'3'9 


DATA  CARDS 

CARD  I 

COL  1-10  FLIGHT  ID 

11-20  FLI5HT  DATE 

25  CLOCK  (1«A'C*2«PNS> 

30  *>'JT  A  1  IN  THIS  COLUNN  TO  INDICATE  USE  OF  KENNEDY  TAPE 

75-30  EXTRAPOLATION  BOUND  (15  FORMAT  -  DEFAULT  400) 

#  . 

CARD  2 

VCOCHAN  NAMELIST  -  VCO  CALIBRATIONS 
(SEE  VCOEF  IN  KNOLL ID) 

CARD  3 

JWADJ  NAMELIST  -  JW-LWC  ADJUSTMENT  PROFILES 
(SEE  JW AD J  IN  KNOLL ID) 

CARD  4  -  (4*N-1) 

N  °ASS  CARDS  IN  FORM  (STARTIN’?  IN  COLUMN  1>  HH-MM-SS-HH-MM-SS 

START  STOP  TIMES 


WORD 


CONTENTS 


I 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 
13 


14 

15-29 


30-44 

45-59 


TIME  OF  RECORD  (IN  SECONDS) 

FLIGHT  IDENTIFICATION 
FLIGHT  DATE 
LATITUDE  (NOT  USED) 

LONGITUDE  (NOT  USED) 

TEMPERATURE  (TRUE  IN  DEGREES  CENTIGRADE) 
TEMPERATURE  (TOTAL  IN  DEGREES  CENTIGRADE) 
DEWPOINT  (IN  DEGREES  CENTIGRADE) 

JW-LWC  (UNADJUSTED  GM/M**3) 

PRESSURE  (MILLIBARS) 

TRUE  AIRSPEED  (M/SEC) 

ICING  COUNT 

CODE  1=STANDBY  RECORD,  2=TRIGGER  RECORD, 
3-SENSING  RECORD 
ALTITUDE  (METERS) 

CHANNEL  COUNTS  FOR  CHANNELS  1-15  OF  THE  SCATTER 
PROBE 

CHANNEL  COUNTS  FOR  CHANNELS  1-15  FOR  THE  CLOUD 
PROBE 

CHANNEL  COUNTS  FOR  CHANNELS  1-15  OF  THE  PRECIP 
PROBE 


ALL  VALUES  ARE  IN  FLOATING  POINT  FORMAT.  REDCORD  IS  59  WORDS 
LONG  WRITTEN  BY  A  FORTRAN  BINARY  WRITE  STATEMENT  (WRITE (UNIT) LIST) 
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SENSE  25  04  13  6043  .1 

EOF  ON  CARO  INPUT  -  OR  INPUT  TAPE 


4.3.2  Program  ICEDMP 


The  purpose  of  this  program  is  to  generate  icing  data 
output  for  inclusion  within  technical  reports.  ICEDMP  pro¬ 
duces  output  efficiently  and  quickly.  It  is  highly  modular 
in  nature  which  allows  changes  to  be  made  quickly  with  a 
minimum  of  additional  programming.  In  addition,  ICEDMP  uti_ 
lyzed  the  XEROX  hard  copy  device  at  the  AFGL  Computer  Center 
in  order  to  facilitate  the  inclusion  of  quality,  reduced  size 
(8%  x  11)  copy  within  technical  reports. 

ICEDMP  uses  the  raw  PMS-1D  data  tape  for  input  and  pro¬ 
duces  three  distinct  output  files?  one  contains  axial  probe 
data;  the  second  contains  processed  data.  Additionally,  the 
standard  output  file  lists  information  singular  to  the  parti¬ 
cular  flight  being  processed. 

The  following  pages  provide  operating  instructions  and 
sample  output. 


69 


4. 3. 2.1  ICEDMP  operating  Instructions  ana  sample  uucpu 

COMMAND  DECK 

DPS  I , T300,TP1 , CM50000.  ID  «  ID  NAME 

ATTACH,L60, 1CEDMPBIN, ID*  KAPLANF , N R=  1 . 

ATTACH, TAPES, VCOCALS, ID*  KAPLANF, MR=1. 

VSN,TAPE1 *PHSXXX. 

R£0UEST,TAPE1 ,S,HI ,MT, NORING. 

FILE ( TAPE1 ,RT*U,IT*K,MRL*1 024, HBL«1024,Ri'1  ,IFS*105) 

LDSEJ, PRESET* ZERO. 

EXITED) 

REWIND  TAPE4,TAPE5,TAPE0,TAPE7. 

C0PY,TAPE4. 

COPY, TAPES. 

COPY.TAPEA. 

COPY,TAPE7. 

7/8/9 


DATA  CARDS 
ID  CARD 

CC  1-10  FIT  ID  IN  FORM  FIT  EXX-XX 
11-20  FLT  DATE  IN  FORM  DD  NON  YR 

21-20  PUS  TAPE  NUNDER  IN  FORM  PMSXXX 

52-59  PNS  ON  TINE  IN  FORK  HH  NH  SS 

NAMELIST  VCOEF  -  USEING  ARRAY  C<3,13) 

USED  TO  OVERIDE  VCO  CALIBRATION  FILE  -  SEE  KN0LL1D 

NAMELIST  JUADJ  -  USEIN6  L  AND  ARRAYS  HTEIO),  XA(10)  AND  SLA(IO) 
USED  TO  HAKE  JU-LUC  CORRECTIONS  -  SEE  KNOLL  ID 

OPTION  CARD 

CC  5  ICLK  -  0  USE  A/C  CLOCK 
1  USE  PMS  CLOCK 

10  IVEL  -  0  USE  TRUE  AIRSPEED 

1  USE  CALCULATED  AIRSPEED 
11-15  IAVE  -  AVERAGING  INTERVAL  IN  15  FORMAT 

20  PASSAVE  -  IF  SET  TO  1  IGNORE  IAVE  AND  GIVE 

PASS  AVERAGES 

25  1TMP  0  USE  VCO  TEMPERATURE  CALIBRATION 

1  USE  VCO  TEMPERATURE  CALIBRATION  AND 
TRUE  AIRSPEED 

20-30  NRTS  -  NUMBER  OF  RECORDS  TO  SKIP 
.0  -  SKIP  NO  RECORDS  (DEFULT) 

N  -  SKIP  N  RECORDS  BEFORE  PR0CCESSIN6  ANY  DATA 

79-80  CLOUD  PARTICLE  TYPE  (DEFAULT  *  RAIN) 

PASS  CARDS  -  AS  MANY  AS  DESIRED  IN  TIME  INCREASING  ORDER 
CC  1-0  PASS  START  TIME  IN  FORM  HHMMSS 

8-13  PASS  STOP  TIME  IN  FORM  HHMMSS 

20  WHEN  SET  TO  1  USE  THE  JU-LUC  FOR  ICE  THICKNESS 
CALCULATION ( DEFAULT  AXIAL  FOR  EACH  PASS) 


0/7/8/9 
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4. 3. 2.1  Program  ICEDMP  (cont'd) 

AXIAL  PROBE  OUTPUT  (TAPE  5)  -  See  Appendix  23B 
PLIGHT  ID 

FLIGHT  DATE  TOP  OF  EACH  PAGE 

AVERAGING  INTERVAL 


Z  TIME  -  START  TIME  OF  AVERAGING  INTERVAL 

CHANNEL  -  RAW  AVERAGED  AXIAL  PROBE  COUNTS 
COUNTS 

CONC  -  AXIAL  PROBE  NUMBER  DENSITY  IN  CUBIC  CENTIMETERS 


CONC 


15  COUNTS A 
A  SCATVOL 


r  COUNTS-, 
LCM**3  J 


COUNTS ±  *  NUMBER  OF  PARTICLES  FOR  CHANNEL  i 

AXIAL 

SCATVOL  =  [CROSS  SECTIONAL  AREA  (M**2) xAIRSPEED (M/SEC) ]*l.E-6 
DBAR  -  MEAN  DIAMETER  OF  AXIAL  PARTICLES 


15 

1  (COUNTS.) (DI AM.) 

i=l  1  1 

DBAR  =  -  -  y5 - 

I  COUNTS.  (TOTAL  COUNTS) 

i=l  1 


DIAM±  =  CHANNEL  ADJUSTED  CENTER  DIAMETER 


DBARV  -  AXIAL  PROBE  VOLUME  AVERAGED  PARTICLE  DIAMETER 


1/3 

=  ( (ALWC*SCATVOL)/ (0.52 359 8E-3*TOTCTS) )  '  *1.E3 


DBARV 


4. 3.2.1  Program  ICEDMP  (cont'd) 

ALWC  -  AXIAL  LWC  (G/M**c) 

SCATVOL  -  AXIAL  VOLUME  (M**3) 

TOTCTS  -  TOTAL  COUNTS 

0.5235988  =  u/6  x  p  (WATER  DENSITY  l.E-3  G/MM**3) 

DMAZ  &  DO  -  95%  AND  50%  AXIAL  MEDIAN  VOLUME  DIAMETER 

The  diameter  at  which  95%  (50%  of  the  water  is  contained 
in  channels  less  than  this  size;  5%  (50%)  of  the  water  is 
contained  in  channels  greater  than  this  size.  Consult  KNOLL ID 
documentation . 

ALWC  -  AXIAL  LIQUID  WATER  CONTENT  (G/M**3) 

„  15 

ALWC  »  £  p  l  CONC .  DIAM .  3 

6  i=l  1  1 

S-VOL  -  AXIAL  SAMPLING  VOLUME  (CM** 3) 

S-VOL  =  CROSS  SECTIONAL  AREA  (M**2) xAIRSPEED (M/SEC) xl.E6 
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4.3. 2.1  Program  ICEDMP  (cont'd) 


CLOUD  PROBE  OUTPUT  (TAPE7)  -  See  appendix  23C 
FLIGHT  ID 

FLIGHT  DATE  TOP  OF  EACH  PAGE 

AVERAGING  INTERVAL 


Z  TIME 

CHANNEL 

COUNTS 

ID— C 

1D-C 


-  START  TIME  OF  AVERAGE  INTERVAL 

-  RAW  AVERAGED  CLOUD  PROBE  COUNTS 

-  CLOUD  PROBE  NUMBER  DENSITY  IN  LITERS 
15  COUNTS i 

=  i=15  V0Li  *  1,E'3 


COUNTS ^  ~  CLOUD  CHANNEL  COUNTS 

VOL^  =  CLOUD  CHANNEL  VOLUME  (M**3) 

_  CHANNEL 

CROSS  SECTIONAL  AREA (M**2) xAIRSPEED(M/SEC) 


DBAR 

DMAX  SEE  ABOVE  (USE  CLOUD  INFORMATION) 

DO 

CLWC  -  CLOUD  LIQUID  WATER  CONTENT  (g/M**3)  (SEE  ABOVE) 

TLWC  -  TOTAL  LIQUID  WATER  CONTENT  (g/M**3) 

TLWC  *  ALWC  +  CLWC  -  (LWC  CONTRIBUTION  OF  ANY  OVERLAPPING 

CLOUD  CHANNELS) 

PCT  CLD  -  PERCENT  CONTRIBUTION  OF  THE  CLOUD  LWC  TO  THE  TOTAL 
LWC 


PCT  CLD 


_ CLOUD  CONTRIBUTION  ln0 

AXIAL  &  CLOUD  CONTRIBUTION 
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4. 3. 2.1  Program  ICEDMP  (cont'd) 


DEWP 


EPT 

EPT 


CONC 

ALWC 

JW-LWC 


-  DEWPOINT  (C) 

IF  FROST  >  0  THEN  DEWP  =  FROST 
ELSE 

DEWP  =  9 . 84E-4*FROST2+l . 1305*FROST-0 . 012 

-  EQUIVALENT  POTENTIAL  TEMPERATURE  (KELVIN) 

=  PT[LWs/CpTc] 

WHERE  T  =  TEMP  °C 
TK  =  TEMP  °K 
D  =  DEWP  °C 
DK  =  DEWP  °K 
P  =  PRESS (Mb) 

o  q  c  n 

=  TK(1000/P)  POTENTIAL  TEMPERATURE  (K) 

=  4. 187(734. -.5J.TK)  LATENT  HEAT  OF  VAPORIZATION  (JOULES) 
*  .622 (e/P-e)  SATURATED  MIXING  RATIO 

=  VAPOR  PRESSURE  (Mb)  AT  TK 

=  .461  ( JOULES/ °K  MOLE) (SPECIFIC  GAS  CONSTANT  FOR 
WATER  VAPOR) 

=  1.005  (JOULES/g  °K)  SPECIFIC  HEAT  OF  DRY  AIR  AT 
CONSTANT  PRESSURE 

=  DK  +  (Dz (T-D/(DZ-1543. 307T) 

ISENTROPIC  CONDENSATION  TEMPERATURE  (°K) 

-  AXIAL  NUMBER  DENSITY  -  (SEE  ABOVE) 

-  AXIAL  LWC  (SEE  ABOVE) 

-  JOHNSON-WILLIAMS  LWC  (g/M**3) 


JW-LWC 


(CALIBRATED  VALUE) 200 
TAS  OR  CAS 


When  desired  JW-LWC  may  be  adjusted  based  on  an  ALTITUDE 
profile  -  consult  KNOLL1D  documentation  for  more  information. 
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4. 3. 2.1  Program  ICEDMP  (cont'd) 

lD-C  -  CLOUD  NUMBER  DENSITY  (SEE  ABOVE) 

♦ICE-TK  -  ICE  THICKNESS  DERIVED  FROM  THE  AXIAL  LWC  PER 
NAUTICAL  MILE  (USES  THE  DENSITY  OF  WATER) 

ICETK  =  ALWC*. 185325 

*  AN  ALGORITHM  HAS  BEEN  DEVELOPED  TO  DETERMINE  WHEN  WE  ARE 
NOT  IN  AN  ICING  STATE.  WHEN  THIS  OCCURS  THE  WORDS  "NO  ICE" 
REPLACE  THIS  CALCULATED  VALUE. 

TO  BE  IN  AN  ICING  STATE  ALL  OF  THE  FOLLOWING  MUST  BE  TRUE: 

1)  TRT  <  -.5°C 

2)  JW-LWC  >  .01  g/m**3) 

3)  NOT  PASS  AVERAGING  ONLY 

A)  IF  ANY  OF  THE  FIVE  SECONDS  AFTER  INITIALLY  REACHING  6000 
COUNTS  (FOR  EACH  CYCLE)  ARE  CONTAINED  IN  THE  AVERAGING 
INTERVAL 

OR 

B)  ANY  VALUES  WERE  BELOW  6000  OR  GREATER  THAN  9818  (IN 
HEATING  CYCLE) 

OR 

C)  6000  <=  ALL  COUNTS  DURING  AVERAGE  <=9818  TEST  AVERAGE 
SLOPE  DURING  PERIOD.  IF  2  POINT  AVERAGE  OR  LESS  IGNORE 
THIS  TEST,  ELSE  BREAK  INTO  5  POINT  MULTIPLES  (4  OR  3 
SECOND  SEGMENT  IF  NECESSARY)  AND  USE  A  DERIVATIVE  ROU¬ 
TINE.  IF  THE  AVERAGE  DERIVATIVE  >=25  THEN  ICEING. 

4)  PASS  AVERAGING  ONLY  -  THERE  MUST  BE  AT  LEASE  ONE  ICING  CYCLE 

ICING  RATE  METER  -  PASS  AVERAGING  ONLY 

CYS/NT  =  NUMBER  OF  ICING  CYCLES  PER  NAUTICAL  MILE 

=  _ CYCLE  COUNT*  18 5 3. 25 _ 

(#  OF  SECONDS  OF  PASS) * (AIRSPEED) 

CM/NT  =  CYS/NT*. 05  =  THICKNESS  OF  ICE  BUILDUP  PER  NAUTICAL 
MILE 
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2.1  Program  ICEDMP  (cont'd) 


PROCESSED  OUTPUT  (TAPE  6)  -  See  appendis  23D 
FLIGHT  ID 

FLIGHT  DATE  TOP  OF  EACH  PAGE 

AVERAGING  INTERVAL 


Z  TIME 
ALT 
ALT 
P 

TAS  or 
CAS 

TAS 

CAS 

P 

TRT 

MCHSQ 


DELP 

IAS 
TOT  T 

TRT 

TRT 

FROST 


-  START  TIME  OF  AVERAGING  INTERVAL 

-  STANDARD  PRESSURE  ALTITUDE  (METERS) 

-  44307.69  -  11872(P)0*190284 
«  PRESSURE  (Mb) 

STRAIGHT  CALIBRATED  VCO 

-  TRUE  AIRSPEED  OR  CALCULATED  AIRSPEED  (KNOTS) 

-  STRAIGHT  CALIBRATED  VCO 

=  /1516. 4 ( TRT+  2  7  3 . 16) MCHSQ  +  3/(TRT+273. 16) /P 
=  PRESSURE  (Mb) 

=  TRUE  TEMPERATURE  (C) 

=  MACH  SQUARED  NUMBER 

5  (  x | DELPj  0.2857143  1 
=  DELTA  PRESSURE  (Mb) 

=  ( (0.001865IAS) -0.060149) *IAS+3. 96965 
=  INDICATED  AIRSPEED  (CALIBRATED  VCO  IN  KNOTS) 

-  TOTAL  TEMPERATURE  (C) 

STRAIGHT  CALIBRATED  VCO 

-  TRUE  TEMPERATURE  (C) 

=  (TOT  T+273.16)/(l+0. 1992*MCHSQ) -273.16 

-  FROST  POINT/DEW  POINT  (C) 

STRAIGHT  CALIBRATED  VCO 


FFFF  l  TTTTT  EEEEE  888 
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EPT  CONC  ALWC  JLWC  1D-C  1CE-TK  ICING  RATE  METER 
K  NCM-3  GM-3  GM-3  Nl-1  CM/NT  CYS/NT  CM/NT  PCT  SENSE 


4.3.3  Program  ICEDERIV 

One  facet  of  ice  detector  analysis  is  a  study  of  how 
the  ice  detector  VCO  values  vary  with  the  amount  of  water 
observed  by  the  JW-LWC  and  the  PMS1D  AXIAl  devices. 


At  first,  the  analysis  concerns  itself  with  the  rate  at 
which  Lhe  VCO  readings  vary;  in  other  words  the  derivative 
of  the  ice  detector  VCO  values.  The  first  problem  with  nu¬ 
merically  calculating  a  derivative  is  in  choosing  the  neighbor¬ 
hood  over  which  it  is  to  be  computed.  The  second  problem  is 
that  a  derivative  calculation  is  extremely  sensitive  to  minor 
fluctuations  in  the  neighborhood  in  which  it  is  computed. 

The  solution  to  the  first  problem  is  to  calculate  the 
derivatives  over  various  intervals  and  deduce  the  best  one 
to  use.  The  second  problem  requires  a  smoothing  function  to 
be  applied  to  the  data  and  then  utilize  a  derivative  method 
to  realize  the  rates  of  change. 

In  APPLIED  ANALYSIS,  by  C.  Lanczos  (Prentice  Hall  1961) 
a  method  to  smooth  the  data  and  to  compute  the  derivatives 
is  given  on  pages  321-324. 

The  method  calculates  a  smoothed  polynominal  of  degree 
2K  (where  K  is  the  number  of  points  on  each  side  of  the  point 
at  which  the  derivative  is  to  be  computed)  and  then  develops 
its  derivative.  The  general  formula  for  the  derivative  of 
f  at  x  with  K  neighbors  on  each  side  for  equally  spaced  one 
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4.3.3  Program  ICEDERIV  (cont’d) 

second  data  is  then 

+K  K 

fMx)  =  Z  at  (x+a) /s  1  a7) 
a=-K  a=l 

The  above  formula  works  for  K  t  2. 

This  method  requires  the  point  which  the  derivative  is 
to  be  computed  to  be  in  a  neighborhood  of  K  points.  Thus  for 
the  first  two  and  last  two  points  the  general  formula  fails. 

The  way  to  derive  these  values  is  to  compute  a  least 
square  second  degree  fit  and  use  its  normal  equations: 

f*(l)  =  (-21f  (l)+13f  (2)+17f  (3)  -9f  (4))/20 
f 1 (2)  *  (-Ilf (l)+3f (2)+7f (3)+f (3))/20 

The  last  two  points  (xmax,  xmax-1)  are  comparably  calculated 

fJ(xmax)  =  (21f (xmax-13f (xmax-1) -17f (xmax-2) +9f (xmax-3) )/20 
f1 (xmax-1)  =  (Ilf (xmax) -3f (xmax-1) -7f (xmax-2) -f (xmax-3) )/20 

ICEDERIV  was  originally  written  to  accept  an  ICEEX 
output  tape,  and  for  each  data  point  calculate  the  deri¬ 
vative,  using  2,  3,  4,  5,  6,  and  7  neighboring  points. 

The  table  printed  by  ICEDERIV  is  broken  at  every  new  icing 
cycle  and  contains  the  ICE  count,  time  and  derivatives. 


4.3.3  Program  ICEDERIV  (cont’d) 

Since  a  comparison  is  to  be  made  with  JW-WJC  or  PMS  ID 
AXIAL  these  other  two  parameters  were  added  to  the  table. 
These  data  sources  are  also  subject  to  one  second  ''noise'' 
and  were  smoothed.  A  five  point  smoothing  formula  is  given 
in  APPLIED  ANALYSIS  pages  316-320.  The  reasoning  used  for 
five  points  and  not  a  higher  number  is  that  these  LWC  devices 
give  "INSTANT"  readings  and  as  such  they  do  give  verifying 
answers  from  second  to  second.  A  longer  smoothing  interval 
would  hide  the  data  behavior  while  a  short  period  would  not 
sufficiently  smooth  the  sample  "noise". 

For  any  given  point  x^  the  smoothing  formula  is : 

Smooth  Ui)=l:‘3£(si-2)+12f<xi‘1,+17f<Jti>+12f<*i+1>-3f(xi+2]/35 
for  the  first  two  points  f(x1),  f(x2): 

£ smooth  <xl>  "  * (Xl)3/35 
^smooth  (x2>  “  f(x2)-4,f(x1)2/5-i*f(x1)l/7 


where: 


A'ffXj)  =  f  (x4)-3f  (x3)  +  3f  (x2)-f  (X]L) 

A'ffXj)  =  f  (x5)-4f  (x4)+6f  (x3)=4f  (x2)+f  (X;i) 
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Program  ICEDERIV  (cont'd) 


for  the  last  two  points  f  (x^j , 

£  smooth ^xxmax)  *  f  <xxmax»  ^xmax' 3/35 

f smooth (xxmax-l>  '  f  (j!xmaX-l,+4’£  txxmax>  2/5-4‘f  (xxmax> 
where: 

£<xxmax>  “  f(x»nax,-3f(xXmai[-l)+3f<xxmax-2,-£<x)!max-l) 
£(xxmax>  =  £<smax,'4£(Bmax-l>+6ftxxmax-2>-,£(xxmax-3)+£(xmax- 

All  the  above  equations  in  addition  to  the  liquid  water 
calculating  module  (for  the  JW-LWC  and  AXIAL  probe  devices) 
from  KNOLL1D  were  encoded  into  ICEDERIV. 

An  additional  feature  of  the  program  is  the  calculation 
of  total  pass  LWC.  Subroutine  SIMPSON  calculates  the  total 
water  content  measured  during  every  icing  cycle.  SIMPSON 
uses  two  separate  methods.  SIMPSON's  and  the  trapezoidal 
rules  of  numerical  integration  are  used  (see:  "Introduction 
to  Numerical  Analysis: by  F.  B.  Hildebrand,  pages  91-95). 

Totals  by  both  methods  are  calculated  for  the  following:  Jtf-LWC 
JW-LWC  smoothed  data,  PMSlD  LWC,  and  PMSlD  LWC  smoothed  data. 
Answers  are  presented  at  the  end  of  the  output  for  each  cycle. 
Operating  instructions  and  sample  output  follow. 


4. 3. 3.1  Program  ICEDERIV  operating  instructions  s.  Sample  Output 
CONTROL  CARDS 

DPSI,CM100000,T400,NT1.  ID#  ID  NAME 

ATTACH  ,  CRT ,  CRTPLOTS  ,MR=1 ,  SN=SHARED . 

LIBRARY , CRT . 

REQUEST , TAPE  3  9 ,  *Q . 

*  VSN ,  TAPE  1  =LY CXXX/NT . 

REQUEST ,  TAPE1 ,  NT ,  E ,  NORING  . 

LGO,PL=99999999. 

EXIT (U) 

DISPOSE, TAPE 3 9, FM. 

REWIND ,  TAPE 2 ,  TAPE 3 ,  TAPE 4  . 

COPY, TAPE 2 
COPY, TAPE 4 
COPY, TAPE 3 
7/8/9 
6/7/8/9 

*  OUTPUT  FILE  FROM  PROGRAM  ICEEX 

DATA  CARDS  (AS  MANY  AS  NEEDED  IN  TIME  ORDER) 

CC  DESCRIPTION 

1-2  MIN-PMS  PROBE  MINIMUM  CHANNEL  #  TO  USE 

(12  FORMAT  -  DEFAULT  ZERO) 

.  3-4  MAX-PMS  PROBE  MAXIMUM  CHANNEL  #  TO  USE 

(12  FORMAT  -  DEFAULT  15) 

*  11-16  START  -  PASS  START  TIME  IN  FORM  HHMMSS 

*  18-23  STOP  -  PASS  STOP  TIME  IN  FORM  HHMMSS 

*  31-35  FCYC  -  FIRST  PASS  CYCLE  (15  FORMAT) 

*  36-40  LCYC  -  LAST  PASS  CYCLE  (15  FORMAT) 

46-50  PLIT  -  PASS  IDENTIFIER  (A5  FORMAT) 

*  either  or  both  pass  times  and  cycle  numbers  may  be  used 


4. 3. 3.1  Program  ICEDERIV  Sample  Outputs  (cont'd) 


A.  Standard  Output  File  for  every  cycle 

1.  second  by  second  listing  of  the  following  parameters 
ICE  COUNT 

2  PT  RATE 

*  JW-LWC 

*  JW-LWC  (SMOOTH) 

*  PMS1D-LWC 

*  PMS1D-LWC  (SMOOTH) 

PMS  MVD 

DEWPOINT  DEPRESSION  (TRUE  TEMP  -  DEWP) 

PRESSURE 
ALTITUDE 
TRUE  TEMPERATURE 
TOTAL  TEMP 

2.  Simpson's  and  trapezoidal  rule  calculation  for  each 
of  the  four  LWC’s  above  (*) 

3.  mean  values  of  DEWPOINT  DEPRESSION  through  TOTAL  TEMP 
in  the  list  above.  Also  the  DEWPOINT  DEPRESSION  and 
TRUE  TEMP  standard  deviations. 


4. 3. 3.1  Program  1CEDERIV  sample  outputs  (cont'd) 

4.  Calculated  least  square  fit  lines  and  correlations 
of  the  following  pairs: 

JW-LWC  -  PMS-LWC 

RATE  -  JW-LWC 

RATE  -  JW-LWC  (SMOOTH) 

RATE  -  PMS1D-LWC 

RATE  -  PMS1D-LWC  (SMOOTH) 

B.  TAPE 2  Output  File 

Cycle  summary  data  file  containing  one  line  of  data  for 
each  cycle  above. 

C.  TAPE 3  Output  File 

Pass  summary  data  file  containing  averages  for  every 
cycle  within  a  pass. 

D.  TAPE 4  Output  File 

Contains  a  one  line  summary  of  least  square  fit  line 
coefficients  for  every  cycle 

NOTE:  B  and  C  above  contain  most  of  the  parameters  found  in 
A-l 

Additionally  ICEDERIV  produces  graphic  output  on  105mm 
microfiche.  The  following  plots  are  produced. 

A.  FOR  EVERY  CYCLE  OF  A  PASS 

1.  Log  normalized  number  density  vs  channel  size 

2.  Scatter  plot  of  PMS  lD-LWC  (SMOOTH)  vs  channel  size 

3.  Scatter  plot  of  JW-LWC  (SMOOTH)  vs  channel  size 


4.3. 3.1  Program  1CEDERIV  sample  outputs  (cont’d) 

4.  time  plot  containing:  rate,  PMS  1D-LWC  (SMOOTH),  and 
JK-LWC  (SMOOTH) 

5.  rate  vs  PMS1D-LWC  (SMOOTH)  scatter  plot 

6.  rate  vs  JW-LWC  (SMOOTH)  scatter  plot 

B.  FOR  EVERY  PASS 

1.  scatter  plot  of  JW-LWC  vs  PMS1D-LWC 

2.  Log  normalized  number  density  vs  channel  size 
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4.3.4  Programs  PLTEXTRACT  &  PLTEXTRACT2 

Program  PLTEXTRACT  produces  pen  plots  which  contain  icing 
data  on  the  top  half  of  the  paper  and  two  lines  of  VCO  plots 
on  the  bottom  half.  The  VCO  data  is  triggered  by  icing  cycles 
VCO  buffers  are  filled  only  when  in  an  icing  cycle.  The  buf¬ 
fers  are  then  plotted  and  flushed  when  one-half  hours'  worth 
of  data  has  been  stored. 

PLTEXTRACT  was  written  as  a  means  of  verifying  the  extrac 
tion  program  (see  ICEEX) .  PLTEXTRACT  plots  the  extracted 
data  on  the  same  scale  as  PLTICE.  This  enables  the  PLTEXTRACT 
to  be  overlayed  on  the  PLTICE  graph  and  visually  check  t  ie 
extraction  process.  Processing  was  performed  on  the  E79-49 
06  DEC  79  flight  (the  benchmark  data  base)  with  no  difference. 

The  usefulness  of  the  PLTEXTRACT  plot  is  more  than  simply 
data  verification.  VCO  information  or  various  LWC  parameters 
can  be  plotted  at  the  bottom  and  LYC  scientists  can  use  this 
information  to  help  them  derive  empirical  relationships . 

In  addition,  PLTEXTRACT2  was  written  to  disregard  the 
presence  or  absence  of  icing  cycles  and  only  plot  two  VCO 
parameters.  This  is  an  important  tool  which  allows  the  JW 
device  adjustments,  made  by  LYC  scientists  to  be  checked. 
PLTEXTRACT2  uses  the  same  KN0LL1D  output  tape  as  input  as 
PLTEXTRACT. 

Operating  instructions  and  sample  output  are  provided 
on  the  following  pages. 


4.3,4. T  PLTEXTRACT  Operating  Instructions 


CONTROL  CARDS 

DPSI,CM130000,T500,NT2.  ID#  ID  NAME 

ATTACH,  LGO,PLEXTRACTBIN,  ID=KAPLANF  ,MR=1 . 

ATTACH  ,  PEN ,  NEWOFFPEN . 

LIBRARY, PEN. 

REQUEST, TAPE39 ,  *Q. 

VSN,TAPE1=TAPE  NO/NT.  (ICEEX  PRODUCED  OUTPUT  TAPE) 
REQUEST  ,  TAPE1 ,  NT ,  NORING ,  E . 

VSN,TAPE2*TAPE  NO/NT.  (KNOLL ID  PRODUCED  OUTPUT  TAPE) 
REQUEST  ,  TAPE2 ,  NT ,  NORING  ,  E . 

LGO. 

7/8/9 

DATA  CARDS 

CARD  1 

COL  1-10  FLIGHT  ID  (FLT  EXX-XX) 

11-20  FLIGHT  DATA  (XXXXXXX) 

80  set  to  1  when  tape2  available 

CARD  2 

COL  1-12  PLOT  #1  LITERAL 

CARD  3 

COL  1-12  PLOT  #2  LITERAL 

CARD  4  (FREE  FORMAT  IN  FOLLOWING  ORDER) 

PLOT  #1  ID  NUMBER 

v  PLOT  <2  ID  NUMBER 

PLOT  #1  MINIMUM  VALUE 
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Plot  #1  MAXIMUM  VALUE 

Plot  #2  MINIMUM  VALUE 

Plot  #2  MAXIMUM  VALUE 

6/7/8/9 

When  tape  2  (KNOLL ID  data)  is  not  available  change  NT2  to 
NTl  on  the  job  card. 

Remove  the  tape  request  cards  for  tape  2  and  set  COL  80 
to  zero  of  data  card  1. 
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D-A133  557  DEVELOPMENT  OF  THEORETICRL  MODELS  AND  ANALYSES  OF  THE 
MICRO-PHYSICS  OF  CL..  (U>  DIGITAL  PROGRAMMING  SERVICES 
INC  WALTHAM  MA  L  E  BELSKV  ET  AL.  31  MOV  82 
UNCLASSIFIED  AFGL-TR-83-0123  F19628-81-C-0141  F/G  4/2 


2/3 

NL 


F/G  4/2 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANOARDS-1963-A 


PLOT  ID  CODES 

AVAILABILITY* 


NUMBER 

DESCRIPTION 

TAPE  1 

TAPE 

1 

TRUE  TEMPERATURE 

YES 

YES 

2 

TOTAL  TEMPERATURE 

YES 

NO 

3 

DEWPOINT/FROSTPOINT 

YES 

YES 

4 

JW-LWC 

YES 

YES 

5 

PRESSURE 

YES 

YES 

6 

TRUE  AIRSPEED 

YES 

YES 

7 

HEIGHT 

YES 

YES 

8 

ASSP  PROBE  LWC 

NO 

YES 

9 

CLOUD  PROBE  LWC 

NO 

YES 

10 

PRECIPITATION  PROBE  LWC 

NO 

YES 

11 

TOTAL  PROBE  LWC 

NO 

YES 

12 

LOG  ASSP  PROBE  NT 

NO 

YES 

13 

LOG  CLOUD  PROBE  NT 

NO 

YES 

14 

LOG  PRECIPITATION  PROBE  NT 

NO 

YES 

15 

LOG  TOTAL  PROBE  NT 

NO 

YES 

16 

ASSP  PROBE  DO 

NO 

YES 

17 

CLOUD  PROBE  DO 

NO 

YES 

18 

PRECIP  PROBE  DO 

NO 

YES 

19 

TOTAL  PROBE  DO 

NO 

YES 

20 

ASSP  PROBE  DO  TIMES  LWC 

NO 

YES 

21 

CLOUD  PROBE  DO  TIMES  LWC 

NO 

YES 

22 

PRECIP  PROBE  DO  TIMES  LWC 

NO 

YES 

23 

TOTAL  PROBE  DO  TIMES  LWC 

NO 

YES 

TAPE 2  (KNOLL ID  OUTPUT  TAPE)  IS  USED  IF  COL  80  OF  ID  CARD 
IS  SET  -  ELSE  DEFAULT  IS  TAPEl  (ICEEX  OUTPUT  TAPE)  WHICH 
LIMITS  AVAILABILITY  TO  FIRST  SEVEN  OPTIONS. 


4. 3. 4. 2  PLTEXTRACT 2  Operating  Instructions  An3  Output 

CONTROL  CAROS 

CLASS, CN130000 ,T20D ,NT1 .  ID  i  ID  NAME 

ATTACH ,  LCO  ,  PLTEXTRACT2B  I N ,  1D-KAPLANF f  MR:=1  _ 

ATTACH , PEN , OF FL INEPEN , SN-SHARED . 

LIBRARY, PEN. 

REQUEST, TAPE39.*Q . 

•VSN.TAFE2.TAPENO/NT  ( KNOLL  ID  PRODUCED  OUTPUT  TAPE) 

•  REQUEST, TAPE2, NT, NOR  INC, E . 

LCO. 

7/8/f 

DATA  CARDS 
CARD  I 

COL  1-10  FL1CHT  ID  <FLT  EXX-XX  ) 

11-20  FLIGHT  DATE  (XX-XXX-XX) 

CARD  2 

COL  1-12  PLOT  81  LITERAL 
CARD  3 

COL  1-12  PLOT  82  LITERAL 
CARD  4 

COL  <  FREE  FORMAT) 

PLOT  81  ID  NUMBER 

PLOT  82  ID  NUMBER 

PLOT  81  MININUN  VALUE 

PLOT  81  MAXIMUM  VALUE 

PLOT  82  MINIMUM  VALUE 

PLOT  82  MAXIMUM  VALUE 

4/7/B/9 


4. 3. 4. 2  PLTEXTRACT2  Operating  instructions  And  Output  (cont'd) 


plot _ ui _ cam 

number  PtSCB-IPTiaN 

1  TRUE  TEMPERATURE 

2  NOT  USED 

3  DEWPOINT/ FROSTPOINT 

4  JV-LVC 

5  PRESSURE 

4  TRUE  AIRSPEED 

?  HEICHT 

0  ASSP  PROBE  LVC 

9  CLOUD  PROBE  LVC 

10  PRECIPITATION  PROBE  LVC 

11  TOTAL  PROBE  LVC 

12  LOG  ASSP  PROBE  NT 

13  LOG  CLOUD  PROBE  NT 

14  LOG  PRECIPITATION  PROBE  NT 

15  LOG  TOTAL  PROBE  NT 

14  ASSP  PROBE  DO 

17  CLOUD  PROBE  DO 

10  PRECIP  PROBE  DO 

17  TOTAL  PROBE  DO 

20  ASSP  PROBE  DO  TIMES  LVC 

21  CLOUD  PROBE  DO  TIMES  LVC 

22  PRECIP  PROBE  DO  TIMES  LVC 

23  TOTAL  PROBE  DO  TIMES  LVC 


NOTE  : 


TAPE2  (KNOLL1D  Output  Ttpi)  is  used.  Column  SO 
of  ID  etrd  must  bs  sst . 


4.4  PMS-2D  Data  Processing 

DPSI  designed  and  developed  a  set  of  programs  to  analyze 
the  data  obtained  from  the  PMS-2D  Knollenberg  devices.  The 
details  of  this  design  were  described  in  reports  from  pre¬ 
vious  contracts.  A  new  program  for  automatic  particle  typing 
developed  during  this  contract  is  described  in  the  next  sec¬ 
tion. 

The  2D  particle  display  system  is  a  valuable  source  of 
data  to  LYC.  Data  was  collected  from  two  independent  systems 
(one  each  on  the  C130E  and  Learjet)  and  recorded  on  magnetic 
tape  using  a  nine-track  Pertec  recorder.  Only  the  C130  tapes 
are  processed  on  the  AFGL  CDC  6600  computer.  PMS-2D  data  were 
collected  aboard  a  Convair  NASA  990  aircraft  and  those  data 
were  processed  via  the  2D  programs  described  below. 

These  data  tapes  contain  two  types  of  records,  composed 
of  "fast"  (long)  and  "slow"  data.  The  fast  data  records 
contain  the  2D  particle  image  slices.  One  dimension,  the 
columns,  is  represented  by  the  32  diode  array.  The  other 
dimension,  the  rows,  is  achieved  by  taking  readings  over 
time  so  that  a  two  dimensional  picture  of  the  shadow  is  ob¬ 
tained.  The  speed  of  the  aircraft  is,  theoretically,  kept 
at  a  predetermined  constant  speed  based  on  the  resolution 
diameter)  of  the  2D  probes  and  the  frequency  of  the  clock 
polling  the  sensors.  That  is,  the  C130  aircraft  flew  at  ap¬ 
proximately  200  knots  in  order  for  the  25y  diameter  diode  sen¬ 
sor  to  produce  a  square  image  of  each  particle  at  the  given 
sampling  rate  of  4  MHz.  Likewise,  the  2D  system  installed  a- 
board  the  NASA  Convair  990  was  operated  at  true  airspeeds  of 


4.4  PMS-2D  Data  Processing  (cont'd) 


of  approximately  466  knots  in  order  for  the  60y  diode  sensor 
to  obtain  a  square  image  at  the  same  4  MHz  sampling  rate. 
However,  because  true  air  speed  of  the  aircraft  cannot  be 
kept  constant,  the  data  are  modified  slightly  by  the  2D  data 
processing  programs  to  maintain  the  square  image  representa¬ 
tion  of  the  particles.  Thus  the  2D  device  gives  a  picture  of 
the  particle (s)  as  subsequent  readings  are  placed  together. 

The  slow  data  records  are  written  to  tape  once  every 
ten  seconds.  These  records  contain  VCO  and  analog  information. 
In  addition,  selected  ID  data  is  multiplexed  in  the  2d  buffer 
and  also  recorded  in  the  slow  data  records.  The  exact  infor¬ 
mation  contained  in  these  records  is  different  for  each  air¬ 
craft. 
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4.4.1  Program  PMS2D 


Program  PMS2D  takes  the  2D  automatic  typing  algorithm  de¬ 
veloped  by  Dr.  Hunter  of  ADAPT,  Inc.  and  packages  them  in  a 
manner  which  allows  maximum  use  and  ease  of  operation  by  LYC 
personnel . 

These  routines  take  the  2D  scans  and  converts  them  to 
normalized  integer  values  called  "hols".  Other  routines  in¬ 
terpret  these  hols  by  converting  them  to  "ANT"  vectors  for 
use  by  eigenvector  analysis.  The  routines  are  copyrighted 
by  Dr.  Hunter  and  will  not  be  described  here.  Essentially 
what  DPSI  has  done  is  to  write  driver  routines  to  process 
2D  data  tapes.  As  originally  delivered  by  Dr.  Hunter  the 
routines  could  only  process  at  most  300  particles  per  comp¬ 
uter  run.  The  new  program  has  been  written  to  allow  process¬ 
ing  of  300  particles  per  data  record.  Experience  has  shown 
that  there  are  only  approximately  a  dozen  acceptable  particles 
per  data  record  hence  no  particles  should  be  lost  during  the 
processing  of  a  run.  Operating  instructions  and  sample  output 
are  provided  on  the  following  pages. 
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4.4.1  PROGRAM  PMS2D  Operating  Instructions  &  Sample  Output 


COMMAND  DECK 


COMMAND  DECK 

MlLLR,CM115000tNT1,T500.  141 

VSN,TAPE1=KNE370/NT. 

REQUEST ,TAPE 1  ,PE  ,L  ,NR , NT . 

ATTACH,P,PMS2D, 1D*MILLERP,MR=1 . 

FTN5. 1 *P ,B*LG0 ,L=0 ,PL*999999 . 

EILE(TAPE1 tRT«U,BT*KfHRL»5574fM8L*5576fRBls1 ,BFS*560) 
LDSET.FILES-TAPE1 ,PRESET*ZERO. 

LGO. 

EXIT (U) 

7/8/9 


MILLERP 


INPUT  DATA 

NAMELIST/IOPTS/ 

NAMELIST/7YPDA7/ 

ID  CARD 

1-10  FIT  ID  (FIT  E80-16) 
11-20  FIT  DATE(DDHONYR) 
21-30  KNE  TAPE  MMKNE370) 
PASS  TINES(MAXIMUM  20) 

1-6  START (HHMMSS) 

7-10  STOP (HHMMSS) 

7/8/9 
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4.5  Special  Analysis  and  Processing  Techniques 
4.5.1  PROGRAM  CNVRT 

Program  CNVRT  is  a  data  processing  program  written  to 
aid  LYC  scientists  in  the  analysis  of  data  from  NCAR  CLOUD 
DATA  TAPES.  These  tapes  are  recorded  with  524  60-bit  words 
per  record.  Each  word  contains  3  20-bit  integers.  These 
parameters  are  converted  according  to  a  given  table  and  dis¬ 
played  in  a  format  suitable  to  LYC  scientists. 

Operating  instructions  and  sample  output  are  provided 
on  the  following  pages. 


4. 5. 1.1  Operating  Instructions  for  Program  CNVRT 


MILLER, T100,NT2.  1416 

ATTACH , I , CNVRT , ID=KAPLANF . 

ATTACH , TAPES , IS=KAPLANF . 

VSN , TAPE3=LYCB0B/NT . 

REQUEST , TAPE 3 , NT , NORING , L , PE . 

FTN, I=I,PL=99999. 

REWIND, TAPE 3. 

LGO. 

EXIT (U) 

DISPOSE, OUTPUT, PR, *C . 


MILLERP 


NCAR  ELECTRA  STRATEX-I  17JUNT6  FIRST  TIME  ON  THIS  FILE  11  15  55  THIS  FILE  I 

S  ALL  OR  PART  OF  TU'E  PERIOD  H  15  55  TO  16  39  40  DESCRIPTION  OF  RECORD  —  115  PARAMETERS  WERE  SAVE 
0  AT  THEIR  RESPECTIVE  RATES.  THIS  REPRESENTS  786  SAStPLES/PROGRAM  CVCLE  WHERE  A  CTCLE  IS  1.000  SEC 
THE  2  CYCLES  OF  786  SAMP/CYC  .  1572  WERE  THEN  SCALED  INTO  20  BIT  INTEGERS  AND  PACKED  3  SAMPLCS/WOR 
0  INTO  524  60  BIT  PACKED  WORDS  —METHOD  OF  SCALING —  A  9IAS  AO II)  WAS  ADDED  TO  EACH  SAMPLE 
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4. 5. 1.1  PROGRAM  CNVRT  Sample  Output 
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4.5.2  PROGRAM  STEVERD 


LYC  scientists  are  continually  searching  for  new  sources 
of  cloud  physics  or  icing  data.  In  early  1983  LYC  began  re¬ 
ceiving  icing  data  from  West  German  scientists.  This  data 
is  formatted  so  that  the  first  two  records  contain  flight  in¬ 
formation  such  as  date,  time,  total  records  and  comments. 

During  each  flight  eleven  parameters  were  collected  and 
written  to  the  tape  (records  3-13,14-24...).  Each  record 
contains  600  words  and  pertains  to  one  parameter.  Each  word 
in  the  record  represents  data  collected  over  a  one  second 
interval. 

This  program  reads  and  formats  these  data  for  LYC 
scientists. 

The  order  of  these  parameters  on  the  tape  is: 
time,  altitude,  air  speed,  static  temperature,  dewpoint,  dif¬ 
ference,  humidity  potential,  LWC,  backscattering  coefficient, 
degree  of  icing,  icing  rate,  and  ice  warning  signal. 

The  German  data  tapes  are  processed  via  the  CDC  job 
shown  below: 

DATA  TAPE  — >  SYSTEM  UTILITY  COPYCF  — >  OUTPUT  DATA  FILE  - 

PROGRAM  BUFFTNS  — >  OUTPUT  DATA  FILE  — >  PROGRAM  STEVERD  - 
FORMATTED  LISTING 

NOTE:  Program  COPYCF  is  batched  via  a  card  deck 

Program  BUFFTNS  and  STEVERD  are  run  interactively 

Operating  instructions  and  sample  output  for  these  programs 
follow. 


4. 5. 2.1  PROGRAM  STEVERD  OPERATING  INSTRUCTIONS 


COPYCF  UTILITY  JOB 

MILLR,NT1.  1416  MILLERP 

REQUEST, GDATA16 ,*PF. 

VSN , TAPE 1=LYCGRM/NT . 

REQUEST , TAPEl , PE , NT , EB , NORING , L . 

FILE (TAPEl, RT=U,BT=K,MRL=10000,MBL=10000,RB=1,BFS=1005) 
COPYCF ( TAPEl , GDATA1 6 ) 

CATALOG , GDATA16 , , ID=KAPLANF . 

7/8/9 

6/7/8/9 


PROGRAM  BUFFTNS 

ATTACH , B , BUFFTN5, ID=KAPLANF 
ATTACH , TAPE  3 , GDATAl 6 , ID=KAPLANF 
FTN,I=B,PL=999999 
LGO 

PROGRAM  STEVERD 

ATTACH , S , STEVERD , ID=KAPLANF 

FTN5,I=S,PL=999999 

LGO 

DISPOSE, OUTPUT, PR, 191 
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4.5.3  Program  ATTCHFLES 


A  procedure  was  written  that  eliminates  the  need  to  in¬ 
dividually  attach  files  every  two  weeks  so  that  they  will 
not  expire.  Every  time  a  new  program  is  written  or  an  old 
program  becomes  obsolete  it  can  be  added  or  deleted  by  using 
this  procedure.  In  addition,  this  procedure  can  be  easily 
modified  so  that  it  will  also  catalog  all  files  to  the  disk. 

Operating  instructions  and  a  listing  of  the  program  may 
be  found  in  the  next  section. 


4. 5. 3.1  PROGRAM  ATTCHFLES  (cont'd) 


TOP  OF  FILE 

MILLR/T100.  1416 

REQUEST, A, *PF. 

ATTACH , B , BATCHPLT , I D=KAPLANF . 

REWIND ,  A ,  B . 

COPY, A, B. 

UNLOAD, A, B. 

REQUEST, A, *PF. 

ATTACH , B , BUFFTN5 , ID=KAPLANF . 

REWIND , A , B . 

COPY, A, B. 

UNLOAD, A, B. 

REQUEST, A, *PF. 

ATTACH , B , CNVRT , ID-KAPLANF , 

REWIND, A, B. 

COPY, B, A. 

UNLOAD, A, B. 

REQUEST, A, *PF. 

ATTACH , B , COPYLYCBOB , ID=KAPLANF . 

REWIND, A, B. 

COPY, B, A. 

UNLOAD, A, B. 

REQUEST, A, *PF . 

ATTACH , B , ICEDERIVBIN , ID=KAPLANF . 

REWIND, A, B. 

COPY, B, A. 

UNLOAD, A, B. 

REQUEST,A,*PF. 

ATTACH ,  B ,  ICEDMPBIN ,  ID«=KAPLANF . 

REWIND, A, B. 

COPY, B, A. 

UNLOAD, A, B. 


MILLERP 


4.5.4  PROGRAM  KNOLL1D  AND  HDATSTAT 


In  order  to  assist  in  the  analysis  of  icing  data  collected 
during  the  flight  of  11  DEC  79  (E79-50) ,  DPSI  has  prepared  a 
modified  version  (cycle  18)  of  KNOLLlD  to  produce  an  output 
file  (tape  11)  for  additional  processing.  This  file  contains 
scatter  probe  counts  for  channels  1-15  (one  channel  per  re¬ 
cord)  ,  with  VCO  data  (temp,  dewpoint,  and  relative  humidity) 
contained  in  the  16th  record.  This  sequence  is  repeated  for 
each  second  of  the  pass. 

Subsequently,  program  HDATSTAT  was  written  to  process 
these  scatter  probe  and  VCO  data.  The  program  produces  two 
types  of  output:  first,  the  total  count  for  channels  1-15  ar •» 
provided  for  each  second  of  the  pass;  and  then  the  average 
value  for  each  channel  is  completed  for  the  pass  (standard 
deviation  and  variance  are  also  completed  for  each  channel) . 

In  addition,  the  average  values  for  temperature,  dewpoint, 
and  relative  humidity  for  the  pass  are  calculated  after  the 
channel  1  summary  is  printed.  Since  these  VCO  values  occur 
only  once  per  second,  repetition  of  these  averaged  values  after 
each  subsequent  channel  summary  would  be  meaningless. 

Operating  instructions  for  KNOLLlD  (cycle  18)  and  oper¬ 
ating  instructions  and  sample  output  for  HDATSTAT  are  provided 
in  the  next  section. 


4. 5.4.1  PROGRAM  KNOLL1D  Operating  Instructions 


COMMAND  DECK 

MILLR,CM65000,TPl,T100.  1416  MILLERP 

ATTACH , LGO , KNOLL 1DB IN , I D= KAPLAN, MR= 1 , CY= 1 8 . 

ATTACH, TAPE 8 ,VCOCALS , ID=KAPLANF ,MR=1 . 

REQUEST, TAPE11 ,*PF. 

VSN , TAPE1=PMS  317. 

REQUEST , TAPEl , S ,HI ,  MT . 

FILE(TAPEl,RT=U,BT=K,MRL=1024 ,RB=1 ,BFS=105) 

LDSET , PRESET=  ZERO , F I LES=TAPE 1 . 

LGO. 

EXIT(U) 

REWIND, TAPE 3. 

REWIND, TAPEl 1. 

COPY,TAPE3. 

CATALOG , TAPEl 1, HDATA43TD , ID=MILLERP . 

7/8/9 
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4. 5. 4.1  KN0LL1D  operating  instructions  (cont'd) 


DATA  CARDS 
CARD  1  -  ID  CARD 


COL  1-10  FLT  XYR-NN 

COL  12-20  DD  MON  YR 

COL  21-30  INPUT  TAPE  NUMBER 

COL  31-40  PLOT  TAPE  NUMBER 

COL  52-59  PMS  ZERO  SECONDS  -  HH:MM:SS 

COL  70  IOSKP ( 0=STANDARD  OUTPUT; 1=CREATE  TAPE  11  DATA  FILE) 


CARD  2  -  NAMELIST  SCOEF 


ARRAY  S  ( 5)  OF  SOUNDING  COEFFICIENTS  (DEFAULTED  INTER¬ 
NALLY)  . 


CARD  3  -  NAMELIST  VCOEF  -  USED  TO  OVERRIDE  VCO  CALIBRATIONS 

ARRAY  C(3,13)  CONTAINING  THE  VCO  CALIBRATION  COEFFI¬ 
CIENTS  IN  THE  FOLLOWING  ORDER:  INTERCEPT,  FIRST  DEGREE 
AND  SECOND  DEGREE  COEFFICIENTS. 

VCO'S  ARE  IN  THE  FOLLOWING  ORDER: 

1  INDICATED  AIRSPEED 

2  TEMPERATURE 

3  EWER 

4  TWCI-4 

5  DEWPOINT/FROSTPOINT 

6  LWC-JW 

7  MAGNETIC  HEADING 

8  PRESSURE  KISTLER 

9  TRUE  AIRSPEED 

10  ICEING  RATE 

11  TWCI-1 

12  TWCI-2 

13  TWCI-3 


CARD  4  -  NAMELIST  JWADJ 

CONTAINS  HEIGHT  PROFILES  FOR  A  JW-LWC  ADJUSTMENT. 

ELEMENTS  ARE: 

L  NUMER  OF  LEVELS  (DEFAULT  0-NO  CORRECTION) (MAX  10) 

HT  HEIGHT  OF  LEVELS  IN  KM'S  (HT(1)  GT  HT(2)  GT...GT  HT (L) ) 
XA  ORIGIN  OF  LEVEL  (ONE  PER  LEVEL) 

SLA  SLOPE  FROM  LEVEL  (I)  TO  LEVEL  (1  +  1)  (L-l  SLOPES 
REQUIRED 
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4. 5.4.1  KNOLLlD  operating  instructions  (cont’d) 


CARD  5  -  OPTION  CARD  (ALL  15  FORMAT) 

COL  1-5  NFHEADR  -  SET  TO  1  WHEN  CONDENSED  FORMAT  HEADER 
CARD  DESIRED 

6-10  ICLK  1  =  USE  A/C  CLOCK 

2  =  USE  PMS  CLOCK 

11-15  IDAT  1  =  USE  FINALIZED  DATA  LITERAL 

0  =  USE  PRELIMINARY  LITERAL 

16-20  IPLT  1  =  PLOT  TAPE  PRODUCED 

0  *=  NO  PLOT  TAPE 

21-25  N2PROBE  2  -  TWO  CLOUD  PROBES 

3  =  TWO  PRECIP  PROBES 

26-30  ITMP  0  =  TEMPERATURE  DETERMINATION  BY  VCO  PROFILE 

1  ■=  TEMPERATURE  DETERMINATION  BY  STANDARD 

ATMOSPHERE 

2  *  TEMPERATURE  DETERMINATION  BY  RADIOSONDE 

PROFILE 

31-35  JVCO  -  NUMBER  OF  VCOFIX  CARDS  DESIRED  (0-10) 

36-40  NSKP  >0  IS  THE  NUMBER  OF  SJWD-OF-FtLES  TO  SKIP 

BEFORE  PROCESSING  DATA 

<0  IS  THE  NUMBER  OF  RECORDS  TO  SKIP  BEFORE 

PROCESSING  THE  DATA 

41-45  IDEK  1  *=  PRODUCE  A  PUNCHED  DECK  CONTAINING  THE 

TIME,  HEIGHT,  LWC  AND  DO  VALUES, 

0  *  NO  DECK 

46-50  MXLINES  -  NUMBER  OF  PASS  LITERALS  DESIRED  FOR  THE 

CONDENSED  OUTPUT  FORMAT  (0-15) 

51-55  INTRP  1  =  USE  INTERPOLATION 

0  ■  NO  INTERPOLATION 

56-60  IDMZ  1  «  SUMMARY  FILE  OF  DO, LWC  AND  Z  VALUES 

PRODUCED 

0  =  NO  SUMMARY  FILE  PRODUCED 

61-65  IVEL  1  =  USE  TRUE  AIRSPEED 

0  «=  USE  CALCULATED  AIRSPEED 

66-70  I  FORM  1  •  CONDENSED  OUTPUT  FORMAT  PRINTED  AND  THE 

PLOT  TAPE  USES  UNMELTED  BARWIDTH’S 
0  =  NORMAL  OUTPUT  AND  NORMAL  PLOT  TAPE 

71-75  I  OUT  1  =  NO  STANDARD  OUTPUT 

0  ■  STANDARD  OUTPUT  PRODUCES 

76-80  ISCAT  1  ■  USE  .9  FACTOR  IS  SCATTER  CALCULATIONS 

0  *  .9  FACTOR  NOT  USED 


4. 5. 4.1  KNOLLlD  operating  instructions  (cont'd) 


CARD  6  NEW  FORMAT  HEADER  CARD  (IF  NFHEADR  .EQ.  1) 

CENTERED  LINE  (A80)  TO  BE  PRINTED  AT  THE  TOP  OF  EACH 
NEW  FORMAT  OUTPUT. 

CARD  7  TYPE  LITERAL  LINE 

IF  MXLINES  >  0  MXLINES  CARDS  ARE  REQUIRED  HERE  IN 
CENTERED  A80  FORMAT.  THESE  LINES  ARE  PRINTED  ON 
NEW  FORMAT  OUTPUT;  BENEATH  THE  INTERVAL  PARTICLE 
TYPE. 

CARD  8  ONWARDS 

A)  ANY  DATA  CARDS  REQUIRED  BY  SWITCHED  SET  ON  OPTION 
CARD  JVCO  0  IMPLIES  VCO  PROFILES  IN  HERE 
ITEMP=2  IMPLIED  RADIOSONDE  PROFILES  HERE 

B)  TYPE,  EDIT,  HTOX,  XTOD  CARDS  INTERSPERED  IN  HERE 

TYPE  (15  MAXIMUM) 

COL  1-4  TYPE 

COL  6-13  HH:MM:SS  (START  TIME) 

COL  16-23  HH:MM:SS  (STOP  TIME) 

COL  25-26  CLOUD  TYPE 

COL  27-28  PRECIP  TYPE 

COL  31-35  AVERAGING  INTERVAL  (15) 

COL  45-64  LEFT  JUSTIFIED  PASS  LITERAL  FOR  HEADER  OF  NEW 
FORMAT  (A20) 

EDIT  (5  MAXIMUM) 

COL  1-4  EDIT 

COL  6-13  HH:MM:SS  (START  TIME) 

COL  16-23  HH :MM: SS  (STOP  TIME) 

COL  26,28,30  PROBES  TO  BE  EDITED  (1,2,3) 

COL  31-54  CHANNELS  TO  BE  EDITED  (13  FORMAT,  8  MAXIMUM) 

HTOX  (NO  MAXIMUM) 

COL  1-4  HTOX 

COL  6-7  TYPE  CODE  (ODD  NUMBER) 

COL  9-10  EQUATION  NUMBER 

COL  12-13  ARGUMENT  TO  BE  CHANGED  (1=M,  2=B,  3=BREAKPT) 

COL  15-30  NEW  VALUE  (F15.0) 

XTOD  (NO  MAXIMUM)  SAME  AS  HTOX 

COL  12-13  ARGUMENT  TO  BE  CHANGED  (l=CO,  2=EX,  3=BREAKPT) 


1  ^  j 


4. 5. 4. 2  PROGRAM  HDATSTAT  Operating  Instructions  and  Sample  Output 
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4.6  Snow  Data  Processing  -  CDC  6600 
4.6.1  Program  CNVRTASCII  to  CDC 

There  have  been  problems  using  the  CDC  built-in  data 
(ASCII)  conversion  software.  DPSI  has  developed  a  program 
which  will  read  an  ASCII  file  and  perform  a  table  look  up 
in  order  to  convert  the  data  the  CDC  display  code. 

Program  CNVRTASCII  to  CDC  incorporates  the  binary  buffer- 
in  I/O  technique  preventing  the  CDC  operating  system  from 
assisting  in  the  conversion.  Three  types  of  output  are  pro¬ 
vided;  formatted  display  code,  octal  bytes,  and  record  by 
record  output  of  the  display  file. 


130 


4.6.1  Program  CNVRTASCII  to  CDC  (cont'd) 


(A) 


(B) 


JOB  NAME,NT1*,T100. 

ATTACH , I , CNVRTASCI ITOCDC , ID=AOBERTSK . 


/**FTN,I*I,SL,R=3. 

4  VSN, TAPE1*  INT* . 

REQUEST, TAPEl,HD*,S , VS , NT* , NORING . 
LGO. 

EXIT(u) 

<  CATALOG, TAPE 2,  , ID* 

CATALOG , TAPE3 ,  , ID* 

CATALOG , TAPE4 ,  , ID* 

6/7/8/9  ” 


aaaa  PPPPPPP 


*  7  or  9  truck,  and  bit  density  may  change  for  each  specific  job. 

(A)  If  the  source  ASCII  data  is  on  disk,  replace  with: 

ATTACH,  TAPE1, _ ,  ID _ . 

(B)  Rearrange  as  required,  for  just  a  printout,  replacement: 

REWIND,  TAPEx. 

COPY , TAPEx , OUTPUT . 
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4.6.2  Program  MAKECRRELTAPE 


Program  MAKECRRELTAPE  is  a  data  conversion  program  which 
will  convert  a  snow  weight,  7 track,  OS/8  tape  generated  by 
the  PDP-8  to  a  9track, ASCII  format  supplied  to  CRREL.  Since 
this  program  uses  tapes,  it  executes  under  batch.  Excluding 
the  media  change,  the  apparent  difference  between  the  input 
and  output  of  the  program  can  be  seen  in  section  4. 6. 3.1. 


The  input  data  tape  containing  384  characters  per  record 


4.6. 2.1  Program  MAKECRRELTAPE  Operating  Instructions 


JOB  NAME,TPl,NTl,T100.  aaaa  ppppppp 

ATTACH, I , MAKECRRELTAPE, ID=ROBERTSK . 

FTN,I=I,SL,R»3. 

VSN,  TAPE1= _ . 

REQUEST , TAPEl , S , MT , NORING . 

FILE ( TAPEl , RT=u , BT=K,MRL=10000 ,MBL=10000 , RB=1 ,BFS=1005) 

VSN,  TAPE  2= _ /NT. 

REQUEST, TAPE 2 , NT, VS , RING . 

LGO. 

7/8/9 

6/7/8/9 


TAPEl=INPUT=PDP-8  OS/8  7 TRACK  TAPE 
TAPE 2=OUTPUT=CRREL9 TRACK  ASCII  TAPE 


4.6.3  Program  MAKE1MINSN0RATE 


The  techniques  employed  for  snow  fall  analysis  requires 
compatible  data  sets  with  respect  to  time.  The  CRREL  ASCME 
data  is  recorded  as  a  one  minute  average.  Therefore  to  pre¬ 
sent  an  accurate  representation  of  fall  velocity,  it  becomes 
necessary  to  calculate  AFGL  snow  rates  as  one  minute  averages. 

Program  MAKE  1MINSN0 RATE  has  been  designed  to  provide  one 
minute  averaged  snow  rates  and  an  instantaneous  value  (snow 
rate  as  taken  from  the  nearest  whole  minute) .  This  program 
operates  under  INTERCOM  and  accepts  as  input  an  AFGL  format 
#2  snow  weight  data  file  (see  section  4 . 6.  3.  l)Output  consists 
of  three  data  files;  averaged  snow  tape,  instantaneous  snow 
rate,  and  snow  rates  for  each  2.809  second  period. 
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4. 6. 3.1  Program  MAKE 1MINSN0RATE  Operating  Instructions 

ATTACH , I , MAKE 1MINSNORATE , ID=ROBERTSK . 

FTN ,1=1, SL , R=  3 . 

ETL,100. 

ATTACH, TAPEl,  _ *,ID= _ . 


TAPE1=INPUT=AFGL  SNOW  WEIGHT  DATA* 

TAPE  2=OUTPUT=AVE RAGED  1  MINUTE  SNOW  RATE  ** 

TAPE 3=OUTPUT= INSTANTANEOUS  1  MINUTE  SNOW  RATE  ** 
TAPE4 =OUTPUT=SNOW  RATE  FOR  EACH  2.809  SECONDS 


AFGL  DATA  IN  FORMAT  #2 


**  AFGL  DATA  IN  FORMAT  #3 


AFGL  FORMAT  #2 


COL 

30 

+ 


COL 
1 
+ 

HH  MM  SS  DATE 
+  XXX. XX  +  XXX. XX 
+  XXX. XX  +  XXX. XX  +  XXX. XX 

DATA  SUPPLIED  TO  CRREL  ON  A  9  TRACK  ASCII  TAPE,  RECORD  LENGTH, 
3-60  BIT  WORDS 


1  17 

MM/YY  SS  . xxxxxxx 


28  39  50  57 

+  i  4-  4- 

. XXXXXXX  . xxxxxxx  . xxxxxxx 

X  X  X  X 

An  X  in  the  decimal  point  columns  indicates  that  the  ASCME 
device  was  imperative;  or  out  for  maintenance.  When  CRREL 
tapes  are  processed  initially,  the  X's  are  replaced  with 
'-1.0000000'.  This  will  simplify  the  post  processing  stream 
programs. 


CRREL  DATA  FORMAT 
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4.6.4  Program  COPYASCME 

At  times  CRREL  (Cold  Regions  Research  Laboratory)  supplies 
LYC  with  data  tapes  containing  field  observation  values  from 
their  ASCME  (Airborne  Snow  Mass  Concentration  Measurements) 
device.  The  ASCME  data  is  used  at  LYC  for  comparison  with 
LYC  snow  rates  and  for  fall  velocity  analysis. 

In  order  to  process  these  data  tapes,  they  must  be  read 
and  converted  to  CDC  display  code.  Program  COPYASCME  is  a 
quick  and  simple  program  which  will  read  each  file  (usually 
there  is  more  than  one  file)  from  a  CRREL  tape  and  allow  them 
to  be  catalogued. 

The  CRREL  data  tapes  have  been  formatted  on  a  9  track 
magnetic  tape  written  in  ASCII  at  800  BPI  (Bits  Per  Inch) . 

This  format  could  change,  but  only  the  "REQUEST"  job  card 
would  then  require  modification. 

Presently,  COPYASCME  will  extract  five  files  from  the 
CRREL  tape,  if  more  files  are  required,  change  the  "PROGRAM" 
statement  within  the  COPYASCME  program.  Only  one  data  card 
is  necessary  when  operating  this  program;  it  defines  the 
number  of  files  to  transfer. 


4.6.4. 1  Program  COPYASCME  Operating  Instructions 


JOB  NAME,NT1*,T100. 

ATTACH , I , COPYASCME , ID=ROBERTSK . 

FTN , 1=1 , SL, R=3 . 

VSN ,  TAPEl= _ /NT. 

REQUEST , TAPE 1 ,  HD* , S , VS* , NT* , NORING . 


LGO. 

EXIT (u) 

CATALOG ,  TAPE  2 , _ ,  ID= 

CATALOG ,  TAPE3 , _ ,  ID= 

CATALOG ,  TAPE4 , _ ,  ID= 

CATALOG ,  TAPE  5 , _ ,  ID= 

CATALOG ,  TAPE  6  f _ ,  ID= 

7/8/9 


xxx  *■  number  of  files  to  copy, column  1-3 
6/7/8/9 


ppppppp 


7  or  9  track,  bit  density,  and  recording  format  could  be  changed 
at  any  time. 


TAPE1=INPUT=CRREL  DATA  TAPE  (see  Appendix  2) 
TAPE  2) 

;  \foutput-cdc  formatted  files 

TAPE6J 
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4.6.5  Programs  CORRSNOWRATEDATA/CORRSNOLOGSCALE 

In  order  to  complete  analysis  of  the  snow  data  collected 
by  both  LYC  and  CRREL,  DPSI  has  developed  two  data  correlation 
programs;  CORRSNOWRATEDATA  and  CORRSNOLOG SCALE .  Both  programs 
operate  in  an  identical  fashion,  however,  the  former  generates 
a  linear  interpretation  of  the  data,  whereas  the  latter  is  a 
logarithmically  based  comparision. 

Both  programs  have  been  designed  to  generate  Tektronix 
CRT  plots  in  an  interactive  environment  (INTERCOM) .  However, 
with  a  few  modifications,  output  can  easily  be  sent  to  any 
of  the  plotting  devices  available. 

The  ASCME  tapes  contain  data  for  three  devices,  and  thus 
three  plots  are  created.  At  the  conclusion  of  each  plot  "FRAME" 
is  called  allowing  the  operator  time  to  produce  a  hardcopy  of 
the  plotted  image. 

It  is  suggested  that  the  19  inch  Tektronix  be  used  for 
final  copies  as  the  resolution  is  preferable  to  the  twelve 
inch  CRT’s. 

All  output  generated  by  both  correlation  programs  is 
displayed  on  the  plot  and  includes  such  parameters  as: 

-  standard  deviation  of  x  and  y 

-  correlation  and  convenience 

-  the  first  degree  best  fit  equations 
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4.6.5  Programs  CORRSNOWRATEDATA/CORRSNOLOGSCALE  (cont'd) 

The  correlation  is  produced  on  a  point  by  point  compar¬ 
ison  of  the  two  data  sets,  and  therefore  identical  times  must 
be  present.  If  a  time  discrepancy  is  found,  "INCONSISTENT 
CLOCK  TIMES"  is  printed.  The  data  files  could  be  edited 
using  "NETED"  as  required. 


4. 6. 5.1  Programs  CORRSNOWRATEDATA/CORRSNOLOGSCALE  Operating  Instruc¬ 
tions  and  Sample  Output 


ATTACH ,  I , _ * ,  ID=ROBERTSK . 

FTN, 1=1 . 

ATTACH ,  TAPEl , af gldata ,  ID= _ . 

ATTACH,  TAPE  2,  crreldata  ,ID= _ . 

ATTACH, TEK,TEXSIM,CY=2**. 

ATTACH , LIB  f  LIBNEWX59  51 , ID=WARD , MR=1 . 

LIBRARY, TEK, LIB. 

LGO 

*  EITHER  CORRSNOWRATEDATA  OR  CORRSNOLOG SCALE 

**  CY=Z  FOR  ALL  TEKTRONIX  EXCEPT  SYNCRONOUS  TERMINAL  DOWNSTAIRS 
AT  THE  CDC  CENTER,  USE  CY=1. 

TAPE1=INPUT=AFGL  DATA,  FORMAT# 3  (see  Appendix  #1) 
TAPE2=INPUT=CRRL  DATA  (see  Appendix  #2) 


4.6.6  Program  PLOTASCME 


During  the  initial  processing  of  the  CRREL  ASCME  data , 
LYC  scientists  require  time  VS  ASCME  plots.  Therefore  an 
INTERCOM  interactive  plotting  program,  PLOTASCME,  was  written 
to  produce  three  ASCME  plots  on  the  Tektronix  CRT  terminal. 


4.6.6. 1 


Program  PLOTASCME  Operating  Instructions  and  Sample  Output 


ATTACH , I , PLOTASCME , ID=ROBERTSK . 
FTN, 1=1 . 

ATTACH , TAPE 1 ,  , I D= 

ATTACH, TEK , TEKSIM, CY=2 . 

LIBRARY, TEK. 

LGO 


TAPE1*=INPUT=CRREL  DATA  (see  section  4. 6. 3.1) 
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4.6.7  programs  SMOOTHAFGL/SMOOTHCRREL 

During  snow  post  processing  development  it  became  necessary 
to  smooth  the  AFGL  snow  rates  and  CRREL  ASCME  data  as  a  tech¬ 
nique  to  remove  inherent  noise.  To  expedite  such  activities, 
DPSI  has  written  two  programs,  both  of  which  run  interactively, 
SMOOTHAFGL  and  SMOOTHCRREL . 

Both  programs  operate  identically.  The  former  requires 
the  input  file  (TAPEl)  to  be  in  AFGL  format  #3  and  the  latter 
requires  a  CRREL  data  file  (see  section  4. 6. 3.1) 

Provisions  have  been  made  to  readily  alter  the  smoothing 
neighborhood.  In  both  programs,  variable  "SMOOTH"  will  de¬ 
fine  the  number  of  points  taken  for  each  averaging  pass. 

NOTE-  this  number  must  be  odd  for  the  program  to  function  prop- 


4. 6. 7.1  Programs  SMOOTHAFGL/SMOOTHCRREL  Operating  Instructions 


ATTACH ,  I , _ * ,  ID=ROBERTSK . 

FTN , I— I . 

ATTACH , TAPE 1 ,  , I D= 

LGO 

TAPE1=INPUT=AFGL  OR  CRREL  DATA  FILE 
TAPE2=OUTPUT=SMOOTHED  DATA  OUTPUT 

EITHER 

SMOOTHAFGL 

OR 

SMOOTHCRREL 


4.6.8  Program  FALLVEL 


31 


V 


The  end  product  of  the  LYC  snow  analysis  is  the  generation 
of  a  snow  fall  velocity  plot.  This  plot  is  calculated  by  com¬ 
paring  the  ratio  of  AFGL  snow  rates  to  that  of  the  CRREL  ASCME 
^ CRREL ^  is  then  presented  on  a  time  axis. 

FALLVEL  runs  interactively  under  INTERCOM  and  generates 
three  frames,  one  for  each  ASCME  device.  Output  is  presented 
on  the  Tektronix  graphics  terminal. 
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4. 6. 8.1  Program  FALLVEL  Operating  Instructions  &  Sample  Output 


ATTACH , I , FALLVEL , ID=ROBERTSK . 

FTN ,1=1 . 

ATTACH , TAPEl ,  * , ID= 

ATTACH, TAPE 2, _ ,  ID= _ . 

ATTACH , TEK , TEKS IM , CY= Z . 

ATTACH , LIB , LI BNEWX 5551, I D=WARD , MR= 1 , 
LIBRARY, TEK, LIB. 

LGO 


TAPE1=INPUT=AFGL  DATA* 
TAPE  2=INPOT=CRREL  DATA 


AFGL  FORMAT  #3  (see  section  4. 6. 3.1) 
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5.  Theoretical  Cloud  Physics  Modeling 

5.1  Program  CUMOD  X2946:  One-Dimensional,  Mathematical/Physical  Model 

As  part  of  past  research  efforts,  LYC  scientists  developed 
a  mathematical-physical  model  of  the  life  cycle  and  precipi¬ 
tation  process  of  an  isolated  warm  cumulus  cloud  (Silverman 
&  Glass  1973).  The  model  integrates  the  vertical  equation  of 
motion  and  mass  continuity,  first  lawof  thermodynamics  ,  ard  the 
equations  of  continuity  of  water  vapor  and  liquid  hydrometeors. 
All  equations  are  formulated  in  one-dimensional  space. 

The  dynamic  interaction  between  the  cloud  and  the  en¬ 
vironment  is  modelled  by  two  entrainment  terms:  turbulent  en¬ 
trainment  representing  horizontal  mixing  at  the  sides  of  the 
cloud;  and  dynamic  entrainment  representing  the  systematic 
inflow  or  outflow  of  air  required  to  satisfy  mass  continuity. 
However,  the  unique  feature  of  this  model  is  the  handling  of 
the  microphysical  processes.  Starting  with  a  cloud  condensa¬ 
tion  nucleus  spectrum,  the  formation/growth  of  cloud  droplets 
and  raindrops  are  modeled  in  detail  for  67  logarithmically 
spaced  size  classes  that  include  a  range  of  particle  radii 
from  2  to  4040y. 

Ih  addition,  LYC  scientists  investigated  the  feasability 
of  reducing  the  amount  of  computer  memory  and  processing  time 
required  to  perform  these  detailed  microphysical  and  dynamic 
calculations  by  devising  a  numerical  method  of  simulating  the 
microphysical  process  by  a  set  of  parameterized  expressions. 


5.1.1  Current  Status  of  Program  CUMOD  X2946 


A  rigorous  microphysical  model  such  as  CUMOD  X2946  could 
be  applied  not  only  to  the  analysis  of  the  life  cycle  of  a 
warm  cumulus  cloud,  but  also  to  the  detailed  study  of  the 
formation ,  growth ,  and  evaporation  process  of  the  cloud  drop¬ 
let  spectra  and  their  role  in  the  precipitation  process. 

the  original  version  of  CUMOD  X2946  (cycle  22)  has 
been  modified  to  produce  cycle  25.  As  part  of  the  modeling 
experimentation  during  late  1982,  this  latest  version  of  the 
model  was  re-created  in  an  non-update  form  (CUMOD  ,  cycle  15) 
in  preparation  for  execution  on  the  CRAY-1  at  the  Air  Force 
Weapons  Laboratory,  Albequrque,  NM.  This  CRAY  compatible 
version  of  CUMOD  is  called  WARMCU,  cycle  14. 

The  command  deck  for  a  production  run  of  CUMOD  X2946 
(cycle  25) ,  along  with  input  data  and  sample  output  are  pro¬ 
vided  in  the  following  pages. 


153 


5.1.2  CUMOD  X2946  Job  Cards,  Input  Data,  and  Sample  Output 


COMMAND  DECK 

MILLR,CM150000 ,T500.  1416 

RFL, 150000. 

ATTACH , OLDPL , CUMODX2  9  4  6 , ID=ROBERTSK , MR= 1 . 
UPDATE (F) . 

FTN, I=COMPILE, PL=999999 , PD=6 , OPT . 

LDSET (PRESET=ZERO) 

LGO ( , , , OUTPUT) 

EXIT (U) . 

REWIND, TAPE 8. 

COPY (TAPE80UTPUT) 

DISPOSE, OUTPUT, *PR,C. 

7/8/9 


MILLERP 


UPDATE  CARDS  AS  NECESSARY 


FLORIOA  OCEAN  408  NUCLEI 


FOO  -  . 86561 7024533E+01 


SEND 
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71  SOUNDING  FLA  STD  MOO  BY  A/C 
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SEED  TITLE  CARO 


1 


VAPOR(G/G)  CLOUD(G/G)  RAIN(G/G)  DcwAT(G/G)  DGWAT(G/G)  TOTCN(G/G)  WlTDE(G/G)  DCVAP(G/G) 

3.8513E-03  9.1511E-06  8.7188E-31  0.  0.  9.1511E-06  0.  7.5536E-03 

BRKNO  (GM-1)  -  0.  CLOUD  TOP  -  850.  METERS 


5.2  PROGRAM  CUSLIDE 


This  program  utilizes  the  output  tape  from  CUMOD  X2946 
to  produce  CRT  plots  on  microfiche.  The  resultant  plots  are 
divided  into  two  categories:  Key  meteorological  parameters 
involved  in  cloud  development  and  cloud  droplet  spectra  analysi 
The  former  group  consists  of  summaries  of  precipitation,  LWC 
(cloud,  precipitation  and  total) ,  velocity,  relative  humidity, 
buoyancy,  and  radar  reflectivity;  while  the  latter  provides  a 
three-dimensional  display  of  the  droplet  spectra  within  the 
cloud  at  user  specified  heights  and  intervals  (IPR0G=1) . 

As  an  alternative,  a  modified  output  file  (IPR0G=2)  exists 
which  contains  contour  plots  without  symbols  for  the  cloud 
development  parameters;  altitude  axes  in  kilometers  instead 
of  meters;  plot  titles  with  larger  characters,  and  deletes 
the  printing  of  the  flight  identification  data  on  the  plots. 
This  format  was  chosen  for  usage  in  reports  and  slide  or 
overhead  presentations. 

Appropriate  command  deck  card  input  data  and  a  list  of 
plots  produced  follow. 


5.2.1  PROGRAM  CUSLIDE  Command  Deck  and  Input  Data 


COMMAND  DECK 

MILLR,CM120000 /T200/PK.  1416  MILLERP 

PAUSE  PLEASE  MOUNT  DISK  LYCP07 
MOUNT , VSN=LYCP  0  7 , SN=LYCP  0  7 . 

ATTACH , TAPE 8 , DATA, ID=MILLERP , MR=1 , SN=LYCP07 . 

ATTACH , S , CUSLIDE , ID=MI LLERP , MR= 1 . 

FTN , I=S , R=  3 , SL , PL=9  9  9  9  9  9 , OPT . 

REQUES  T , TAPE 3  9 , * Q . 

ATTACH , CRT , CRTPLOTS . 

LIBRARY, PRESET=ZERO. 

LGO. 

EXIT (U) . 

7/8/9 


INPUT  DECK 
CARD  1 

PLTARR ( 1 )  USER  SPECIFIED  ID(AlO) 

PLTARR(2)  PROB.  NO.  &  ID 

CARD  2 

IDENT  FLIGHT  ID, DATE, SELECTED  INPUT  VARIABLES ( 8A10) 

CARD  3 

IPROG , OPTS , NCYCLE, DELZ , TZERO, RIMAX , RCMAX (15,1011, 
15,3F10.0,2F5.2) 

CARD  4 

LEVTOP , LEVBOT (215) 
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5.2.2  List  of  Plots  Produced  by  Program  CUSLIDE 


1)  RAIN  INTENSITY  (MM/HR)  VS  TIME (Minutes) 

2)  CUMULATIVE  RAINOUT  (MM)  VS  TIME  (Minutes) 

3)  ALTITUDE  (METERS)  VS  TOTAL  LWC  (GM/M3)  VS  TIME  (Minutes) 

4)  ALTITUDE  (METERS)  VS  CLOUD  LWC  (GM/M3)  VS  TIME  (Minutes) 

5)  ALTITUDE  (METERS)  VS  PRECIPITATION  LWC  (GM/M3)  VS  TIME 

(Minutes) 

6)  ALTITUDE  (METERS  VS  VELOCITY  (M/SEC)  VS  TIME  (Minutes) 

7)  ALTITUDE  (METERS)  VS  RELATIVE  HUMIDITY  (%)  VS  TIME 
(Minutes) 

8)  ALTITUDE  (METERS)  VS  BUOYANCY  (M/SEC2)  VS  TIME  (Minutes) 

9)  ALTITUDE  (METERS)  VS  LOG 10  (RADAR  REFLECTIVITY)  VS 
TIME  (Minutes) 

10)  DROPLET  SPECTRA 

HEIGHT  (METERS  M.S.I.)  VS  RADIUS  (MICRONS)  VS 
TIME  (Minutes) 


5.3  CRAY-1  and  CDC  Interface  Procedures  and  Documentation 

LYC  increased  its  cloud  modeling  research  effort  by 
developing  and  applying  time-dependent ,  numerical  models 
which  integrate  the  numerous  and  complicated  microphysical 
and  dynamic  characteristics  of  the  atmosphere.  Such  computer 
programs  require  extensive  amounts  of  both  CPU  memory  and 
execution  time  in  order  to  model  these  processes.  Unfortune- 
ately,  the  two  CDC  6600  computers  available  at  AFGL  do  not 
adequately  meet  these  needs  in  a  cost  effective  manner.  In 
order  to  meet  the  new  computer  needs  of  LYC,  use  of  the  CRAY-1 
facilities  at  Kirtland  AFB  (AFWL) ,  Albequerque,  NM  was  in¬ 
itiated.  To  utilize  the  computer  systems  at  AFWL  from  AFGL, 
DPSI  became  familiar  with  the  three  distinct  communication 
options  available. 

-  Direct  link  to  CRAY 

-  IBM  to  CRAY 

-  SAMMET  (military  CDC  communication  network) 

Direct  link  to  CRAY-1  or  linking  with  IBM  are  the  simplest 
procedures.  They  require  a  terminal,  modem,  long  distance 
telephone  line  and  knowledge  of  CRAY-1  and  IBM  procedures. 

Each  of  these  coimtunication  procedures  only  allows  for  the 
execution/editing  of  programs  and  viewing  (via  CRT  terminal) 
of  output,  and  neither  provides  for  the  transferral  of  pro¬ 
grams  or  data  files  between  the  two  sites.  Another  drawback 
of  these  options  is  the  fact  that  only  one  user  may  access  the 
CRAY-1  or  IBM  computer  at  a  given  time. 


5.3  CRAY-1  and  CDC  Interface  Procedures  and  Documentation  (cont'd) 

A  second  alternative,  SAMNET,  will  allow  multiple  users 
access  to  the  AFWL  CDC  6600  computer.  In  order  to  use 
SAMNET,  a  variety  of  commands  and  procedure  must  first  be 
understood.  DPSI  has  been  extensively  involved  in  their  ident¬ 
ification  and  implementation  for  the  LYC  applications. 

DPSI  helped  AFGL  test  and  evaluate  a  variety  of  procedures 
for  interfacing  with  the  CRAY-1.  These  are  to  be  published 
as  an  AFGL  technical  memorandum. 

An  interesting  approach  to  job  execution  at  AFWL  (via 
SAMNET)  allows  users  to  batch  jobs  into  the  AFGL  6600.  The 
CDC,  in  turn,  activates  SAMNET  and  transmits  the  job  control 
cards  to  the  AFWL  6600,  from  which  control  can  be  forwarded 
to  the  CRAY.  When  the  job  terminates  on  the  destination  host, 
it  follows  the  flow  path  back  to  AFGL  and  output  is  printed 
on  ther  terminal  (Central,  AC,  AB)  specified  by  the  job  cards. 


6. 


LYC  Post  Processing:  PDP-8/E  Computer 


The  PDP-8/E  computer  at  LYC  consists  of  the  following 
equipment 

1)  32K  (MOS  memory  and  CPU) 

2)  Twin  Dectapes  (non-interruptable) 

3)  Tektronix  I/O  Screen/Keyboard 

4)  Input  or  output  %"  magnetic  tape 

5)  Dec-writer 

6)  Bootstrap  loader  in  read-only  memory  (ROM) 

7)  Centronics  printer 

8)  High  speed  paper  tape  reader 


The  TU-10  magtape  drive  is  primarily  used  as  an  output 
device  and  for  creating  backups  (Program  FUTIL)  of  Dectapes. 
Using  a  7  track  magnetic  tape  it  is  also  possible  to  transfer 
data  to  the  CDC  6600  computer.  For  example,  once  PDP-8  snow 
scale  data  has  been  cataloged  on  the  6600  (transferred  via 
tape)  it  can  then  be  sent  to  one  of  the  LYC  Tektronix  4052 
computers  via  a  telephone  modem. 


6.1  OS/8  operating  system 

OS/8  is  a  sophisticated  operating  system  designed  for 
the  PDP-8/E.  Besides  the  monitor  facilities,  OS/8  includes 
a  library  of  powerful  system  programs  which  allow  the  user  to 
develop  programs  using  Fortran  IV  and  assembly  language. 

To  bring  up  OS/8,  the  PDP8/E  must  be  bootstrapped: 

Turn  the  power  key  to  the  "on"  position.  A  DECTAPE  with 
the  OS/8  operating  system  (see  below)  must  be  mounted  on  DEC 
tape  drive  unit  0,  write  lock,  REMOTE.  The  switch  register 
must  then  be  set  to  7470g.  Lower  then  raise  the  "HALT/sinGLE 
STEP*'  keys.  Press  "LOAD  ADDR" ,  "EXT  LOAD  ADDR" ,  "CLEAR"  and 
"CONT" .  The  DECTAPE  will  spin  and  the  terminal  will  respond 
by  typing  a  monitor  dot  OS/ 8  is  running  the  keyboard 

monitor . 

Any  system  or  user  program  can  be  run  from  this  point. 

At  any  time  if  a  program  "bombs",  OS/8  can  instantly  be  re¬ 
started  by  this  boot  proceedure.  Detailed  instructions  on 
the  OS/8  operating  system  may  be  found  in  the  "OS/8  Handbook" 

A  list  of  the  DECTAPES  and  their  titles  that  are  used  on 
the  Cloud  Physics  PDP8/E  is  given  in  table  6,1.  The  contents 
of  each  tape  are  included  in  a  file  on  the  DECTAPE  and  can  be 
accessed  by  typing: 


6.1  OS/8  operating  system  (cont*d) 


TAPE  NO 


TITLE 


1 

2 

3 

A 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 


*  SYSTEM  WORK  TAPE 

*  OS/ 8  BACK  UP 

*  QWIK4  AND  PLOT 

*  MAINDECS  DIAGNOSTICS 

*  MORT ' S  TAPE 

*  BACK  UP 

*  A/C  BACK  UP 
RTX/8  DEVELOPMENT 

*  FORTRAN  IV  SYSTEM 

*  BASIC 

*  FLAP  TAPE 
RTX/8  TAPE 
LWC  TAPE 
FORTRAN  PROGRAMS 
TKPLOT 

WORKING  STORAGE  1 
WORKING  STORAGE  2 
WORKING  STORAGE  3 
WORKING  STORAGE  4 

*  BACK  UP 


*  indicates  the  tape  contains  the  OS/8  operating  system  and 
can  be  bootstrapped 


Table  6.1:  Lab  DEC tape  Listing 


FORTRAN  IV 


This  chapter  will  describe  the  FORTRAN  IV  system  avail¬ 
able  on  the  LYC  PDP-B/E  computer .  The  steps  necessary  to 
create  and  run  a  FORTRAN  program  are  described  in  the  next 
few  sections,  however  it  is  assumed  the  reader  has  prior 
knowledge  of  FORTRAN  programming. 

OS/8  FORTRAN  IV  provides  full  standard  ANSI  FORTRAN  IV 
under  the  OS/8  operating  system.  The  FORTRAN  IV  package 
requires  a  hardware  environment  consisting  of  a  PDP-8/E  with 
32k  of  memory,  a  console  terminal  and  a  mass  storage  device. 
The  system  is  designed  to  employ  a  KE8-E  Extended  Arithmetic 
Element,  FPP-8/E  Floating  Point  Processor,  and  any  bulk 
storage  or  peripheral  I/O  device. 

The  FORTRAN  system  is  highly  optimized  with  respect  to 
memory  requirements,  and  an  overlay  feature  is  included  that 
can  permit  programs  requiring  up  to  300k  of  virtual  storage. 
The  library  functions  permit  the  user  to  access  a  number  a 
laboratory  peripherals,  to  evaluate  a  complete  set  of  tran¬ 
scendental  functions,  to  manipulate  alphanumeric  strings,  and 
to  output  to  the  Tektronix  graphics  terminal  or  the  calcomp 
incremental  plotter. 


6.2.1  FORTRAN  system  components 


OS/8  FORTAN  IV  is  a  system  of  five  programs,  (table  6.2) 
Each  program  is  a  necessary  part  of  the  system  and  must  be 
run  in  proper  sequential  order. 

A.  TECO  -  text  editor 

B.  F4  -  FORTRAN  IV  compiler 

1.  PAS Si 

2.  PASS 2 

3.  PASS 20 

4.  PASS 3 

C.  RALF  -  Relocatable  Assembler 

D.  LOAD  -  Relocatable  binary  loader 

E.  FRTS  -  FORTRAN  Run-time  system 

Table  6.2:  OS/8  FORTRAN  IV  system 

A  FORTRAN  IV  program  written  by  the  user  is  called  a 
source  program,  to  distinguish  it  from  the  various  object 
programs  generated  by  the  OS/8  FORTRAN  IV  system,  (see  figure 
6.1).  Source  programs  are  prepared  on-line  by  means  of  the 
system  editor,  TECO. 


components  ( con t ' d ) 


DECwrlltr 


6.2.1  FORTRAN  system  components  (cont'd) 

Once  a  source  program  has  been  prepared,  it  is  supplied 
as  input  to  the  FORTRAN  IV  compiler,  F4 ,  which  translates 
each  FORTRAN  statement  into  one  or  more  RALF  (Relocatable 
Assembly  Language, Floating  Point)  statements  and  produces  an 
output  file  containing  an  assembly  language  version  of  the 
source  program. 

Compilation  errors  are  printed  on  the  DECwriter  in  two 
letter  codes,  (see  "OS/8  Handbook,  pages  8-14  for  F4  error 
messages) .  If  a  source  listing  has  been  requested,  the  errors 
are  printed,  after  the  line  in  error  in  plain  English. 

The  RALF  assembly  language  output  produced  by  the  com¬ 
piler  must  be  assembled  by  the  RALF  assembler.  During  assembly 
each  RALF  Assembly  Language  Statement  is  translated  into  one 
or  more  instruction  for  either  the  PDP-8/E  computer  or  the 
Floating  Point  Processor  and  an  output  file  is  created  con¬ 
taining  a  relocatable  binary  version  of  the  assembly  language 
input,  (gee  figure  6,2) 

The  relocatable  binary  file  produced  by  the  RALF  assembler 
is  a  machine  language  version  of  a  single  program  or  subroutine. 
This  file,  called  a  RALF  module,  must  be  linked  with  its  main 
program  (if  it  is  a  subroutine)  and  with  any  other  subroutines 
or  functions,  including  subroutines  from  the  system  library, 
that  it  requires  in  order  to  execute.  The  OS/8  FORTRAN  IV 
loader,  LOAD,  accepts  a  list  of  RALF  module  specifications  from 
the  DECwriter  and  builds  a  loader  image  file  containing  a 
relocated  main  program  linked  to  relocated  versions  of  all  sub¬ 
routines  and  library  components  that  the  mainline  requires  to 
execute,  (see  figure  6.3) 
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6.2.1  FORTRAN  system  components  (cont'd) 

The  loader  image  file  is  an  executable  core  load,  complete 
except  for  run-time  I/O  specifications.  It  may  be  saved  on 
the  DECtapes  and  run  at  any  time.  The  loader  also  provides 
for  an  optional  core  load  map  that  indicates  memory  allocation 
of  the  individual  routines  loaded.  The  overlay  feature  of  the 
loader  permits  certain  user  defined  modules  of  a  program  to 
be  stored  in  the  loader  image  file  during  execution  and  read 
into  memory  only  as  needed,  which  effectively  provides  a  ten¬ 
fold  increase  in  maximum  program  size.  (See  figure  6.3) 

The  loader  image  file  is  read  and  executed  by  FRTS,  the 
FORTRAN  RUN-TIME  SYSTEM.  FRTS  configures  an  I/O  supervisor 
to  handle  any  FORTRAN  input  or  output  reguest.  The  run-time 
system  assigns  I/O  device  handlers  to  the  I/O  unit  numbers 
referenced  by  the  FORTRAN  program,  allocates  I/O  buffer  space, 
and  also  diagnoses  certain  types  of  run-time  errors;  Such  as, 
I/O  errors,  numeric  underflow/overflow,  hardware  malfunctions, 
etc.  Run-time  errors  are  indicated  at  the  DECwriter;  fatal 
errors  cause  the  program  to  abort,  however. 

The  system  provides  complete  error  traceback  to  identify 
the  full  sequence  of  FORTRAN  statements  that  terminated  in 
the  error  condition. 

The  compiler,  assembler,  loader  and  run-time  system  each 
accept  standard  OS/8  command  decoder  option  specifications. 
"OUTPUTl ,  OUTPUT2 ,  OUTPUT3  <  INPUTl ,  ...» 1 NPUT9/0PTI  ON  ” 

Options  are  alphnumeric  charactei  r.  preceeded  by  a  slash. 
They  can  be  anywhere  in  the  command  Dine,  but  usually  placed 
on  the  right. 


6.2.2  TECO-Text  Editor 


OS/8  TECO  is  a  powerful  text  editing  and  correcting  pro¬ 
gram  that  runs  under  the  OS/8  operating  system.  TECO  may  be 
used  to  edit  any  form  of  ASCII  text  such  as  program  listings, 
manuscripts,  correspondence  and  the  like.  Since  TECO  is  a 
character-oriented  editor  rather  than  a  line  editor,  text 
edited  with  TECO  does  not  have  line  numbers  associated  with 
it,  nor  is  it  necessary  to  replace  an  entire  line  of  text  in 
order  to  change  one  character. 


Writing  a  FORTRAN  program  begins  with  the  text  editor. 
TECO  is  called  by  typing: 

.R  TECO  4 

and  responds  with  an  asterik  ,  the  OS/8  prompt.  Editing 
with  TECO  is  actually  very  simple.  All  commands  are  one  or 
two  letters  with  or  without  arguments.  The  following  table 
briefly  illustrates  some  commonly  used  TECO  commands. 


ERdev : f i lnam . ex 


EWdev: f ilnam.ex 


Itext 


Open  for  READ  (input)  'f ilnam.ex' 
on  device  'dev' 

Open  for  WRITE  (output) ’ f ilnam.ex' 
on  device  'dev' 

Clear  the  text  buffer,  then  read 
the  next  page  from  input 

Write  the  text  buffer  to  output 
Read  the  next  page  from  input 

Move  the  character  pointer 
(cp)  to  the  beginning  of  the 
next  line 

Type  the  content  of  the  text 
buffer  from  the  cp  to  the  be¬ 
ginning  of  the  next  line. 

Insert  mode.  All  subsequent  chara¬ 
cters  are  placed  before  the  cp 
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TECO-text  editor  (cont'd) 


Stext  Search  for  string  'text', 

cp  positioned  after  search  string 

d  move  the  cp  to  the  beginning  of 

the  text  buffer 

FStextl$text2$  search  for  sting  'textl'  and 

replace  with  'text2' 

EX  write  text  buffer  to  output. 

copy  remaining  input  to  output, 
close  Input/Output,  exit  to 
OS/8  monitor. 

Table  6.3;  TECO  Command  Summary 

note:  all  commands  are  separated  by  one  'ALT  MODE’  key 

(which  echos  as  a  dollar  sign)  and  a  command  strina  is 
terminated  bv  two  consecutive  altmodes. 

For  more  information  regarding  TECO,  see  the’’0S/8  Handbook1 
pages  2-132. 


6.2.3  F4  FORTRAN  IV  compiler 

The  FORTRAN  IV  compiler  accepts  one  FORTRAN  source  language 
program  or  subroutine  as  input,  examines  each  FORTRAN  state¬ 
ment  for  validity,  and  produces  a  list  of  errors  plus  a  RALF 
assembly  language  version  of  the  source  program,  along  with 
an  optional  source  listing.  A  job  which  requires  more  than 
one  module  (i.e.  subroutine)  must  have  each  program  compiled 
separately,  then  combined  using  the  loader.  F4  terminates 
by  chaining  to  the  RALF  assembler  automatically  unless  it  was 
requested  to  exit  to  the  monitor  after  compilation.  The  com¬ 
piler  is  called  by  typing: 

.R  F4* 

F4  responds  with  the  Command  Decoder  prompt,  the  asterik  . 
The  file/option  specification  line  is  entered  as: 

dev : RALF . RL , dev : LIST . LS<  dev : FTCODE . FT/ (options ) 

terminated  with  a  carriage  return.  The  compiler  accepts 
four  options.  Any  run-time  options  recognized  by  RALF, LOAD 
or  FRTS  may  be  specified  to  the  compiler. 

OPTION  OPERATION 

A  return  to  keyboard  monitor  after  compilation. 

RALF  is  not  called 

F  produce  a  RALF  listing  on  'dev:LIST.LS ' 

N  do  not  include  error  traceback  faciltities 

(decreases  memory  requirements) 

Q  optimize  cross-statement  subscripting  during 

compilation. 

Table  6,4:  F4  Compiler  Options 
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6.2.3  F 4  FORTRAN  IV  compiler  (cont'd) 


Any  errors  detected  by  the  compiler  are  reported  on 
the  DECwriter  in  short  form  (see  "OS/8  Handbook" ,  pages  8-14) 
If  a  listing  is  requested,  the  errors  are  printed  in  English. 


RALF  -  relocatable  assembler 


The  RALF  assembler  accepts  one  RALF  assembly  language 
program  or  subroutine  as  input  and  produces  a  relocatable 
binary  file,  called  a  RALF  module,  as  output.  An  optional 
listing  of  the  assembled  input  file  is  also  available.  RALF 
terminates  an  assembly  by  exiting  to  the  keyboard  monitor 
unless  it  was  requested  to  chain  to  the  loader. 

RALF  honors  three  options: 

OPTION  OPERATION 

G  After  assembly,  chain  to  the  loader,  then  to 

the  Run-Time  system 

L  After  assembly,  chain  to  the  loader,  but  not 

the  Run-time  system 

T  If  a  listing  file  has  been  specified  in  the 

command  line,  surpress  the  listing  and  produce 
only  the  symbol  table 

Table  6.5:  RALF  Options 

RALF  is  called  by  typing: 

.R  RALF 4 

and  responds  with  an  asterik.  The  command  line  appears  as: 
dev : RALF. RL, dev: LIST. LS< dev :RALF1.RA, . . . , dev: RALF 9 .RA/ (options) 


6.2.4  RALF  -  relocatable  assembler  (cont'd) 


If  more  than  one  Input  is  specified,  they  are  combined 
into  one  module  and  assembled  as  if  it  were  one  file. 

When  an  error  is  detected,  the  error  code  (see  "OS/8 
Handbook"  pages  5-38  for  RALF  error  codes)  and  the  line  in 
error  are  printed  on  the  DECwriter.  If  a  listing  was  re¬ 
quested,  the  error  code  appears  before  the  line  in  error. 

For  more  details  on  the  RALF  assembler  and  how  to  code 
Floating  Point  programs,  see  the  "OS/8  Handbook",  chapter  5 
and  pages  8-15  to  8-20. 


I 
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6.2.5  LOAD  -  relocatable  loader 

The  OS/8  FORTRAN  IV  loader  accepts  up  to  128  RALF  modules 
along  with  any  necessary  library  components,  to  form  a  loader 
image  file  that  may  be  loaded  and  executed  by  the  run-time 
system.  This  is  accomplished  by  replacing  the  relative 
starting  location  of  each  section  (module)  with  an  absolute 
core  address.  Absolute  addresses  are  also  assigned  to  all 
entry  points  defined  in  the  input  modules.  Once  all  RALF 
modules  and  library  components  have  been  assigned  to  some 
portion  of  memory  and  linked,  absolute  addresses  are  assigned 
to  the  relocatable  binary  text  and  to  the  externals. 

The  overlay  feature  of  the  loader  facilitates  running 
programs  which  are  too  large  to  fit  into  memory,  allowing 
programs  which  require  300K  of  memory  to  run  in  less  than 
32K  actual  core  memory. 

An  overlay  is  a  set  of  subroutines  stored  on  a  bulk 
storage  device.  When  any  subroutine  in  an  overlay  is  called 
by  the  mainline  or  other  subroutine,  the  entire  overlay  is 
read  into  memory,  where  it  generally  replaces  another  overlay 
of  equivalent  size. 

Overlays  are  variable-size  portions  of  memory  reserved 
for  specific  sets  of  overlays.  FORTRAN  IV  permits  eight 
levels,  0-7.  Level  0  is  always  present  in  memory,  and  con¬ 
tains  the  mainline,  common  blocks,  PDP8  mode  code,  and  library 
components . 

Levels  1  thru  7  each  may  contain  up  to  16  overlays,  one 
of  which  is  core  resident  at  any  given  time  during  execution. 
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6.2.5  LOAD  -  relocatable  loader  (cont'd) 

As  execution  begins,  overlay  MAIN  is  loaded  into  level  0  and 
started.  Other  overlays  are  read  into  memory  when  one  of 
their  constituent  subroutines  is  called.  No  two  overlays 
from  any  given  level  are  ever  co-resident 

To  call  the  loader,  type: 

.R  LOAD  + 

Load  will  respond  with  an  asterik.  The  command  string  appears 
as: 

dev : IMAGE .  LD , dev : MAP . LS<dev : PROG 9 . RL , . . . , dev : PR09 .RL/ (options) 

IMAGE. LD  is  the  loader  image  output  file.  MAP.LS  is  the  loader 
symbol  map  output  file.  Possible  run-time  options  a  a: 

OPTION  OPERATION 

C  .  Continue  the  current  line  of  input  on  the  next 

line.  The  command  decoder  only  permits  nine 
input  files  per  line.  This  option  circumvents 
this. 

G  Treat  the  current  line  of  input  as  the  last  and 

chain  to  the  run-time  system. 

L  Accept  the  single  input  file  on  this  line  as  an 

alternate  library,  used  in  place  of  the  system 
library. 

0  Close  the  current  level,  and  open  the  next  se¬ 

quential  level  for  input. 

S  If  a  symbol  map  has  been  requested,  include  system 

symbols. 

U  Ignore  the  rules  governing  subroutine  calls  be¬ 

tween  overlays 

TABLE  6.6:  LOAD  OPTIONS 
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6.2.6  FRTS  -  FORTRAN  run-time  system 

The  OS/8  FORTRAN  IV  run-time  system  reads,  loads,  and 
executes  a  loader  image  file  produced  by  the  loader.  It  also 
configures  a  software  I/O  interface  between  the  FORTRAN  IV 
program  and  the  OS/8  operating  system,  then  monitors  program 
execution  to  direct  I/O  processes  and  identifiescertain  types 
of  run  time  errors.  The  run-time  system  is  called  automati¬ 
cally  to  load  and  execute  the  loader  image  file  produced  by 
the  loader  whenever  the  "G"  option  is  specified  to  the  loader. 

The  run-time  system  is  able  to  accept  file  I/O  specifi¬ 
cations.  This  allows  the  user  to  write  a  source  program  which 
refers  to  an  I/O  device  as  an  integer  constant  or  variables. 
This  program  may  be  compiled, assembled  and  loaded  into  an 
image  file.  This  image  file  may  be  run  any  number  of  times^ 
each  time  specifying  different  physical  I/O  devices.  Thus 
logical  unit  #8  may  refer  to  the  DECwriter  in  one  run,  and 
the  line  printer  in  another  run. 

To  call  the  run-time  system,  type: 

.R  FRTS 4 

FRTS  calls  the  command  decoder  and  responds  with  an  asterisk. 
The  run-time  system  accepts  two  classes  of  input.  (1)  the 
load  module  to  be  executed.  (2)  Run-time  file  assignments. 

To  define  the  image  file,  type: 

♦dev: IMAGE. LD/ (options) 
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FF.TS  -  FORTRAN  run-time  system  (cont'd) 


Possible  options  to  FRTS  are: 


OPTION  OPERATION 

H  Halt  after  loading,  but  before  executing  the 

program.  Pressing  'cont'  switch  starts  the 

program . 


E  Ignore  the  following  run-time  errors 

a.  Illegal  subroutine  call 

b.  Reference  an  external  in  an  overlay 
other  than  in  the  form  'JSR  EXTERN' 

c.  Reference  to  an  undefined  symbol 


TABLE  6.7:  Run-Time  System  Options 


Once  the  image  file  has  been  defined,  FRTS  returns  with 
another  asterisk  to  accept  I/O  specifications.  Four  out  of 
nine  possible  I/O  unit  numbers  are  initially  assigned  by  FRTS. 

I/O  UNIT  INTERNAL  HANDLER 

paper  tape  reader 
paper  tape  punch 
line  printer 
DEC writer 

user  defined 


TABLE  6.8:  FORTRAN  I/O  Unit  Assignment 
To  associate  a  device  with  a  unit  number  type: 
dev:/n 

where  "n”  is  the  unit  number,  and  "dev:" 
non-directory  device  (LPT: ,MTA0: ) 


is  any  name  of  a 


FRTS  -  FORTRAN  run-time  system  (cont'd) 

To  define  a  file  structured  data  file  type: 

dev: file. ex/n  for  previously  created  files 
dev: file. ex</n  for  non-existent  files. 

In  any  case,  only  one  file  or  device  specification  is 
permitted  on  each  line,  and  no  more  than  six  directory  devices 
files  may  be  created  by  the  Fortran  program.  A  specification 
terminated  with  an  ALTMODE  starts  the  Fortran  program. 


6.3  Data  Processing  Programs 
6.3.1  Program  PLOT. PA 

Program  PLOT  was  written  to  aid  LYC  scientists  in  the 
calibration  procedure  of  aircraft  VCO  data  and  for  general 
data  plotting  where  best  fit  (1st  and  2nd  degree)  equations 
can  be  calculated.  Input  is  accepted  through  the  DEC-Writer 
terminal,  output  is  displayed  on  the  Tektronix  CRT. 

Program  PLOT  will 

(a)  plot  an  x-y  table  on  the  Tektronix  plotter;  the 
table  is  input  at  the  Decwriter.  If  desired, 
the  x  and/or  y  values  can  be  modified  by  logging 
them  in  order  to  produce  a  linear-linear,  a  linear- 
log,  a  log-linear  or  a  log-log  plot. 

(b)  generate  the  x  values  automatically  after  a  speci¬ 
fication  of  the  first  x  and  x-step  is  given 

(c)  allow  the  user  to  select  the  low  and  high  values 
of  x  and  y  to  be  used  on  the  plot 

(d)  plot  each  point  with  a  ' +  the  user  can  choose  v'hether 
or  not  to  connect  the  plotted  points  with  a  line 

(e)  label  the  plot  with  an  alphabetic  description 

(f)  plot  a  least  square  best  fit  curve  to  the  data 
(first  or  second  degree) ,  and  p^int  out  the 
fitting  function 

(g)  generate  a  table  of  deviations  of  the  least  square 
fit  to  the  original  data,  and  calculate  the  RMS 
error . 

(h)  allow  the  user  to  modify  the  data,  limits  and 
descriptions  and  replot,  with  a  new  least  square 
fit,  without  having  to  retype  the  entire  x-y 
table. 
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6.3.1  Program  PLOT. PA  (cont'd) 


In  order  to  execute  program  PLOT,  and  perform  the  many 
optional  capabilities,  the  following  step-by-step  procedure 
should  be  executed: 


1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 


10. 


11, 


12. 


Mount  Dec tape  5  on  unit  0,  write  enabled  and  remote 
Switches  set  to  7470 

Press  ADDR-load,  EXT. -addr-load.  Clear,  Cont 
Computer  responds  with  a  dot 
Type  R  PLOT 

Within  22  seconds  the  computer  responds  with 
PLOTTING  PROGRAM 
IS  X  ON  LOG  SCALE? 

If  the  x-data  is  to  be  logged  base  10,  answer  Y  otherwise  N 
IS  Y  ON  LOG  SCALE? 

If  the  y-data  is  to  be  logged  base  10,  answer  Y  otherwise  N 
IS  DELTA  X  CONSTANT? 

If  the  x-data  (before  taking  logs)  is  equally  spaced, 
answer  Y,  otherwise  N 

If  delta  x  is  constant,  the  computer  will  request  X-START 
and  DELTA-X,  as  X-START? 

Type  the  starting  value  and  return, 

DELTA  X? 

Type  the  difference  between  successive  x-values  and  return 
The  computer  responds  with 
ENTER  X, Y  TABLE 

or  ENTER  Y  TABLE  (for  constant  delta-x) 

For  each  value  input  the  computer  will  first  print  the 
index  number  (1,2,3, ,. .etc)  When  only  y-values  are  being 
input ,  the  user  will  type  the  y-value  followed  by  a 
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€.3.1  Program  PLOT. PA  (cont'd) 


return;  if  both  x  and  y  values  are  being  input,  the 
user  will  type  the  x-value,  followed  by  a  comma,  followed 
by  the  y-value  and  a  return. 

13.  When  all  values  have  been  inputted,  type 

1E35  return 

14.  The  computer  will  respond  with 

LIMITS  OF  X  ARE  xxxxx,  xxxxx 
which  indicate  the  low  and  high  values  of  the  x- table. 
Then  it  prints 

TYPE  LIMITS  TO  USE 

The  computer  expects  two  numbers,  separated  by  a  comma, 
the  first  number  is  the  value  of  x  at  the  left  end  of 
the  plot;  the  second  is  the  value  of  x  at  the  right  end 
of  the  plot.  There  will  automatically  be  10  divisions 
in  this  range.  If  the  first  limit  typed  is  unacceptable 
(larger  than  the  limit  found)  or  if  the  second  limit 
typed  (smaller  than  the  limit  found)  is  unacceptable, 
step  14  is  repeated. 

15.  Step  14  is  done  for  Y 

16.  The  computer  will  print 

CONNECTING  LINE? 

If  a  line  is  desired  between  points  in  the  order  they 
were  typed  in,  answer  Y.  An  answer  of  N  will  eliminate 
the  connecting  line. 

17.  The  computer  will  respond 

TURN  ON  PLOTTER,  TYPE  DESCRIPTION 
Type  a  description  to  be  printed  on  the  bottom  of  the 
plot.  This  description  is  limited  to  67  characters.  A 


6.3.1  Program  PLOT. PA  (cont'd) 

should  type  the  entire  description  over  again.  Make 
sure  plotter  is  on,  and  "on-line"  before  pressing  return. 

18.  The  plot  will  be  made.  At  its  conclusion,  the  computer 
will  ask 

LEAST  SQUARE  FIT? 

If  a  curve  fitting  the  data  is  desired  answer  Y;  other¬ 
wise  answer  N  and  proceed  to  step  22, 

19.  The  computer  responds 

DEGREE  (1  OR  2)? 

Degree  1  is  a  straight  line  fit;  degree  2  is  a  quadratic 
fit.  Respond  either  1  or  2,  Any  other  response  will 
force  step  19  again. 

20.  The  fitting  curve  will  be  plotted.  The  least  square  poly 
nominal  will  be  printed 

(a)  at  the  top  of  the  plot,  and 

(b)  on  the  Decwriter 

21.  The  computer  then  asks 

DEVIATIONS? 

If  the  table  of  x,  y,  calculated-y,  deviation-y  and  RMS 
is  desired,  answer  Y;  otherwise  answer  N 

22.  The  computer  will  ask 

REPLOT? 

If  the  data  is  to  be  modified,  eliminated,  or  extended, 
or  if  another  plot  is  to  be  made  changing  either  x  or  y 
limits,  or  the  description,  or  the  fitting  function, 
respond  with  Y;  otherwise  respond  N  and  proceed  at  step  6 

23.  The  computer  will  ask 

NUMBER? 

If  there  are  no  changes  go  to  step  23(c) . 
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6.3.1  Program  PLOT. PA  (cont'd) 


(a)  The  computer  expects  the  index  number  of  the  data 
point  being  changed.  It  will  accept  an  index  num¬ 
ber  one  greater  than  the  length  of  the  table  (omit 
the  number  corresponding  to  1E35)  with  the  assump¬ 
tion  that  the  table  is  being  extended.  The  table 
may  be  extended  repeatedly  by  entering  an  index 
number  equal  to  one  more  than  the  previous  maxi¬ 
mum  number.  If  the  number  typed  is  invalid,  step 
23  is  repeated.  After  the  number  is  accepted,  the 
computer  will  type  X,Y= 

Respond  with  the  value  of  x,  a  comma,  the  value  of 
y  and  return.  Step  23  will  be  repeated.  The  data 
being  typed  will  be  logged  according  to  responses 
made  in  steps  7,8 

CAUTION:  DO  NOT  ENTER  1E35  FOR  THE  LAST  NUMBER. 

The  program  knows  how  many  points  there  are  in  the 
table. 

(b)  If  a  value  of  x,y  is  to  be  eliminated,  after  its 
index  number  has  been  accepted,  and  X,Y  =  has  printed 
type 

1E35,1E35  return 

(c)  VThen  all  changes  have  been  made,  answer  NUMBER?  with 
0  return 


The  procedure  restarts  at  step  14 


ERRORS 


Alphabetic  errors;  If  an  error  is  discovered  before  re¬ 
turn  is  pressed,  and  the  information  was  alphabetic,  as 
in  steps  7,8,9,16,17,18,19,21,22,  type  the  rubout  key; 
the  computer  will  respond  with  a  carriage  return;  retype 
the  entire  response,  beginning  with  the  first  character. 


6.3.1  Program  PLOT . PA  ( cont ' d ) 

2.  Numeric  errors:  If  an  error  is  discovered  before  a  comma 
or  return  is  pressed,  and  the  information  was  a  number, 
as  in  steps  10,12,14  and  23,  type  the  rubout  key;  the 
computer  will  not  respond  (automatic  feature  of  the  27 
bit  floating  point  package) ;  then  type  the  entire  number 
over  again,  beginning  with  the  first  character.  CAUTION: 
If  two  numbers  were  being  inputted,  separated  by  a  comma, 
the  rubout  will  erase  the  CURRENT  NUMBER  ONLY. 

3.  Incorrect  responses  to  Y  or  N:  The  computer  will  specie 
ficially  look  for  the  Y.  Any  character  other  than  Y  will 
be  treated  as  N  except  a  carriage  return,  which  should 
be  used  only  after  a  character  has  been  typed. 

RESTART  ADDRESS:  0200 


6.3.2  Program  KNOLlD . FT 


a 


Introduction 

KNOLlD  is  a  PDP-8/E  Fortran  IV  program  which  has  been 
developed  for  the  Cloud  Physics  Branch  to  obtain  quick  and 
informative  printer  display  of  aircraft  data  from  the  lab 
PDP-8/E  computer. 

KNOLlD  will  accept  input  from  either  the  PMS1D  Kennedy 
recorder  or  the  PDP-8/E  TU-10  flight  tape.  The  output  device 
is  the  line  printer  and  an  example  may  be  found  in  figure  6.3 
Note  however,  the  example  contains  one  'block'  of  data,  where 
as,  KNOLlD  prints  two  per  page. 

Operation 

The  KNOLlD  program  is  currently  stored  on  the  "Fortran 
System"  DECtape  (#9) ,  the  source  code  is  on  tape  14,  "For¬ 
tran  programs".  To  run  KNOLlD,  tape  #9  must  be  mounted  on 
unit  0  and  the  operating  system  booted  up.  The  TU-10  or 
Kennedy  tape  should  be  mounted  on  the  magtape  drive,  set  to 
unit  0,  write  protect,  and  on-line.  In  response  to  the  OS/8 
monitor  dot  on  the  DECwriter,  type  'EXECUTE  KNOLlD. LD.' 

FRTS  will  load  and  start  running  KNOLlD. 

KNOLlD  begins  by  asking  a  battery  of  questions: 

"PDP-8/E  KNOLlD  VX.YY" 

1  "FLIGHT  NUMBER?" 

THE  FLIGHT  NUMBER  IS  ENTERED  IN  THE  FORM  "EYY-NN" 


6.3.2 


KNOLlD.FT  (cont'd) 

2  "TAPE  FORMAT? 

0  -  PMS 

1  -  TU-10  (RAWINS) 

2  -  TU-10  (FPPINS) 

3  -  TU-10  (SEMI-INS)" 

THE  TAPE  FORMAT  MAY  BE  PMS  OR  TU-10  (SELF  EXPLANATORY)  , 
RAWINS  IS  FOR  TAPES  WITH  INS  FORMATTED  LAT/LON.  FPPINS 
IS  DATA  CONVERTED  TO  FORTRAN  FLOATING-POINT,  SEMI -INS 
IS  THE  CURRENT  FORMAT  AND  CONTAINS  LAT/LON  IN  BINARY 
FORMAT-ADJUST  FORM. 

3  "A-D  BUFFER? 

0  »  NEW  (>  APRIL  80) 

1  *  OLD  (<  APRIL  80)* 

IN  APRIL  1980  THE  A-D  BUFFER  FORMAT  WAS  CHANGED.  FLIGHTS 
BEFORE  APRIL  1980  ARE  «0,  OTHERWISE  #1. 

4  "FLIGHT  DATE?" 

THE  FLIGHT  DATE  IS  ENTERED  IN  THE  FORM  "DD-MMM-YY" 

5  "MAG  HEADING  DEVIATION  -  0.0,  CHANGE?” 

THE  MAGNETIC  DEVIATION  IS  PRESET  TO'  ZERO  AND  MAY  BE 
MANUALLY  CHANGED  TO  ACCOMODATE  FLIGHTS  IN  AN  INFINITE 
RANGE  OF  LATITUDES  AND  LONGITUDES.  AN  ANSWER  OF  "Y" 
RESPONDS  WITH: 

6  "CHANGE  TO?" 

THE  USER  NOW  ENTERS  THE  NEW  MAGNETIC  DEVIATION  IN 
DEGREES . 


6.3.2  KN0L1D.FT  (cont'd) 

7  "AVERAGING  INTERVAL  IN  SECONDS?" 

THE  NUMBER  OF  SECONDS  TO  AVERAGE  PER  BLOCK  OF  DATA 
(i.e.  60  SECOND  AVERAGE). 

8  "START, STOP  TIMES  (AS  'HH  MM  SS  HH  MM  SS')?" 

ENTER  START  AND  STOP  TIMES  AS  INDICATED.  TO  PROCESS  AN 
ENTIRE  TAPE,  TYPE:  M00  00  00  99  99  99" 

9  "PARTICLE  TYPE:  HARDWARE  INPUT  OK?" 

IF  HARDWARE  INPUT  OF  PARTICLE  TYPE  (ONLY  TU-10  FLIGHT 
TAPES)  FROM  THE  Cl 30  EVENT  SWITCHES  IS  DESIRED,  TYPE 
"Y" ,  OTHERWISE  MANUAL  RESPONSE  IS  NECESSARY,  AND  ... 

10  "PARTICIPATE  TYPE  (1 , 2 , 3 , 4 , 5) ?“ 

ENTER  THE  NUMBER  1-5  CORRESPONDING  THE  APPROPRIATE 
WEATHER  CONDITIONS. 

1  =  RAIN  (DEFAULT  FOR  PMS  TAPES) 

2  =  WET  SNOW 

3  =  SMALL  SNOW 

4  =  LARGE  SNOW 

5  =  BULLET-ROSETTES 

11  "COMMENTS?" 

TWENTY-FOUR  CHARACTERS  ARE  RESERVED  FOR  ADDITIONAL 
USER  COMMENTS 

TAPE  errors  are  reported  on  the  DECwriter  terminal, 
and  could  be:  parity  errors,  incorrect  record  length  and 
end  of  tape/file.  When  KNOL1D  has  reached  the  stop  time, 
it  will  return  to  question  8. 
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6.3.2  KN0L1D.FT  (cont'd) 


i  ® 


f 

I  . 


FLIGHT  ID 
FLIGHT  DATE 

INITIAL  AIRCRAFT  AND  ELAPSED  SECONDS 

DATA  AVERAGING  CVCLE 

COMMENTS 


<D 


INITIAL  LATITUDE  AND  LONGITUDE 


PARTICLE  TVPEj  HARDWARE  OR  USER  INPUT 


fu)  AVERAGED  VCO'S;  RAH  ANO  CALIBRATED 


SUMMED  SCATTER,  CLOUD  AND  PREC1P  PROBE  COUNTS 


F)  PROBE  CHANNEL  SIZES 

NUMBER/METERS**3  -  MM  BANDWIDTH 


SCATTER,  CLOUD  AND  PREC1P  PROBE  HOUSEKEEPING  STATUS  WORDS 

VARIOUS  SELECTED  PARAMETERS 


Figure  6.4A:  KNOL1D . FT  sample  output  section  definitions 
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Figure  6.4B:  KN0L1D.FT  sample  output 


6.3.3  Program  KNlUTL.FT 


Introduction 

DPSI  has  developed  a  PDP-8/E  quick  look  dump  of  PMS1D 
or  TU-10  flight  tapes.  Program  KNlUTL  will  read  either  tapes 
and  produce  a  raw  VCO  dump  of  user  selected  areas  contained  in 
the  Knollenberg  buffer. 

KNlUTL  accepts  user  input  from  the  DECwriter  terminal 
as  to  which  type  of  tape  to  process  and  what  area  of  the  tape 
to  list.  Currently ,  KNlUTL  will  process  either  a  PMS1D  tape 
or  a  TU-10  flight  tape  (in  any  of  three  formats:  RAWINS, 
FPPINS,  SEMINS) ,  and  dump  VCO,  scatter  probe,  cloud  probe 
or  precip  probe  data. 

Operation 

The  source  of  KNlUTL  is  stored  on  DEC tape  14;  the  image 
file  is  on  tape  9.  To  run  KNlUTL  mount  tape  9  on  unit  0, 
write  enable,  boot-up  the  OS/8  monitor  and  type  'EXECUTE 
KNlUTL. LD',  FRTS  will  load  and  start  running  KNlUTL. 

The  following  list  represents  operator  dialogue  with 
KNlUTL.  Note:  under  score  indicates  user  response. 

(1)  KNlUTL  V  2.20 

ENTER  FLIGHT  ID/DATE?  E80-38  ON  17-DEC- 80 

(2)  TAPE  FORMAT? 

0=PMS  1=TU- 10 (RAWINS)  2=TU-10 (FPPINS)  3 =TU- 10 (SEMINS) 
0 
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Program  KNlUTL.FT  (cont'd) 


(3)  A  TO  D  BUFFER  (not  asked  if  PMS  tape) 

0=OLD  1=NEW 

(4)  OPTIONS  _ ENTER  'T'  OR  *F’ 

DUMP  VCO  DATA?  T 

DUMP  SC A  DATA?  F 

DUMP  CLD  DATA?  T 

DUMP  PRE  DATA?  F 

(5)  TAPE  STARTS  AT  XX: XX: XX 

ENTER  START  AND  STOP  TIMES  (HH  MM  SS  HH  MM  SS) 

XX  XX  XX  XX  XX  XX 

NOTE:  If  a  start  time  of  *'00  00  00”  is  entered,  the 

first  time  on  the  tape  is  substituted.  If  "-1"  is 
entered  KN1UTL  will  restart  itself  (used  if  a  new 
tape  is  to  be  dumped,  with  out  reloading  KN1UTL, 
and  return  to  question  #1. 


6.3.4  Program  JWPLOT . FT 

One  of  the  liquid  water  content  measurement  systems 
used  by  Cloud  Physics  on  board  the  MC130E  aircraft  was  known 
as  the  "JW"  (Johnson  and  Williams) .  While  the  device  was  a 
good  means  of  measuring  LWC,  it  has  problems.  As  aircraft 
pitch  changed  during  climbs  and  decents,  it  produced 
values  below  zero,  DPSI  has  developed  a  program  (JWPLOT) , 
run  on  the  PDP-8/E  computer,  which  will  correct  JW  data  in 
real-time  on  the  Tektronix  CRT.  Running  JWPLOT  is  3  little  dif¬ 
ferent  than  other  PDP-8  Fortran  programs  because  it  requires  a 
DECtape  data  file  used  as  a  temporary  scratch  space.  The  DECtape 
must  initially  be  zeroed.  This  can  be  accomplished  via  the 
PIP  program  (Peripheral  Interchange  Program).  With  tape  #9 
on  unit  0,  and  the  data  tape  on  unit  1,  write  enabled,  type: 

,R  PIP* 

*  A: </Z$  (ALTMODE  PRODUCES  DOLLAR  SIGN) 


The  tape  on  unit  1  is  now  zeroed.  Type 


.R  FRTSI 
* JWPLOT I 

*A: JWPLOT. DF</8$  (ALTMODE) 


This  generates  a  data  file  on  unit  1  called  JV7PLOT.DF, 


6.3.4  Program  JWPLOT.FT  (cont'd) 


The  PDP-8/E  JWPLOT  program  is  constructed  from  eleven 
programs . 


JWPLOT  main  code 

CORE  memory  device  buffer 

RMTAO  PMS  tape  processor 

TEKLIB  Tektronix  library 

CQNTIM  seconds  to  H:M:S  subroutine 

INPUT  accepts  start/stop  times  flight  ID/date 

DISPL  displays  plot  area 

JWCORR  JW  correction  subroutine 

SCALE  scales  x,y  coordinates  to  plot  area 

READT  returns  one  second  of  JW/ALT  data  per  call 

LINEAX  draws  x  or  y  line  axis 


JWPLOT,  CONTIM,  INPUT,  DISPL,  JWCORR,  SCALE  and  READT  are  taken 
directly  from  the  CDC  code.  RMTAO,  CORE  and  TEKLIB  are  common¬ 
ly  used  PDP-8/E  library  routines.  LINEAX  is  a  newer  version 
of  the  line  axis  routine  previously  used  on  the  PDP-8/E. 

Figure  6.4  contains  terminal  dialogue  between  the  operator 
and  JEWPLOT.  Figures  6.5A  and  6.5B  are  sample  plot  output 
of  uncorrected  and  corrected  data,  respectively. 


CALIBRATION  OATE  IS  26-SEP-80*  CHANGE  VALUES?  <V/N> 
V 

ENTER  THE  JW-LWC  CALIBRATION  <A,  B,  C> 

“3  B8/-4N7.  4E-4«  0  0 

3.  68000E+00  7.  40GO0E-04  B.  000O0E-01 

ENTER  THE  PRESSURE  CALIBRATION  <A, 6,0 
1132.  0.  -.  0992,  1649E-6 

1.  13200E+03-9  92080E-02-1.  649OOE-07 

ENTER  THE  TAS  CALIBRATION  <A,B,  C> 

-50  0, .  05/  0  g 

i.  0000OE+01  5.  00OO0E-02  0.  00G00E-01 

ENTER  THE  IAS  CALIBRATION  <A,  B, C) 

jHS2^94i;18661-;9991E25 

£.  1.  86600E-01-9.  99100E-06 

ENTER  THE  TEMP  CALIBRATION  <A, B, C) 

-80.  74. .  0098, ,  232SE-7 

S.  07400E+01  9.  80000E-03  2.  52500E-08 

RAH  DATA  OUTPUT?  <T/F> 


USE  CALCULATED  TAS  <T/F> 

£_ 

INPUT  THE  START  AND  STOP  TIMES  <HH  MM  SS  HH  MM  SS) 
00  00  00  99  99  99 


DEFAULT  VALUES  FOR  HT<KM)  ARE  0.  0/  10.  0 
DO  VOU  HISH  TO  CHANGE  THE  LIMITS?  <V/N) 

V 

ENTER  THE  NEM  LIMITS  <MIN,  MAX) 

0.  9,  5.  0 


DEFAULT  LIMITS  FOR  JW-LWC  ARE  -0.  3,  1.  5 

DO  VOU  WISH  TO  CHANGE  THE  LIMITS?  <Y/N> 

N 

ENTER  FLT  ID  <EXX-XX> 

E8B-04 


ENTER  FLT  DATE  <DD-MMM-VV> 
'B- JAM-80  ' 


Figure  6.5:  JWPLOT.FT  Dialogue 


Program  JWPLOT.FT  (cont'd) 


6. 3.4.1  Subroutine  CONTIM.FT 

CONTIM  is  used  by  JWPLOT  and  will  convert  total  seconds 
to  "HH:MM: SS" . 

CALL  CONTIM  (SECS, TIME) 

where : 

SECS  are  total  seconds 

TIME  is  a  two  element  array  which  will  contain  "HH:MM:SS" 

NOTE:  CONTIM  will  insert  leading  zeros  for  time  parameters 
less  than  ten.  i,e.  12:03:00 

+  t 


6. 3.4. 2  Subroutine  INPUT. FT 

INPUT  is  another  JWPLOT  subroutine  which  will  accept 
from  the  user: 

(1)  start  and  stop  times  (HH  MM  SS  KH  MM  SS) 

(2)  new  height  (km)  limits 

(3)  new  JW-LWC  limits 

(4)  flight  ID 

(5)  flight  date 


This  information  is  returned  to  JWPLOT  through  common  memory. 


6. 3.4. 3  Subroutine  DISPL.FT 


DISPL  is  used  in  JWPLOT  to  draw  the  plotting  area  of 
the  JW  graph.  Flight  ID,  flight  data,  start  and  stop  times 
are  drawn  at  the  top;  x  and  y  axis  are  drawn  and  labelled. 
Whether  the  plot  is  corrected  or  uncorrected  data  is  indi¬ 
cated  at  the  bottom. 

6 . 3 . 4 . 4  Function  JWCORR . FT 

JWCORR  is  used  in  JWPLOT  to  perform  the  JW=LWC  correc¬ 
tion.  By  building  a  table  of  intercept  points  using  the 
Tektronix  cross-hairs  (identifying  true  zero  values) ,  JWCORR 
can  subtract  or  add  the  proper  offset  to  align  the  JW-LWC 
data  point  to  zero. 

6 . 3 . 4 . 5  Subroutine  SCALE . FT 

SCALE  is  used  to  program  JWPLOT  to  scale  JW/HT  data 
points  to  the  plotting  area  and  to  ensure  that  no  points 
fall  beyond  x  and  y  limits.  If  so,  the  points  are  set  to 
extreme  x  and  y  values  and  plotting  occurs  on  the  graph 
boundaries. 

6. 3. 4. 6  Subroutine  READT.FT 

READT  was  developed  from  JWPLOT  to  return  one  second  data 

«  from  a  PMS1D  flight  tape.  READT  uses  the  RMTAO  subroutine 

which  processes  four  second  tape  data;  however  for  every  call 
to  READT,  a  pointer  (KNBASE)  to  the  current  second's  data 
is  incremented  to  the  next  second.  After  four  calls  (four 
seconds  of  data) ,  RSADT  reads  another  tape  record  and  resets 


6. 3. 4. 6  Subroutine  READT.FT  (cont'd) 

KNBASE  to  second  #1.  Errors  are  directed  by  READT  and  generate 
an  error  message  on  the  DECwriter  terminal: 

"TAPE  ERROR  X" 

Where  X  is  the  error  (1-7) 

Error  codes  are  then  returned  to  JWPLOT  for  further 
processing.  (i.e.  END  OF  TAPE  CONDITION) . 

£.3.4.7  Subroutine  LINEAX.FT 

LINEAX  is  a  Fortran  subroutine  which  will  draw  and  label 
an  X  or  y  axis  on  the  Tektronix  CRT.  It  replaces  an  older 
version  also  called  LINEAX.  The  major  differences  being: 

(1)  numbers  plotted  on  the  y  axis  are  at  0°  angles  as  is 
the  x  axis,  and  (2)  axis  max  and  min  limits  are  scaled  down 
automatically  by  LINEAX  to  ensure  a  maximum  field  width  of 
4  characters  for  each  tic  mark  numeric- label . 


6.3.5  Program  SNOPLT . FT 


SNOPLT  has  been  developed  by  DPSI  for  analysis  of  snow 
*  data  collected  during  LYC's  SN0W1A  project  (1981).  SNOPLT 

contains  many  options  for  manipulating  the  data  via  a  smooth¬ 
ing  technique,  zeroing  negative  weights  or  rates,  enlarging 
and  reducing  plot  sizes,  and  specifying  the  smoothing  neighbor¬ 
hood  for  the  'DERIV'  function. 

SNOPLT  requires  as  input  a  data  file  processed  by  the 
OS/8  operating  system  (either  DECtape  or  magtape) .  Snow 
data  contained  on  this  tape  will  be  buffered  into  a  2400 
element  array  (=1:52:21  hours).  Data  is  then  processed  in 
one  hour  increments  (weight  and  rate  plots  as  defined  by  user 
input  options)  and  then  shifted  left  through  the  array  leaving 
,  300  original  data  points.  The  remaining  space  in  the  array 

(2100  points)  is  filled  with  the  tape  data  and  the  next  hour 
processed. 

I  This  method  of  shifting  the  data  left  through  the  array 

is  known  as  "the  mobile  page"  and  simulates  virtual  memory. 

This  was  necessary  in  order  to  properly  implement  the  "Deriv" 
function;  whereas  data  from  the  previous  hour  is  required  in 
I  the  N-point  neighborhood  smoothing. 

Figure  6.6  contains  the  dialogue  for  operating  SNOPLT 
and  figure  6.7  is  a  sample  plot.  A  brief  description  of  each 
I  question  (from  figure  6.6  follows: 

1)  The  frequency  of  each  snow  weight  (Delta  'T')  has  been 
1  timed  at  2.809  seconds.  This  is  a  default,  but  can  be 

f 

!  changed.  Note:  <N>  indicates  that  pressing  carriage 

return  is  the  same  as  typing  'N'  (NO). 
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Program  SNOPLT . FT  (cont'd) 


2)  If  question  one  was  yes, accept  a  new  delta  'T'. 

3)  Not  a  question;  prints  date/ID  of  collection  mission* 

4)  Allows  output  selection  of  printer,  plotter  or  both* 

5)  Page  width  is  default  to  11  inches  -  allows  redefining 
this  to  any  size* 

6)  Defines  how  many  hours  of  data  to  plot  on  each  page% 

7)  Fixed  hour  size  determined  by  questions  five  &  six.  This 
will  adjust  the  frame  size  to  center  the  number  of  hours 
within  a  page*  Adjustable  for  plotting  overlays  of  other 
cbpartment's  snow  data. 

8)  Should  a  snow  rate  plot  be  produced,  enter  T(rue)  or 
F(alse) . 

9)  Minimum  and  maximum  scale  of  snow  weight  plot. 

10)  If  snow  rate  plot  selected  (question  eight) ;  minimum  and 
maximum  scale  of  snow  rate  plot;  and  smoothing  neighborhood 
'N*  for  derivative  function. 

11)  Should  the  snow  weight  be  smoothed  prior  to  plotting; 

T(rue)  or  F(alse) 


6.3.5  Program  SNOPLT.FT  (cont'd) 


12)  If  question  eleven  is  true;  should  the  snow  weight  be 
generated  from  snow  rate  data  (True)  or  by  smooth  func 
tion  (False). 

13)  Should  negative  weights  and  negative  rates  be  set  to 
zero  (True  or  False) . 
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6.3.5  Program  SNOPLT.FT  (cont'd) 


NOTE 

1) 

2) 

3) 

4) 

5) 

6) 

7) 

8) 

9) 

10) 

11) 

12) 

13) 


User  responses  are  identified  by  underscores 


SNOPLOT  —  SNOW  RATE  PLOTTING  PROGRAM 
SNOW  WEIGHT  DELTA  ’T’  ■  2.809 
CHANGE  DELTA  'T*  (Y/N)  <N>?  _ 

Y1  ENTER  NEW  DELTA  'T'  (XX. XXX)  : 


SNOW  MISSION  ID:  S82-04  31- JAN- 8 2 

PLOTTER  OUTPUT  (T/F) : 

PRINTER  OUTPUT  (T/F):  _ 

PAGE  WIDTH  =  11.0,  CHANGE  TO  (CR*=SKIP)  :  __ 

HOW  MANY  HOURS  OF  DATA  PER  PAGE  (XX)?  _ 

ENTER  FIXED  HOUR  FRAME  SI2E  (CR=X.y) ?  ___ 

YOU  WANT  TO  INCLUDE  A  SNOW  RATE  PLOT  (T/F)? 

SNOW  WEIGHT  PLOTS 
ENTER  MIN, MAS:  _ 


Y8  SNOW  RATE  PLOTS 
ENTER  MIN, MAX: 

ENTER  NEIGHBORHOOD  F6R  DERIVATIVE  FUNCTION: 

SMOOTH  WEIGHT  DATA  (T/F) : 

Yll  SMOOTH  WEIGHT  FROM  RATE  DATA  (T/F) :  _ 

ZERO  NEGATIVE  WEIGHTS  (T/F) : 

ZERO  NEGATIVE  RATES  (T/F) : 


Figure  6,7;  SNOPLT  dialogue 


Figure  6.8:  SNOPLT  sample  plot 


Fortran  subroutines 


There  is  a  set  of  Fortran  (also  written  in  RALF  assembler) 
subroutines  used  on  the  PDP8/E  for  processing  aircraft  data. 

(1)  RMTAO :  was  written  in  RALF  assembler  to  convert  PMS  Id 

or  TU-10  flight  tapes  to  Fortran  format 

(2)  CORE:  was  received  through  the  "PDP-8  Software  News" 

Fall  1980  magazine.  CORE  enables  Fortran  IV  to 
execute  the  ENCODE  and  DECODE  functions  not  pos¬ 
sible  before  in  PDP8/E  Fortran  IV 

(3)  SCANT:  was  developed  to  eliminate  the  drudgery  of  pro¬ 

cessing  tape  errors  when  using  RMTAO.  SCANT  also 
implements  the  binary  search  method,  and  thus, 
very  quickly  searches  tapes  for  start  times. 


6. 3. 6.1  Subroutine  RMTAO.RA 


RMTAO  is  a  RALF  coded  subroutine  used  by  Fortran  programs 
to  read  PMS  ID  and  TU-10  flight  tapes.  Since  its  inception, 
many  revisions  have  been  developed  to  make  RMTAO  an  extremely 
versatile  addition  to  the  PDP8  Fortran  library.  RMTAO  features 
include: 


1)  processes  PMS  or  flight  tapes 

2)  seven  diagnostic  error  codes 

3)  five  tape  functions  (expandable 
to  eight) 

4)  TU-10  tapes  include  LAT/LON  and 
event  word  data 

5)  economical  memory  resources 


RMTAO  requires  six  arguments  per  call: 


CALL  RMTAO  (Al,A2,A3,A4rA5, A6) 


where : 


Al  =  276  element  array  tape  buffer 


Second  1 

Second  2  Knollenberg  buffer  see  Appendix  1 
Second  3 
Second  4 
Second  1 

Second  2  Latitude/Longitude  (in  radians) 
Second  3 
Second  4 
Second  1 

Second  2  Event  word  data  (1  through  5)* 
Second  3 
Second  4 


1-64 

65-128 

129-192 

193-256 

257-258 

259-260 

261-262 

263-264 

265 

266 

267 

268 


6. 3. 6,1  Subroutine  RMTAO.RA  (cont'd) 


269-270  Second  1 

271-272  Second  2  SIGN  &  ID  data  word  for  LAT/LON 

273-274  Second  3 

275-276  Second  4 


*  event  word  codes 

1  =  RAIN 

2  =  WET  SNOW 

3  =  SMALL  SNOW 

4  =  LARGE  SNOW 

5  =  BULLET  ROSETTES 


A2  =  tape  format 


0  =  Kennedy  tape 

1  =  TU-10  tape 

INS  data  in  RAW  format 

2  =  TU-10  tape 

INS  in  FPP  format 

3  =  TU-10  tape 

INS  in  semi-FPP  format 


A3  =  A-D  buffer  tape 


0  =  OLD  (BEFORE  APRIL  1980) 
1  =  NEW  (AFTER  APRIL  1980) 


A4  =  tape  error  code 


CODE  MEANING 


0  NO  ERROR 

1  PARITY 

2  END  OF  FILE 

3  INCORRECT  RECORD  LENGTH 

4  BEGINNING  OF  TAPE 

5  DRIVE  OFF  LINE 

6  END  OF  TAPE  (AT  EOT,  RMTAO  AUTOMATICALLY  REWINDS) 

7  BLANK  TAPE  READ 


6. 3. 6.1  Subroutine  RMTAO.RA  (cont'd) 


A5  =  -  number  of  second  to  process  for  READ  (FUNCTION  0) 
OR  +  number  of  records  to  space  forward  or  backward 
(FUNCTIONS  3  and  4) 

A6  =  tape  drive  function 

CODE  MEANING 

READ  A  RECORD  &  CONVERT  TO  FPP  FORMAT 
REWIND  DRIVE,  TURN  OFF  LINE 
REWIND  DRIVE 
FORWARD  SPACE 
BACKSPACE 

UNUSED 

The  source  file  of  RMTAO  is  on  tape  number  14  and 
is  called  RMTAO.RA. 


6. 3. 6. 2  Subroutine  CORE. RA 


The  CORE  subroutine,  received  from  "PDP-8  Software  News 
Fall  1980  magazine,  is  a  RALF  assembler  coded  program  which 
will  create  a  memory  I/O  buffer  accessible  via  Fortran  READ/ 
WRITE  commands.  CORE  requires  only  one  call  to  define  the 
I/O  unit  number.  More  than  one  call  will  produce  unpredic¬ 
table  results.  Calling  format  is: 

CALL  CORE  (u) 

Where:  u  is  the  I/O  unit  to  assign  the  memory  buffer  not 
previously  used  in  system  configuration  (the  LYC 
PDP8/E  system  uses  unit  5) 

A  write  to  the  memory  device  is  like  printing  on  the 
terminal.  A  read  will  transfer  the  buffer  to  the  Fortran 
program  as  though  the  characters  were  typed  on  the  keyboard. 

ENCODE  -  Simulated  feature 

A  number  can  be  written  to  the  memory  device  using  a 
numeric  format  (I,F,D)  then  read  back  using  alphanumeric 
format  (A) . 

example:  CALL  CORE (5) 

NUMBER=12 . 3 
WRITE  (5,10)  NUMBER 
10  FORMAT  (F6.2) 

READ  (5,20)  ALPHA 
20  FORMAT  (A6) 

At  this  point.  Alpha  =  "12.3" 


S.3,6.2  Subroutine  CORE .RA  (cont'd) 


The  ENCODE  feature  is  a  particularly  useful  tool  in  for 
matting  for  plotter  output. 


DECODE  -  Simulated  feature 


Virtually  identical  to  ENCODE,  DECODE  will  convert  an 
alphanumeric  string  to  numeric  format. 


CALL  CORE (5) 
ALPHA=6H  14.67 
WRITE  (5,10)  ALPHA 
10  FORMAT  (A6) 

READ  (5,20)  NUMBER 
20  FORMAT  (F6.2) 


NUMBER  is  now  set  to  14.67 


Subroutine  SCANT. RA 


Every  job  run  on  the  PDP8/E  computer  which  processes 
aircraft  data  tapes  requires  TU-10  drive  control.  Every 
job  must  also  process  tape  error  messages,  should  an  error 
occur  during  execution.  Most  jobs  handle  data  in  a  time 
frame;  that  is,  within  a  start  and  stop  time.  For  these 
reasons  DPSI  has  developed  a  RALF  assembler  routine  (SCANT) 
to  handle  magtape  overhead  control. 


Subroutine  SCANT  interfaces  the  user  Fortran  IV  program 
with  the  RMTAO  (read  magtape  0)  subroutine.  SCANT  features 
include: 

1)  reset,  read  and  scan  functions 

2)  common  memory  communication  block 

3)  stop  time  detection  and  start  time 
not  located 

4)  full  tape  error  diagnostics  in  simple 
English 

5)  aircraft  clock  rollover  algorithm 

All  communication  variables  and  parameters  are  passed 
through  common  'A'.  This  common  block  contains  the  following: 


KNBUF  Knollenberg  tape  buffer 

KNBASE  PTR  to  current  second  in  KNBUF 
FORMAT  tape  format  (PMS /TU-10) 

ATOD  A  to  D  buffer  type 

ERROR  RMTAO  error  code 

SECCNT  number  of  second  to  process  per  record  (1,2, 3,4) 
CURSEC  current  second 

RECORD  current  tape  record 

STARTT  search  start  time 


6. 3. 6. 3  Subroutine  SCANT. RA  (cont'd) 


STOPT  search  stop  tine 

FIRSTT  first  tine  on  tape 

FUNCT  SCANT  function 

TERROR  SCANT  error  code 

Internal  tine  is  always  maintained  in  a  total  seconds 
format  (H*3600  +  M*60  +  S)  as  in  variables  CURSEC,  STARTT, 
STOPT,  and  FIRSTT.  Whenever  tine  is  printed,  as  in  an  error 
message,  it  is  converted  back  to  a  normal  format  (HH:MM:SS) . 
When  a  rollover  condition  occurs,  24  hours  (86400  seconds) 
is  added  to  CURSEC,  therefore  if  a  start  time  greater  than 
23:59:59  is  specified,  it  too  should  include  a  24  hour  roll¬ 
over,  i.e.  25:05:00. 


The  three  functions  RESET,  READ  and  SCAN  of  subroutine 
SCANT  are  selected  through  an  arithmetic  variable,  FUNCT: 


FUNCT.LT. 0 


FUNCT. EQ.O 

FUNCT. GT.O 


RESET 

rewind  the  tape  drive  and  read  record  one. 
FIRSTT  is  set  to  the  first  time  on  the  tape. 

READ 

One  second  of  tape  data  is  read.* 

SCAN 

The  time  in  STARTT  is  searched  for  in  whatever 
direction  necessary  (forward  or  backward) . 


*  One  tape  record  contains  four  seconds  of  data.  When  a 
record  is  physically  read,  KNBASE  is  set  to  second  1,  (=0) ; 
when  a  READ  is  processed,  KNBASE  is  incremented  to  second  2 
(-64) ,  second  3  (-128) ,  second  4  (-192) ,  second  5;  next  record 
is  read,  and  KNBASE  is  reset  to  zero. 


6.3.6. 3  Subroutine  SCANT. RA  (cont'd) 

If  SCANT  detects  a  search  error,  it  reports  it  through 
another  arithmetic  variable,  TERROR. 

TERROR.LT. 0  the  stop  time  (in  STOPT)  has  been  exceeded  by 
the  last  read.  * 

TERROR. EQ.O  no  error 

TERROR. GT.O  the  start  time  (in  STARTT)  has  not  been  located. 

It  could  be  either  before  FIRSTT  or  after  the 
EOF,  EOT  or  BTR  (blank  tape  read) . 

*  A  stop  time  will  result  if  a)  the  stop  time  has  been  passed, 
b)  or  if  an  EOF,  EOT  or  BTR  has  been  detected.  ERROR  will 
contain  the  actual  situation.  See  section  4. 5. 5.1  for  RMTAO 
error  codes. 

After  each  call  to  subroutine  RMTAO,  SCANT  interrogates 
the  tape  error  word,  a  non-zero  value  calls  the  error  message 
processor.  There  are  seven  possible  errors: 

1)  PARITY 

displays  error  message 

"HH:MM:SS  PARITY  ERROR  ON  RECORD:  XXXX" 
the  next  record  is  read 

2)  END  OF  FILE  (EOF) 

displays  error  message 

"HH:MM:SS  END  OF  FILE  AT  RECORD:  XXXX" 

if  SCANT  was  currently  searching  for  a  start  time,  TERR0R=+1 
if  the  EOF  occurred  during  a  read,  TERROR*- l 


6. 3. 6. 3  Subroutine  SCANT. RA  (cont'd) 


3)  INCORRECT  RECORD  LENGTH  (IRC) 
displays  error  message 

"HH:MM:SS  INCORRECT  RECORD  LENGTH  ON  RECORD:  XXXX" 
the  next  record  is  read 

4)  BEGINNING  OF  TAPE  (BOT) 

no  error  message ,  used  for  search  function.  If  BOT  is 
hit  twice  in  a  row  during  a  search,  TERROR5*  +1 

5)  OFF  LINE 

"***  TU-10  OFF-LINE...  FIX  AND  PRESS  RETURN!" 

Correct  the  tape  drive  problem  and  press  the  keyboard 
return  key. 

6)  END  OF  TAPE  (EOT) 

displays  error  message 
"physical  end  of  tape" 

the  end  of  tape  marker  has  been  sensed.  If  SCANT  was 
currently  searching,  TERROR* -ft.  If  the  EOT  occurred  during 
a  read,  TERROR*- 1. 

7)  BLANK  TAPE  READ  (BTR) 

displays  error  message 
"BLANK  TAPE  READ  (HH:MM:SS)" 

The  last  read  took  too  long,  and  therefore  no  data  is  on 
the  tape.  If  SCANT  was  currently  searching  for  a  start 
time,  TERR0R®+1.  If  the  BTR  occurred  during  a  read,  TERROR 
—1. 


